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ABSTRACT 


A  quantitative  analysis  was  carried  out  on  the  performance  of 
turboprop  aircraft  within  a  microburst  windshear.  The  objective  of 
the  analysis  was  to  provide  specific  flight  procedures  for  optimal 
navigation  through  the  windshear.  The  microburst  windshear  model 
used  in  the  analysis  embodied  the  severe  characteristics  of  the 
microburst  encountered  by  Delta  Flight  191  during  an  approach  to 
landing  at  Dallas/Ft.  Worth,  2  August,  1985.  Different  escape 
strategies  were  tested  using  the  flight  performance  characteristics  of 
the  U.S.  Navy's  P-3  "Orion"  and  T-44  "Pegasus"  aircraft.  The  three 
flight  phases  investigated  were  approach  to  landing,  takeoff,  and  the 
low  altitude  ASW  mission.  Results  from  the  analysis  were  coupled 
with  the  pilot’s  view-point  from  which  conclusions  were  drawn.  The 
results  of  the  analysis  support  a  constant-pitch-angle  escape 
procedure.  The  same  procedural  steps  can  be  used  for  both  aircraft 
in  any  configuration  or  situation  with  the  difference  being  the  degree 
of  pitch  to  employ.  The  conclusions  are  in  a  format  for  integrating 
specific  microburst  escape  procedures  within  the  NATOPS  programs 
for  the  P-3  and  T-44.  i 
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I.  INTRODUCTION 


The  microburst  windshear  is  documented  as  a  serious  hazard  to 
aviation.  The  microburst  has  been  a  contributing  factor  in  several 
airline  accidents  in  the  past  two  decades.  The  seriousness  of  these 
events  has  lead  to  a  significant  amount  of  research  on  the  effect  of 
microbursts  upon  airline  transport  aircraft.  However,  little  has  been 
documented  on  the  effect  of  a  microburst  on  light-to-medlum  weight 
turbopropeller  aircraft.  The  U.S.  Navy  operates  a  fleet  of  P-3  “Orion" 
turboprop  aircraft  In  all-weather,  long-range  patrol  missions.  The 
P-3  aviation  community  requested  a  study  of  the  effects  of  a 
microburst  upon  the  aircraft's  dynamic  flight  performance  and  the 
appropriate  escape  procedure  if  one  were  encountered.  Training  of 
the  P-3  pilot  starts  in  the  T-44,  a  Beechcraft  King  Air  twin  turboprop. 
Therefore,  the  ultimate  goal  of  this  thesis  was  to  derive  viable 
microburst  windshear  escape  procedures  which  can  be  incorporated 
into  the  Naval  Air  Training  and  Operating  Procedures  Standardization 
(NATOPS)  Program  for  both  the  P-3  and  the  T-44  aircraft. 

A.  MICROBURSTS  AND  AVIATION 
1.  Mlcroburst  Definition 

The  microburst  is  a  form  of  downburst  windshear. 
Downbursts  are  associated  with  strong  convective  activity  (namely 
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thunderstorms).  Downbursts  are  formed  by  a  cold  air  mass  dropping 

from  within  a  thunderstorm  producing  strong  downdrafts.  The 

downdraft  turns  outward  as  it  approaches  the  earth.  Fujita  is  a 

prominent  pioneer  in  identifying  and  classifying  this  type  of 

windshear  [Ref.  1].  He  defines  the  microburst  as  [Ref.  2:p.  8]; 

A  small  downburst  with  its  outburst,  damaging  winds 
extending  only  4  km  (2.5  miles)  or  less.  In  spite  of  its  small 
horizontal  scale,  an  intense  microburst  could  induce  damaging 
winds  as  high  as  75  m/sec  (168  mph). 

Conclusions  drawn  from  the  Northern  Illinois  Meter ological  Research 

on  Downbursts  (NIMROD)  and  Joint  Airport  Weather  Studies  (JAWS) 

[Ref.  2]  programs  show  that  the  microburst  structure  consists  of  a 

down  flow  with  a  structure  of  one  or  more  vortex  rings  as 

illustrated  in  Figure  1.  The  microburst  may  be  wet  or  dry; 

precipitation  may  or  may  not  be  present  in  the  downdraft.  The 

design  of  a  microburst  structure  was  not  understood,  or  even 

accepted  as  a  viable  phenomenon  prior  to  1982.  Through  the  1980s, 

meterological  researchers  and  airline  accident  investigators  came  to 

the  conclusion  that  the  microburst  weather  phenomenon  is  a 

significant  hazard  to  aviation  and  has  been  a  contributing  factor  in 

many  airline  accidents.  The  accumulation  of  the  NIMROD  and  JAWS 

project  results  and  inflight  data  conclusively  shows  that  ring  vortex 

flow  is  prevalent  in  microburst  windshear  (Ref.  3]. 
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Microburst  wlndshear  is  uniquely  different  from  other  types  of 
windshear.  The  other  prevalent  type  of  windshear  hazard  to 
aviation  is  the  thunderstorm  gustfront  [Ref.  4].  A  gustfront 
windshear  is  associated  with  massive  settling  of  rain-cooled  air 
resulting  from  thunderstorm  activity.  A  sudden  wind-shift  with 
gusty  conditions  that  precedes  a  thunderstorm  is  characteristic  of  a 
gustfront.  The  main  difference  between  a  gustfront  and  a 
microburst  is  the  size  of  the  area  effected.  A  gustfront  is  widespread 


Figure  1.  Wind  streamlines  characteristic  of  a 
microburst. 

across  the  face  of  one  or  more  thunderstorms.  A  microburst  is  small 
and  Isolated  within  convective  activity.  The  local  magnitude  of  the 
horizontal  and  vertical  windshear  is  what  matters  to  an  aircraft 
penetrating  convective  activity.  Although  smaller,  a  microburst  can 


be  much  more  of  a  detriment  to  an  aircraft's  flight  performance. 


2.  The  Effect  Of  Microbursts  On  Aviation 


The  aviation  community  has  always  been  aware  of  the 
potential  danger  of  windshear.  Classification  of  different  types  was 
not  addressed  until  1976,  when  The  Federal  Aviation  Administration 
(FA A)  published  an  advisory  circular  discussing  the  effects  of 
windshear  on  flight  performance.  It  was  updated  in  1979  [Ref.  5] 
where  the  term  '‘downburst"  was  used  to  describe  a  strong 
downdraft  associated  with  the  center  of  a  thunderstorm.  There  was 
considerable  debate  on  the  hazard  of  the  downburst.  Many  argued 
that  the  downdraft  rapidly  weakened  as  it  approached  the  ground. 
The  research  by  Fujita  provided  overwhelming  evidence  to  the 
contrary,  thereby  defining  the  microburst.  By  the  mid-eighties,  the 
term  “microburst"  was  readily  being  used  by  the  aviation 
community.  However,  the  potential  effect  of  a  microburst  upon 
aircraft  performance  was  not  readily  understood  by  aircrew. 

The  danger  of  a  microburst  became  very  apparent  following 
the  crash  of  Delta  Flight  191  (DAL  191)  at  Dallas/Ft.  Worth,  August  2, 
1985.  The  crash  ensued  following  the  penetration  of  a  thunderstorm 
during  the  approach  to  landing.  The  aircraft  involved  was  a 
Lockheed  L-1011,  a  three  engine  heavy  airline  transport  aircraft.  The 
National  Transportation  Safety  Board  (NTSB)  investigating  the  DAL  191 
accident  published  in  their  report  the  following  [Ref.  6:p.  81]: 
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The  National  Transportation  Safety  Board  determines  the 
probable  causes  of  the  accident  were  the  flight  crew's  decision 
to  initiate  and  continue  the  approach  into  a  cumulonimbus 
cloud  which  they  observed  to  contain  visible  lightning;  the  lack, 
of  specific  guidelines,  procedures,  and  training  for  avoiding  and 
escaping  low-altitude  windshear;  and  the  lack  of  definitive,  real¬ 
time  windshear  hazard  information.  This  resulted  in  the 
aircraft's  encounter  at  low  altitude  with  a  microburst-induced, 
severe  windshear  from  a  rapidly  developing  thunderstorm 
located  on  the  final  approach  course. 

Several  pertinent  questions  arise  from  this  conclusion.  Airlines 

routinely  fly  through  "bad"  weather  every  day.  How  do  a  flight 

crew  and  Air  Traffic  Control  (ATC)  authorities  recognize  microburst 

type  weather?  Can  any  aircraft  navigate  successfully  through  a 

microburst  or  is  it  beyond  the  performance  capabilities  of  the 

average  airliner?  Is  there  a  significant  departure  between  normal 

piloting  techniques  and  microburst  escape  procedures?  A  large 

amount  of  research  has  been  done  to  answer  these  questions. 

However,  the  focus  has  been  directed  predominantly  toward  airline 

transport  aircraft. 

The  advent  of  Digital  Flight  Data  Recorders  (DFDft)  provided 
conclusive  evidence  of  the  structure  and  the  potential  ejects  of  a 
microburst  windshear  on  aircraft.  The  DFDR  records  the  last  30 
minutes  of  voice  transcripts,  as  well  as  several  aircraft  performance 
parameters.  Aircraft  body  axis  angles  and  accelerations,  airspeed, 
altitude,  and  time  are  a  few  of  the  recorded  parameters.  These 
records  allow  the  reconstruction  of  the  wind  field  vectors  encountered 
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by  the  aircraft  [Ref.  7].  DAL  191  was  equipped  with  a  DFDR. 
American  Airline  Flight  539  (AA  539)  followed  DAL  191  and  penetrated 
the  same  convective  weather  110  seconds  later  [Ref.  6]  as  depicted  In 
Figure  2.  This  second  aircraft  was  a  McDonnell/Douglas  MD-80,  also 
equipped  with  a  DFDR.  The  availability  of  digital  flight  data  allowed 
a  comprehensive  database  to  be  built  on  the  wind  field  generated  by 
that  microburst. 


Figure  2.  Flight  paths  of  the  two  airlines  through  the 
Dallas/Ft.  Worth  microburst. 


The  center  of  the  thunderstorm  activity  was  located 
approximately  10,300  feet  from  the  approach  end  of  Runway  17L  at 
Dallas/Ft.  Worth  Airport.  The  DFDR  data  show  that  DAL  191  initially 
encountered  a  headwind  followed  by  a  very  abrupt  tailwind  [Ref.  7]. 
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A  significant  downward  airflow  also  accompanied  the  horizontal 
windshear  which  resulted  in  a  loss  of  airspeed  accompanied  by  a 
sharp  nose  down  pitch.  A  very  high  sink  rate  resulted.  The  aircraft 
exited  the  bottom  of  the  thunderstorm  approximately  100  feet  above 
ground  level  (A6L).  Ground  impact  ensued  with  the  aircraft  in  a 
slight  nose  up  attitude  and  at  a  high  airspeed  [Ref  5]. 

In  retrospect,  several  previous  airline  accidents  can  be 
attributed  to  windshear  of  this  nature.  Fujita  had  analyzed  several 
takeoff  and  landing  airline  accidents  involving  windshear  prior  to  the 
DAL  191  crash  [Ref.  2].  However,  the  Dallas/Ft.  Worth  accident 
showed  the  need  to  address  the  microburst  hazard  and  gain  a  more 
thorough  understanding  of  how  to  alleviate  the  danger. 

3.  Reducing  The  Hazard  To  Aviation 

One  approach  to  reducing  the  hazard  of  a  microburst  on 
aviation  is  recognition  and  avoidance.  This  solution  is  very  difficult 
given  available  aviation  weather  avoidance  equipment.  The  primary 
windshear  equipment  available  today  is  the  aircraft  weather  radar 
and  ground  based  Low  Level  Windshear  Alerting  System  (LLWAS). 
On-board  weather  radar  provides  a  picture  of  the  amount  of 
convective  activity  ahead,  but  cannot  recognize  microburst  activity. 
Ground  based  LLWAS  provides  warning  when  wind  speed  sensors 
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located  about  the  airfield  measure  significantly  different  wind 
vectors.  However,  the  microburst  must  occur  at  the  location  of  the 
LLWAS  system.  An  LLWAS  system  was  installed  at  the 
Dallas/Ft.  Worth  airport  at  the  time  of  the  DAL  191  accident.  The 
microburst  occurred  on  the  final  landing  approach  corridor  two  miles 
from  the  center  of  the  airport.  Because  of  the  location  of  the  LLWAS 
sensors  around  the  airport  perimeter,  the  system  did  not  give  any 
wlndshear  warning  prior  to  the  accident.  An  LLWAS  alert  was 
sounded  as  the  thunderstorm  passed  overhead  the  airfield  several 
seconds  after  DAL  191  crashed.  Neither  weather  radar  nor  LLWAS 
provide  real-time  information  on  microburst  wlndshear. 

Very  recent  technology  has  produced  some  viable,  and 
expensive,  means  to  recognize  and  avoid  related  microburst  weather 
activity.  Ground  based  Doppler  Radar  has  been  successfully  proven 
at  Denver’s  Stapleton  airport  [Ref.  8].  On-board  forward  looking 
devices  such  as  Infrared  sensors,  lasers,  and  doppler  radar  are  being 
developed  to  provide  several  seconds  of  warning  prior  to  a 
microburst  encounter.  But  until  avoidance  systems  are  made  readily 
available,  survival  of  a  microburst  encounter  will  be  dependent  on 
the  use  of  a  successful  escape  procedure  by  the  flight  crew. 
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The  second  approach  in  reducing  the  microburst  hazard  is 
developing  viable  escape  procedures.  Research  has  shown  that  quick 
recognition  of  a  microburst  encounter  is  paramount.  The  reaction 
time  of  a  flight  crew  coupled  with  the  type  of  escape  maneuver  will 
govern  the  success  of  escape  [Ref.  9).  The  escape  procedure  the  flight 
crew  executes  depends  on  the  flight  phase  and  configuration  of  the 
aircraft.  Much  research  has  been  performed  concerning  aircraft 
performance  and  escape  procedures  with  the  focus  on  airline 
transport  aircraft  operations  to  be  described  below. 

B.  MICROBURST  ESCAPE  PROCEDURES 

A  large  amount  of  research  has  been  done  on  defining  the  optimal 
escape  maneuvers  for  airline  transport  aircraft.  Much  of  that 
research  is  based  upon  control  theory.  Comparisons  of  different  flight 
strategies  for  microburst  encounters  have  been  addressed  [Ref.  10, 
Ref.  11].  The  total  energy  concept  is  readily  used  to  compare  the 
outcomes  of  different  maneuvers  [Ref.  12,  Ref.  13].  One  very 
important  aspect  in  determining  the  optimal  escape  procedure  is  the 
ability  of  a  flight  crew  to  execute  the  appropriate  maneuver  utilizing 
the  available  flight  performance  information  available. 
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The  FAA  has  generated  an  exhaustive  Wlndshear  Training  Aid 
[Ref.  14]  aimed  at  modern  day  transport  aircraft.  It  addresses  the 
crew  training  requirements  and  suggests  viable  microburst  escape 
strategies.  Analytical  research,  as  well  as  tried  and  tested 
procedures,  have  verified  the  suggestions  presented  by  the  FAA.  All 
the  major  airlines  have  Incorporated  microburst  recognition  and 
escape  into  their  recurrent  training.  The  majority  of  the  civil  flight 
simulators  possess  some  degree  of  windshear  emulating  a  microburst. 
As  of  today,  the  aviation  community  recognizes  the  microburst  as  a 
aviation  hazard  and  has  taken  significant  steps  in  reducing  the 
vulnerability. 

The  different  escape  strategies  fall  within  three  general  categories: 

1)  Optimum  aircraft  performance  through  the  airmass. 

2)  Optimum  aircraft  performance  related  to  an  inertial  reference. 

3)  A  combination  of  the  two. 

Optimum  aircraft  performance  is  defined  as  the  performance 
providing  the  largest  cushion  of  escape  from  a  microburst  encounter. 
The  parameters  measured  for  evaluating  an  escape  strategy's 
outcome,  are  airspeed  and  altitude  gained  (or  lost). 

The  first  category  Includes  maneuvers  that  trade  altitude  for 
airspeed.  The  second  category  uses  measured  data  such  as  body  axis 
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acceleration  and  wlndshear  magnitude  in  calculating  the  optimal 
flight  path  [Ref.  12,  Ref.  10,  and  Ref.  11].  Obviously  the  second 
category  requires  significant  airborne  computing  capabilities  not 
normally  available.  The  FAA’s  Windshear  Training  Aid  recommends 
an  escape  strategy  that  falls  within  the  third  category.  The 
recommendation  is  to  maintain  a  constant  pitch  attitude  with 
maximum  engine  thrust  applied.  In  general,  increasing  pitch 
attitude  toward  z.  15o  pitch  angle  shows  optimal  performance  on  most 
of  the  large  airline  transport  aircraft.  The  constant  pitch  attitude  is 
maintained  with  disregard  to  airspeed. 


C.  EFFECTS  OF  A  MICROBURST  ON  FUGHT  INSTRUMENTS 

How  a  flight  crew  perceives  and  navigates  through  a  microburst 
encounter  Is  strongly  Influenced  by  their  flight  Instrument 
Indications.  The  primary  flight  Instruments  available  to  pilots  are: 

1)  Attitude  Indicator  (AI)  -  a  gyroscoplc/inertlal  stabilized  aircraft 
attitude  indicator.  Modern  design  allows  accurate  attitude 
Information  In  the  highest  levels  of  turbulence. 

2)  Airspeed  Indicator  -  a  pltot/statlc  system  Instrument  which 
displays  airspeed  In  knots.  It  Is  sensitive  to  airmass  pressure 
changes  and  position  error. 

3)  Pressure  Altimeter  -  a  static  system  instrument  which  relates 
altitude,  in  feet,  to  static  pressure.  It  is  sensitive  to  airmass 
pressure  changes. 
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A)  Vertical  Speed  Indicator  (VSJ)  -  a  static  system  Instrument 
which  measures  the  change  in  static  pressure  in  feet  per 
minute.  It  is  unreliable  in  turbulence  due  to  mechanical  and 
pressure  lag. 

All  aircraft  have  these  primary  instruments.  Inertial  navigation  and 
true  airspeed  computers  do  provide  other  flight  performance 
references  usually  through  flight  directors.  However,  the  majority  of 
the  flight  stations  rely  primarily  on  the  basic  flight  instruments. 

Microbursts  are  always  associated  with  some  degree  of  turbulence. 
Also,  pressure  differentials  can  be  expected  within  convective 
activity.  Ground  and  airborne  data  show  that  the  pressure  within 
a  microburst  varies  about  ±  3  millibar  [Ref.  2,  Ref.  3].  This  pressure 
change  equates  to  about  ±  80  feet  in  altimeter  variation. 

The  attitude  indicator  will  be  the  flight  instrument  least  affected 
by  mlcroburst  atmospherics.  This  Instrument  will  always  give 
reliable  aircraft  attitude  information.  Due  to  the  minor  static 
pressure  changes,  the  airspeed  indicator  and  altimeter  can  still  be 
relied  upon  as  performance  Instruments.  The  airspeed  indicator  will 
be  one  of  the  primary  instruments  to  indicate  a  windshear 
penetration.  The  VSI  would  probably  be  erratic  and  unreliable  due 
to  turbulence  and  Instrument  lag. 
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There  exits  a  means  of  measuring  angle-of-attack  (AOA)  in  most 
larger  aircraft.  Transport  aircraft  utilize  AOA  as  a  stall  warning 
device  (stick  shaker).  Navy  aircraft  incorporate  a  means  to  read  AOA 
in  units.  AOA  measuring  devices  may  have  a  high  damping  factor  in 
their  measurement,  resulting  in  substantial  lag  during  turbulent 
conditions.  Therefore,  AOA  indication  may  not  be  a  reliable 
performance  indictor  during  a  microburst  encounter. 

S.  PROBLEM  STATEMENT 

The  P-3  has  no  means  to  recognize  and  then  avoid  microburst 
windshear.  It  must  rely  solely  on  escape.  The  published  guidelines 
for  microburst  escape  are  based  upon  airline  transport  aircraft  which 
significantly  differ  from  the  P-3.  In  this  study,  the  results  of 
different  escape  procedures  based  upon  the  available  flight 
instruments  are  compared  using  the  parameters  of  altitude,  airspeed, 
and  specific  energy.  Characterizing  the  effect  of  a  microburst 
windshear  upon  the  aircraft's  performance  followed  by  the 
implementation  of  the  appropriate  escape  procedure  form  the 
objective  for  the  research.  The  final  conclusions  must  provide  viable 
escape  procedures  that  are  based  upon  the  primary  flight 
instruments. 
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The  microburst  windshear  model  selected  must  conform  to  a 


known,  measured  phenomenon.  The  windshear  experienced  by  DAL 
19i  provided  a  suitable  database  for  emulation. 

A  series  of  equations  of  motion  for  the  aircraft  were  developed 
and  designed  to  be  controlled  through  pitch  angle  and  thrust  inputs. 
These  two  parameters  emulate  the  pilot's  available  flight  control 
inputs. 

Three  flight  phases  critical  to  a  microburst  encounter  were 
considered.  Approach  to  landing,  takeoff,  and  on-station  (P-3  only) 
phases  require  operation  low  to  the  ground  in  all-weather  conditions. 
Viable  escape  procedures  using  attitude,  airspeed,  altitude,  and  angle 
of  attack  were  analyzed  for  all  three  flight  phases. 

Three  questions  were  considered  for  each  flight  phase: 

1.  What  is  the  optimum  microburst  escape  procedure  given 
available  flight  information  to  the  pilot? 

2.  Does  the  optimum  escape  procedure  change  with  gross  weight 
or  available  thrust? 

3.  What  flight  instrument  Indications,  flight  control  feeling,  and 
dynamic  response  would  be  expected  during  the  optimum  escape 
maneuver? 

Viable  microburst  escape  procedures  were  derived  for  the  P-3  and 
T-44  aircraft. 
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11.  MATH  MODELS  AND  EXPERIMENTAL  DESIGN 


This  analysis  is  based  upon  mathematical  representations  of 
microburst  windshear  and  aircraft  performance.  Microburst 
windshear,  inertial  reference,  aircraft  performance,  and  critical  flight 
phases  form  the  major  segments  of  the  total  algorithm  suite.  Most 
of  the  mathematical  simulation  performance  was  carried  out  on  a 
Macintosh  SE®  computer  utilizing  MatLab®  software. 

A.  MICROBURST  WINDSHEAR  MODEL 

A  double  ring  vortex  was  chosen  to  emulate  microburst 
windshear.  Fujita  [Ref.  2:p.  14]  demonstrated  that  the  microburst 
was  a  settling  airmass  with  vortices  generated  near  the  earth.  This 
phenomenon  produces  varying  degrees  of  three-dimensional 
windshear  when  close  to  the  earth’s  surface.  Figure  1  showed  the 
general  concept  of  a  microburst. 

A  strong  microburst  windshear  was  experienced  at  Dallas/Ft. 
Worth  on  August  2,  1985.  A  Delta  Airlines  L-lOll  (DAL  191)  crashed 
while  on  final  approach  during  a  microburst  windshear  encounter. 
An  American  Airlines  MD-80  (AA  539)  performed  a  missed  approach 
110  seconds  following  the  DAL  191  crash.  Both  aircraft  were  equipped 
with  DFDRs  which  register  a  time  history  of  the  aircraft’s 
parameters.  Body  axis  accelerations,  velocities,  and  Euler-angle 
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values  were  used  to  calculate  the  flow-field  winds  generated  from 
the  microburst.  Wingrove  and  Bach  [Ref  7]  provided  the  analytical 
means  to  calculate  the  windshear.  From  these  data,  an  insight  was 
gained  of  the  magnitude  and  characteristics  of  the  windshear. 

A  double-vortex-ring  model  was  chosen  for  the  windshear  model. 
Schultz  [Ref.  15]  devised  a  double- vortex-ring  mathematical  model 
which  closely  approximated  the  windshear  experienced  by  the  AA  539 
flight.  The  aircraft  flew  through  the  microburst  windshear  at  2500 
feet  AGL.  The  original  algorithm  was  developed  and  tested  against 
recorded  flight  data.  The  double-vortex-ring  algorithm  was  modified 
in  the  current  study  with  the  addition  of  source  flows  to  better  fit 
low  altitude  windshear.  This  modified  windshear  model  was  fit  to 
the  DAL  191  flight  data  by  the  application  of  a  parameter  sensitivity 
scheme.  Ring  location  and  vortex  Intensity  were  found  to  differ 
from  the  A  A  539  model.  However,  a  multiple  microburst  structure 
can  be  expected  during  strong  convective  activity  [Ref.  2:p.  35]. 
MatLab*  programing  was  utilized  in  this  study  and  Is  listed  In 
Appendix  A. 

1.  Vortex  Ring  Model  Development 

The  development  of  this  vortex  ring  model  comes  directly  from 
the  work  of  Schultz  [Ref.  15j.  Source  flow  was  integrated  into  this 
algorithm  to  more  closely  match  the  DAL  191  windshear.  The  wind 
component  in  the  y  direction  was  not  considered  in  the  current 
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model.  The  y  component  Is  not  easily  modeled  by  symmetric  vortex 
rings  [Ref.  15]  and  has  no  direct  impact  on  the  longitudinal  dynamics 
of  an  aircraft.  Figure  3  Illustrates  the  geometric  aspects  of  the 
vortex  ring  model.  Imaging  of  the  vortex  rings  leads  to  in  viscid 
ground  effect  upon  the  wind  field.  The  wind  field  becomes  horizontal 


Figure  3.  Vortex  ring  mathematical  representation, 
next  to  the  ground  plane  as  seen  in  nature.  Boundary  layer  effects 
were  neglected  In  this  model.  Bowles  and  Oseguera  have  performed 
research  in  boundary  layer  effects  on  windshear  [Ref.  16].  From 
Figure  3,  the  following  relationships  can  be  seen  to  hold: 
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(2) 


(3) 


where  G,  P,  and  I  are  ground,  primary,  and  image  reference 
respectively.  The  velocity  components  for  each  vortex  were  obtained 
by  differentiating  the  vortex  ring  stream  function  4>: 


y  1  3# 

V  •  ii-— 
X  rr3z 

y  1 

r  r3z 
2  r3r 


where  r  is  the  radial  distance  from  the  z  axis: 


(4) 


r-V?7v?  (5) 

The  stream  function  for  a  vortex  ring  filament  was  obtained 
from  the  evaluation  of  an  elliptical  Integral  [Ref.  15]: 


0  788X 


a254  0. 


(6) 


where  ri  and  ri  represent  the  closest  and  farthest  distances  to  the 
point  of  interest  from  the  ring’s  filament;  k  is  the  ring  filament 
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strength;  and  A  is  a  scaling  term.  Their  algebraic  relationships  are  as 


follows: 


rj-Vz^+Cr-R)^ 


r^- Vz^  +  (r+R)^ 


Xm  ^  ^ 

r  ♦r 
2  1 


where  R  is  the  radius  of  the  filament. 


Equation  6  was  differentiated  to  obtain  expressions  for  3r  and 
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(16) 


Equations  10  and  11  were  substituted  into  Equation  4  to  get  the 
individual  velocity  components  produced  by  each  ring  at  a  given 
point  in  space  (with  respect  to  the  ground  reference). 

A  velocity  damping  factor  (0  was  calculated  for  each  ring  to 
prevent  erroneously  high  values  of  velocity  near  the  filament  cores. 
The  value  ranged  from  zero  to  one;  it  is  zero  when  the  point  of 
Interest  is  at  the  viscous  core  and  approaches  one  at  increased 


distances.  The  algebraic  relationship  for  C  Is  as  follows  [Ref.  15]; 


where  ru  Is  the  closest  distance  from  the  point  of  Interest  to  the  /th 


ring's  viscous  core.  Schultz  obtained  the  weighting  factors,  l  and  p, 
by  a  visual  comparison  of  the  vortex  ring  strength  and  through  a 
parameter  estimation  scheme  [Ref.  15).  A  total  damping  factor  (Z)  is 
obtained  for  the  point  of  Interest  from  the  product  of  the  damping 
factors  for  each  of  the  four  vortex  rings: 


i 
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(20) 


The  low  altitude  windshear  model  required  the  addition  of  two 
point  source  flows  located  10,000  feet  above  and  below  the  ground 
plane  at  the  center  of  the  microburst.  The  radial  flow  is  represented 
by  the  following  expression; 


r  2«R, 


where  A  is  the  source  strength  and  Ri  is  the  radial  distance  from  the 

point  source  to  the  point  of  interest.  The  horizontal  and  vertical 
wind  vectors  contributions  are  estimated  by: 


Vz  aVr  — 


Combining  the  upper  and  lower  (mirror  image)  sources  produce: 


1*^1  ^2, 

W,  ^2fl0,000*h  10,000- h 


In  summary,  wind  velocity  for  a  point  of  Interest  in  space 

referenced  to  earth  was  calculated  in  the  following  manner: 

1)  X  (horizontal)  and  z  (altitude)  positions,  measured  in  feet, 
are  inputs  for  Equation  4.  Calculations  are  made  for  the 
four  vortex  rings  (large  primary,  small  primary,  large 
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Image,  and  small  image).  A  horizontal  wind  velocity  (vxj) 
and  a  vertical  wind  velocity  (vzj)  are  calculated  for  each 
ring. 

2)  A  damping  factor  (^j)  was  calculated  for  each  ring  using 
Equation  19. 


3)  The  induced  velocities  for  each  ring  were  added  and 
multiplied  by  the  total  velocity  damping  factor  (Z): 
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where  Wx  and  Wh  are  windshear  values  (ft/s). 

4)  For  the  low  altitude  model,  the  source  flow  velocities  were 
added  to  the  vortex  velocities  to  give  the  final  horizontal 
and  vertical  windshear  components. 

2.  Windshear  Model  Fit 

The  unmodified  double  vortex  windshear  model  was  applied  to 
AA  539  flight  data  by  Schultz,  using  a  parameter  minimization 
scheme  to  fit  the  double-vortex-rlng  model  to  the  actual  DFDR  data. 
Table  1  lists  the  parameters  varied  and  the  final  results.  Schultz 
concluded  that  the  difference  between  the  model  results  and  recorded 
wind  data  was  an  rms  (root  mean  square)  of  16  ft/s  for  the  entire 
system.  Figure  4  shows  the  comparison  of  recorded  winds  to  the 
model  prediction. 
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The  AA  539  wlndshear  model  was  fit  in  the  current  study  for 
model  verification  only.  The  windshear  that  DAL  191  encountered 
was  the  windshear  to  emulate.  This  windshear  model  was  used  to 
TABLE  1.  2500  FEET  ALTITUDE  WINDSHEAR  MODEL  PARAMETERS 


parameter _ large  rinR  small  ring 


ring  radius,  ft 

8503.3 

1701.7 

core  radius,  ft 

20004.1 

323.9 

ring  circulation,  ft'2/s 

431968.8 

57204.9 

X  position,  ft 

0.0 

50.0 

v  position,  ft 

3350.4 

830.9 

z  position,  ft 

3400.6 

2333.6 

analyze  aircraft  performance  during  takeoff  and  landing  flight 
segments.  Therefore,  the  double-vortex-rlng  model  was  modified  by 
th  parameters  listed  in  Table  2.  The  data  from  DAL  191  showed  that 
TABLE  2.  LOW  ALTITUDE  WINDSHEAR  PARAMETERS 


parameter _ large  ring  small  ring 


ring  radius,  ft 

7000 

1300 

core  radius,  ft 

2004.1 

323.9 

ring  circulation,  ft''2/s 

431968.8 

131571.3 

X  position,  ft 

2500.0 

300.0 

v  position,  ft 

-300.0 

1.0 

z  position,  ft 

3400.6 

800.0 

X  dir.  source  strength 

1355396  ft*2/s 

z  dir.  source  strength 

3049200  ft‘2/s 

the  two  rings  were  not  concentric  about  the  center  [Ref.  3].  The 
small  vortex  ring  was  lower  to  the  ground  and  displaced  In  the  x 
direction  from  the  large  ring.  Point  source  flows  located  10,000  feet 
above  and  below  the  surface  plane  were  added  to  the  double-vortex- 
rlng  model  to  Increase  the  outflow.  A  simple  parameter  sensitivity 
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Figure  4.  Vortex  ring  model  comparison  to  recorded 
flight  data  [Ref.  3]. 
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scheme  was  then  used  to  estimate  the  parameters  that  best  fit  the 
modified  windshear  model  to  the  recorded  DAL  191  flight  data. 

The  sensitivity  scheme  (listed  in  Appendix  B)  used  an  iterative 
process  of  varying  each  parameter  individually,  then  calculating  the 
X  direction,  z  direction,  and  total  system  rms  difference  between 
model  and  recorded  winds.  The  initial  parameter  matrix  was  built 
up  as  an  8x3  matrix.  The  center  column  was  composed  of  the  "best 
guess"  values.  Columns  one  and  three  were  values  obtained  from  an 
interval  surrounding  the  initial  guess  value.  First  guess  parameters, 
including  the  corresponding  Intervals,  were  obtained  by  graphical 
comparison  of  model  and  recorded  data.  The  parameter  matrix  was 
refined  after  each  successive  run  of  the  sensitivity  scheme.  Final 
parameter  resolution  was  less  than  3*  change  for  the  total  system 
rms.  Table  2  lists  the  best-fit  parameters  to  the  DAL  191  flight 
windshear  measurements.  Figure  5  is  a  graphical  representation  of 
the  current  model's  vortex  rings  and  source  flow  point  relative 
locations.  The  wind  vectors  produced  from  the  model  are  shown  in 
Figure  6.  Figure  7  compares  the  model  windshear  with  the 
measured  wind  velocities  from  DAL  191. 

Total  system  rms  was  calculated  as  16.5  ft/s  for  the  microburst 
windshear  model.  The  final  model  shows  a  close  approximation  of  an 
actual  microburst,  serving  as  a  good  windshear  model  for  testing 
aircraft  performance  during  low-level  encounters. 
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Figure  5.  Geometric  relationship  of  the  vortex  rings  and 
point  source  within  the  microburst  model. 


26 


^//////nnnmuw  iiiiiiii 

^///////mnmwwmw 

•^/////m/nnmum 

n  1 1 1 1 1 1 1 

•0»M*s'y/i'///////tllllh\\\\\\  I 
I KWWMJ 
I M  KSWW  \ 

v^yyf///n  /Ml  I  IK^WWiij 
1 1 1  tV«XW\\ 

•  1 1 1  J 

. 

. 


.\\1111 

tuiti 


o 

o 


o 

o 

o 

<o 


o 

o 

o 


. 

. . 


. . 

~»n«a>v>N\\\x\\n  1 1 

“'wwwwwu  1 1 1 1 1  iy////m 

-••"iiiiiiiiiniiuiy/////ii 

I  1  I 

o  o  o 

o  o  o 

o  o  o 

CO  CM  »- 


o 

c 

8 

ffi 

2 

E 

o 

I 

CO 

w 


o 

o 

o 

CM 


o 

o 

o 


o 

o 

o 

(D 


m 


o 

o 

o 


o  w 
o 

CO  Q> 


8 
CD 
o  S 

O  E 

CM  § 
w 

8 

c 

_  CO 

3  To 
2  =6 


o 
o 
o 
o  r* 


(u)  epn)!)|e 


(W)  epninie 


Figure  6.  Vector  plots  ot  the  wind-fields  generated  by 
the  microburst  wlndshear  model. 
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Figure  7.  DAL  191  recorded  winds  compeured  to 
mathematical  model. 


B.  INERTIAL  REFERENCE  FRAME 

The  concept  of  aircraft  specific  energy  was  used  to  compare  the 
outcomes  of  different  windshear  escape  strategies.  The  aircraft’s 
total  energy  is  defined  as  the  sum  of  the  airmass-relative  kinetic 
energy  and  the  inertial  potential  energy.  A  vectorial  relationship  of 
an  aircraft’s  motion  through  the  airmass  is  presented  in  Figure  8.  A 

i  h 
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Figure  8.  Inertial  reference  coordinate  system, 
set  of  coupled,  non-linear  differential  equations  can  be  developed 
from  this  vectorial  relationship.  These  differential  equations  can  be 
solved  numerically  where  specific  energy  can  be  determined  at  any 
point  in  Inertial  space. 
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1.  Total  Energy  Concept 


The  aircraft's  altitude  and  airspeed  at  any  given  point  can  be 
used  to  determined  total  energy  from  the  following  relationship: 

E  =  16mV2  +  mgh  (26) 


where  m  is  mass,  V  is  airspeed,  and  h  is  the  altitude.  The  specific 
energy  (sometimes  called  energy  height),  Es,  is  defined  by; 


Es-  •^♦h 
2g 


(27) 


The  time  rate  of  change  of  specific  energy  is  equal  to  aircraft 
acceleration  plus  rate  of  climb.  Differentiating  Equation  (27)  with 
respect  to  time  gives  specific  power: 


(28) 


Values  for  V,  V,  h,  and  h  can  be  obtained  at  any  instant  of  time 
from  the  relevant  equations  of  motion. 

2.  Equations  of  Motion 

The  relationship  of  an  aircraft's  motion  through  a  moving 
airmass  relative  to  earth  leads  to  the  development  of  six  coupled, 
non-linear  equations  [Ref.  123.  Using  the  coordinate  system  in 
Figure  8,  the  following  equations  of  motion  evolve: 

x-Vcosr^+Wx  (29) 

h-Vsinyj^*Wh  (29) 
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(30) 

(31) 

(32) 
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Note  that  at  and  at  were  set  equal  to  zero  due  to  steady  state 
conditions  being  assumed  for  the  windshear  model.  This  assumption 
is  valid  due  to  the  minimal  time  the  airplane  is  exposed  to  the 
microburst  (30  to  60  sec.).  Note  that  the  equations  are  for  point 
mass  analysis.  No  dynamic  cross-coupling  is  considered. 

The  equations  of  motion  were  solved  using  a  predictor-corrector 
numerical  scheme.  Euler  First-Forward,  Euler  Half-Step  and 
Richardson  Extrapolation  [Ref.  17)  were  combined  in  the  following 
manner: 


.a) 


n 


yn*l/2-yn*fyn 


(2)  0)  At 

^n+1  ^n+I/2  2  ^n+1/2 

y  .2y^2)  (1) 


(34) 

(35) 

(36) 

(37) 


31 


Second-order  accuracy  Is  expected  from  this  numerical  scheme.  By 
providing  input  values  for  aircraft  constants  and  initial  values  for 
the  inertial  reference,  the  equations  of  motion  were  solved  to  obtain 
the  aircraft  performance. 


C.  AIRCRAFT  PERFORMANCE  MODEL 

A  MatLab®  program  was  written  to  take  aircraft  performance 
parameters  and  use  the  set  of  motion  equations  to  measure  aircraft 
response  in  the  windshear.  Initial  values  were  determined  from 
aircraft  performance  parameters  and  the  initial  position.  Iteration  of 
the  set  of  motion  equations  was  based  upon  the  time  interval 
chosen.  Values  for  x,  h,  V,  Vi  Yj,  Y»,  Es,  and  Es  were 

tabulated  after  each  iteration.  The  control  of  the  model  was  taken 
from  a  pilot's  point  of  view.  The  controlling  variables  were  0  (deck 
angle  or  pitch  angle)  and  T  (thrust).  Both  variables  can  be  changed 
between  iterations. 

Certain  aircraft  performance  parameters  were  provided  for 
calculating  lift,  drag,  and  AOA.  The  parameters  were  Clo»  Clo,  Cdo,  S, 

K,  W,  maximum  thrust  available,  and  maximum  AOA.  Cl  and  Cd 


were  calculated  as  follows: 


(38) 
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(39) 


Note  that  the  lift-curve  slope  was  referenced  from  the  longitudinal 
axis  of  the  aircraft  (see  Figure  8).  Lift  and  drag  were  calculated 
using  the  familiar  relationships; 


L-QSC, 


(40) 

'D  (41) 

where  the  dynamic  pressure,  Q,  is  based  upon  equivalent  airspeed: 


D-QSC, 


(42) 


Initial  values  for  airspeed,  thrust,  and  theta  were  chosen  to 
match  the  aircraft  requirements  for  the  phase  of  flight  under 
scrutiny.  Initial  values  for  x  and  h  were  chosen  referenced  to  the 
microburst  wlndshear  center.  The  iteration  of  the  model  was  based 
upon  the  time  step.  At. 

Three  aircraft  were  analyzed  for  their  response  in  a  microburst 
windshear.  They  were  the  U.S.  Navy  P-3  (Lockheed  L-188),  the  U.S. 
Navy  T-44  (Beechcraft  King  Air  H-90),  and  a  generic  3-engine  "heavy" 
airline  transport.  The  P-3  is  a  four-engine  turboprop  with  gross 
weights  in  the  medium  range  (75,000  to  135,000  lbs.).  The  T-44  is 
powered  by  two  turboprop  engines  and  falls  within  the  category  of  a 
“light-twin"  transport.  The  3-engine  heavy  airline  transport  is 
powered  by  turbofans  and  as  the  name  implies,  is  capable  of  high 
gross  weights.  Tables  3,  4,  and  5  delineate  the  performance 
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parameters  of  the  three  aircraft.  The  lift  curve  slope  for  the  P-3  is 
greater  than  the  theoretical  prediction.  This  result  is  due  to  the 
influence  of  the  propeller  induced  flow-field.  (Power-off  CLa=5.7  rad-i.) 


TABLE  3:  P-3  PERFORMANCE  PARAMETERS  [Ref.  18]. 


Performance 

Landing 

Approach 

Takeoff 

On-station 

Parameters 

Configuration 

Configuration 

Configuration 

W  (lbs) 

89,500;  114,000 

89,500;  120,000; 

120,000 

S  (ft"2) 

1300 

135,000 

1300 

1300 

K 

0.05041 

0.05041 

0.05041 

Cdo  4  engine 

0.0567 

0.0551 

0.0213 

3  engine 

0.0630 

- 

- 

Clc  (rad-0 

6.38 

6.38 

6.25 

Clo 

0.800 

0.800 

0.350 

a  max  (rad) 

0.244 

0.244 

0.209 

Tmaxdbf)  4  eng. 

33400 

33400 

33400 

3  eng. 

25050 

- 

- 

q  (rad/s) 

0.0873 

0.0873 

0.0873 

Vref  (ft/s) 

236;  262 

- 

354 

VI  (ft/s) 

- 

204;  214;  229 

- 

V2  (ft/s) 

220;  227;  239 

- 

As  mentioned,  the  AOA  used  in  the  calculations  was  based  upon  the 
longitudinal  axis  of  the  aircraft.  Maximum  thrust  available  and 
gross  weight  were  adjusted  to  meet  the  scenario  requirements. 
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TABLE  4.  3-ENGINE  HEAVY  AIRLINE  TRANSPORT 
PERFORMANCE  PARAMETERS  [Ref.  10]. 


Landing 

Performance 

Approach 

Takeoff 

Parameters 

Configuration 

Configuration 

W  (lbs) 

362,000 

462,000 

S  (ft'2) 

4578 

4578 

K 

0.059 

0.059 

Cdo 

0.108 

0.098 

Clo  (rad-1) 

4.96 

4.96 

Clo 

0.532 

0.532 

®  max  (r^d) 

0.317 

0.314 

Tmaxdbf) 

126,000 

126,000 

q  (rad/s) 

0.0873 

0.0873 

Vref  (ft/s) 

227 

- 

VI  (ft/s) 

- 

238 

V2  (ft/s) 

- 

255 

TABLE  5: 

T-44  PERFORMANCE  PARAMETERS  [Ref.  19]. 
Landing 

Performance 

Approach 

Takeoff 

Parameters 

Configuration 

Configuration 

W  (lbs) 

8,280 

8,280 

S  (ft"2) 

210 

210 

K 

0.0503 

0.040 

Cpo 

0.120 

0.100 

Clo  (rad-1) 

6.24 

6.24 

Clo 

0.587 

0.0523 

a  max  (rad) 

0.244 

0.227 

Tmax(lW) 

3,023 

3,023 

q  (rad/s) 

0.0873 

0.0873 

Vref  (ft/s) 

203 

- 

VI  (ft/s) 

- 

152 

V2  (ft/s) 

- 

202 
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For  each  Iteration  of  At,  the  following  aircraft,  Inertial,  and 


wlndshear  model  values  were  calculated  and  recorded: 

1)  X  -  Distance  from  microburst  center  in  feet. 

2)  h  -  Absolute  altitude  in  feet. 

3)  V  -  Equivalent  airspeed  in  knots. 

4)  e  (theta)  -  Pitch  angle  in  degrees. 

5)  a  (alpha)  -  Angle-of-attack  in  degrees. 

6)  Ya  “  Flight  path  in  degrees  referenced  to  the  airmass. 

7)  Yi  -  Flight  path  in  degrees  referenced  to  the  earth. 

8)  Vg  -  Ground  speed  in  knots. 

9)  ROC  -  Rate  of  climb  in  ft/sec.  referenced  to  the  earth. 

10)  Thrust  -  instantaneous  thrust  in  Ibf . 

11)  Es  -  Specific  energy  In  feet. 

12)  Es  (Esdot)  -  Time  rate  change  of  specific  energy  in  ft/sec. 

13)  Wx  -  Horizontal  wind  speed  In  ft/sec. 

14)  Wh  -  Vertical  wind  speed  in  ft/sec. 

From  the  collection  of  data,  plots  were  made  of  altitude,  theta,  alpha, 
airspeed,  specific  energy,  and  time  with  comparison  to  distance  from 
microburst  center. 
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Four  basic  escape  maneuvers  were  examined  for  critical  phase  of 


flight.  Constant  airspeed,  constant  altitude,  constant  pitch  angle  (6), 
and  constant  angle-of-attack  (a)  escape  maneuvers  were  evaluated 
quantitatively  (engineering  standpoint)  and  qualitatively  (piloting 
standpoint). 

D.  CRITICAL  FLIGHT  PHASES  MODELED 

There  are  three  phases  of  flight  in  which  a  microburst  encounter 
becomes  critical  for  airplane  survival.  Approach  to  landing,  initial 
takeoff  and  climb  out,  and  low-altitude  maneuvering  are  the  three 
phases  of  flight  considered  in  this  analysis.  All  three  phases  occur  at 
an  altitude  of  less  than  1000  feet  AGL.  The  low-altitude 
maneuvering  flight  phase  is  applicable  to  the  P-3  community.  The 
P-3  routinely  operates  200  feet  above  the  water  in  all  weather 
conditions  during  Anti-Submarine  Warfare  (ASW)  missions. 

1.  Landing  Approach 

The  landing  approach  scenario  is  based  upon  a  3  degree 
glideslope  descent  to  landing.  In  this  scenario,  the  microburst  center 
is  placed  10,300  feet  from  the  end  of  runway.  The  simulation  starts 
500  feet  (x=-500  feet)  before  the  microburst  center  on  glideslope 


37 


(h*566  feet).  This  situation  closely  represents  the  scenario  of  DAL  191. 
The  aircraft  is  exposed  to  the  windshear  immediately  at  start,  t=0. 
The  simulation  is  run  for  a  programmed  length  of  time  or  until 
ground  impact  occurs.  Five  aircraft/weight  combinations  were 
analyzed  for  the  approach  to  landing  scenario.  They  were  a  P-3  at 
89,500  pounds  gross  weight,  a  P-3  at  114,000  pounds  gross  weight,  a  3- 
engine  heavy  airline  transport,  a  T-44,  and  a  P-3  at  89,500  pounds 
gross  weight  with  one  engine  shut  down. 

Starting  aircraft  parameters  were  such  that  the  given  pitch 
angle  (theta)  and  thrust  would  maintain  a  30  glldeslope  at  target 
approach  airspeed  (Vref).  Pitch  angle  and  thrust  were  maintained 
until  a  loss  of  airspeed  equated  to  Vref  minus  20  knots.  At  such 
time,  the  aircraft  performed  one  of  the  following  programmed  escape 
procedures: 

1)  constant  airspeed  -  Maximum  thrust  was  applied  and  a  pitch 
angle  set  to  Oo,  This  pitch  angle  was  maintained  until  the 
airspeed  equaled  Vref.  Pitch  angle  was  then  adjusted  to 
maintain  the  airspeed  at  Vref  1  5  knots. 

2)  constant  altitude  -  Maximum  thrust  was  applied  and  a  pitch 
angle  was  set  to  obtain  a  positive  rate  of  climb.  This  pitch 
angle  was  maintained  until  the  target  altitude  (altitude  at 
which  the  maneuver  began)  was  established.  Pitch  angle  was 
then  adjusted  to  maintain  the  target  altitude  by  ±  20  feet. 
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3)  constant  theta  -  Maximum  thrust  was  applied  and  a  pitch 
angle  (theta)  was  set  and  maintained.  Specifically,  theta 
values  of  5o,  go,  lOo,  and  15®  were  used. 

4)  constant  alpha  -  Maximum  thrust  was  applied  and  a  pitch 
angle  was  adjusted  to  get  a  given  angle-of-attack  (alpha).  AOA 
values  of  12,  15,  and  20  units  were  used  for  the  P-3  model. 

The  above  escape  maneuvers  were  constrained  by  certain  limits. 

Pitch  rate  (q)  was  set  at  5o/sec.  Thrust  application  rate  for  the  P-3 

model  was  set  at  0.5  maximum  thrust/sec.  Thrust  application  rate 

for  the  generic  3-englne  heavy  and  T-44  aircraft  was  set  at  0.2 

maximum  thrust/sec.  This  rate  accounts  for  engine  "spool-up"  time. 

Maximum  pitch  angle  was  limited  not  to  exceed  maximum  AOA. 

Appendix  C  contains  the  MatLab#  program  used  for  the  approach  to 

landing  scenario. 

2.  Takeoff 

The  takeoff  scenarios  primarily  explored  the  effects  of 
microburst  wlndshear  on  the  takeoff  performance  with  penetration 
at  liftoff.  However,  some  analysis  was  performed  on  the  microburst 
center  distance  from  liftoff  point.  Generally,  the  simulation  began 
with  the  aircraft  lifting  off  1200  feet  from  the  microburst  center 
(x=-1200  feet)  at  the  appropriate  liftoff  speed  (Vi).  Initial  pitch  angle. 
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theta,  was  such  to  achieve  takeoff  safety  speed  (V2)  at  50  feet,  no 


windshear.  Maximum  thrust  was  used  for  all  cases.  Execution  of 
the  particular  escape  maneuver  was  begun  when  the  rate  of  climb 
(ROC)  is  less  than  or  equal  to  zero  or  the  airspeed  fell  below  V2.  Five 
aircraft/weight  combinations  were  looked  at.  They  were  a  P-3  at 
90,000  pounds,  114,000  pounds,  and  135,000  pounds  gross  weight,  as 
well  as  a  3-englne  heavy  airline  transport  and  a  T-44. 

Four  escape  methods  were  considered  for  a  microburst 
encounter  at  takeoff.  When  the  ROC  was  less  than  or  equal  to  zero, 
or  the  airspeed  fell  below  V2,  one  of  the  following  programs  was 
executed. 

1)  constant  airspeed  -  If  airspeed  was  less  than  V2  at  initiation, 
theta  was  reduced  to  0®.  If  airspeed  was  greater  than  V2  at 
Initiation,  theta  was  increased.  In  both  cases,  theta  was 
manipulated  to  maintain  airspeed  =  V2  ±  5  knots  once  V2  was 
achieved. 

2)  constant  altitude  -  Theta  was  varied  to  maintain  altitude  ±  20 
feet  about  the  target  altitude.  The  target  altitude  was  the 
altitude  at  which  the  escape  maneuver  began. 

3)  constant  theta  -  Theta  was  held  at  a  programmed  constant 
value  throughout  the  maneuver.  Theta  values  used  were  5®, 
10®,  and  15®. 

4)  constant  alpha  -  Theta  was  varied  to  maintain  a  constant 
AOA  value. 
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For  the  above  maneuvers,  a  pitch  rate  of  5o/sec  was  used  and 


maximum  thrust  was  maintained.  Also,  theta  was  reduced  any 
time  critical  AOA  was  exceeded. 

Two  variants  of  the  takeoff  scenario  were  analyzed  for 
performance  sensitivity.  Both  used  the  120,000  pound  P-3.  For  one, 
Viof  was  increased  18  knots  to  145  knots.  This  equated  to  an  increase 
of  the  rotate  speed  (Vr)  to  145  knots.  The  second  variant  moved  the 
microburst  center  from  1200  feet  to  2000  feet  from  the  point  of 
liftoff.  This  resulted  in  the  aircraft  gaining  airspeed  and  altitude 
before  encountering  the  severe  horizontal  windshear.  The  second 
scenario  closely  simulated  a  microburst  encounter  after  takeoff. 
MatLab«  programming  for  the  takeoff  scenario  is  listed  in 
Appendix  D. 

3.  P-3  On-Station 

The  on-station  microburst  encounter  scenario  used  a  P-3  at 
120,000  pounds  gross  weight  in  a  4-engine  loiter  configuration.  No 
wing  flaps  or  landing  gear  were  extended.  Only  aircraft  reaction  and 
performance  were  analyzed.  No  specific  escape  maneuver  was  used. 
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This  scenario  tried  to  emulate  an  autopilot  maintaining  altitude 


through  pitch  authority  and  the  flight  crew  controlling  thrust. 

The  microburst  encounter  began  with  the  P-3  at  200  feet  AGL 
and  5000  feet  from  the  microburst  center  (x=-5000  feet).  The  initial 
airspeed  was  210  knots  (prescribed  loiter  airspeed).  Theta  was  varied 
throughout  the  encounter  to  maintain  altitude.  Thrust  was 
maintained  at  loiter  power  until  40  knots  of  airspeed  was  lost.  At 
this  point,  maximum  thrust  was  applied  while  still  maintaining 
altitude.  The  scenario  was  ended  at  5000  feet  on  the  opposite  side 


of  the  microburst  center. 


III.  RESULTS  AND  DISCUSSION 


The  microburst  and  Inertial  reference  math  models  were 
validated  by  comparing  actual  flight  data  to  model  data.  The 
inertial  model  was  further  tested  for  the  stability  of  the  coupled 
differential  equation  scheme.  Any  weakness  noted  was  considered 
when  the  results  were  analyzed.  Each  critical  phase  of  flight  was 
studied  for  aircraft  response  when  performing  a  particular  escape 
maneuver.  Mathematical  results  were  combined  with  other 
observations  to  build  a  foundation  for  conclusions. 

A.  WINDSHEAR  AND  INERTIAL  REFERENCE  MODEL  VALIDATION 

The  wlndshear  math  model  was  developed  from  fitting  a 
vortex/source  flowfield  to  recorded  flight  data.  The  wind  field 
recorded  by  the  AA  539  DFDR  (Digital  Flight  Data  Recorder)  showed  a 
definite  vortex  flow.  However,  the  flight  data  from  the  DAL  191  DFDR 
showed  a  different  vortex  flow  arrangement.  In  addition,  a  strong 
outflow  at  the  surface,  and  Increased  vertical  sink,  required  the 
addition  of  source  flows  to  the  windshear  model.  This  use  of  point 
source  flows  led  to  a  much  closer  fit  of  the  original  vortex  model  to 
the  wind  field  DAL  191  experienced.  The  source  flows  have  no  range 
damping  terms.  Therefore,  the  DAL  191  emulation  windshear  model 
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becomes  invalid  at  distances  greater  than  6000  feet  (x=6000  feet) 
from  the  microburst  center.  All  analysis  was  easily  done  within  this 
distance  limit. 

It  is  important  to  note  that  the  aircraft  performance  results  are 
insensitive  to  the  exact  modeling  of  a  particular  windshear.  The 
rapid  change  of  a  headwind  to  a  tailwind  is  the  governing  factor  in  a 
microburst.  Secondary  are  the  vertical  down  drafts  that  can  be 
experienced  above  100  feet  AGL  (above  ground  level).  The  windshear 
math  model  contains  both  characteristics  to  the  same  degree  as  the 
windshear  experienced  by  DAL  191. 

Validation  of  the  inertial  reference  model  was  scrutinized  for 
proper  aircraft  response  to  changing  conditions  and  the  effect  of  the 
time  step  size.  Stability  of  the  solutions  obtained  from  the 
differential  equation  numerical  scheme  was  the  greatest  concern. 
Results  generated  from  the  inertial  model  were  studied  for  the  light 
P-3  under  stable  and  turbulent  conditions.  Also,  a  comparison  of  the 
3-engine  transport  model  was  made  with  the  DAL  191  DFDR  data 
during  the  final  seconds  of  the  fateful  flight. 

The  first  approach  was  to  chccK  the  response  of  a  light  P-3 
initially  stabilized  on  a  3o  descent  path  with  no  windshear.  The  input 
came  from  a  subroutine  that  would  not  vary  pitch  angle  (theta)  or 
thrust.  No  excursions  were  noted  in  descent  path,  airspeed,  or 
alpha,  as  shown  in  Figure  9.  The  time  step  used  was  I  second. 
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The  second  validation  was  to  observe  the  stability  of  the  inertial 
model  solution  with  different  size  time  steps.  A  time  interval  of  0.5 
seconds  was  the  target  time  step  for  running  all  analyses.  Therefore, 
time  steps  of  1,  0.5,  and  0.1  second  were  investigated.  The  light  P-3 
on  approach  to  landing  was  again  used.  This  time,  the  windshear 


Figure  9.  Light  P-3  on  a  stabilized  descent  path, 
model  was  incorporated  to  induce  changes  in  the  flight  environment. 
Aircraft  control  was  through  theta  and  thrust  inputs.  Control 
feedback  was  set  to  vary  theta  to  keep  within  *  100  feet  of  the 
descent  path  and  to  vary  thrust  to  maintain  Vr«f  within  i  5  knots. 

The  results  are  depicted  in  Figure  10.  Note  the  close  correlation 
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between  the  0.1  and  0.5  second  time  step  cases,  indicating  that  a  0.5 
second  time  step  provides  stability  and  adequate  accuracy  for  the 
desired  results. 


Figure  10.  The  effect  of  time  step  size  on  the  inertial 
reference  calculations. 

The  final  validation  was  comparing  the  actual  aircraft  response  to 
model  predictions.  DAL  191  DFDR  data  were  used  to  compare  an  L-lOll 
flight  path,  airspeed,  and  alpha  data  to  that  of  the  3-englne  heavy 
airline  transport  model.  Inputs  to  the  aircraft  control  were  the 
recorded  theta  and  thrust  of  DAL  191.  The  windshear  experienced  by 
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DAL  191  was  also  incorporated  for  this  comparison.  Figure  11 
graphically  displays  the  close  characterization  of  the  actual  flight 
path  by  the  model.  Airspeed  and  alpha  do  not  reflect  actual  aircraft 


Figure  11.  Comparison  of  aircraft  response  in  windshear 
between  the  DAL  191  Flight  L-1011  and  the  3-englne 
transport  model. 

response  as  closely  as  the  P-3  model,  but  do  follow  the  same  general 
response.  The  lag  In  airspeed  and  alpha  Is  probably  due  to  the  model 
Cl  and  Cp  equations  being  only  of  second  order. 
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B.  WEAKNESSES  AND  OTHER  INDEPENDENT  VARIABLES 

Some  weaknesses  in  the  math  models  and  unobserved 
independent  variables  exist  in  the  total  analysis.  The  most  notable 
are: 

1)  Cl  and  Cd  value  errors  at  angle  of  attack  near  stall. 

2)  Unobserved  effects  of  dynamic  pressure  and  AOA  upon  engine 
thrust. 

3)  Unobserved  effects  of  rain  and  turbulence  on  aircraft 
aerodynamic  performance. 

4)  Unobserved  effects  of  coupled  longitudinal  and  lateral  dynamic 
modes  excited  by  turbulence  and  pilot  induced  oscillations. 

None  of  the  weaknesses  or  unobserved  variables  is  believed  to 

significantly  Impact  the  results  of  the  analysis.  The  overall  concept 

was  to  compare  outcomes  of  different  microburst  escape  strategies 

given  the  same  parameters. 

A  weakness  in  the  aircraft  model  is  the  ability  to  predict  the 
effects  of  flow  separation  at  low  airspeed  and  high  angle  of  attack. 
These  effects  would  be  an  Increase  in  drag  and  a  decrease  in 
generated  lift.  Table  look-up  or  higher  order  equations  could  be  used 
for  the  Cp  and  Cl  expressions.  The  model  does  limit  the  maximum 

theta  not  to  exceed  maximum  angle  of  attack  for  that  aircraft. 
Keeping  this  weakness  in  mind,  the  results  obtained  from  the  math 
models  are  valid  when  determining  the  "best"  escape  maneuver. 
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Engine  thrust  is  somewhat  effected  by  airspeed  and  AO  A.  The 
turbojet  and  turbofan-type  engines  can  be  significantly  effected  while 
the  turboprop  is  effected  to  a  much  lesser  degree.  Only  nominal 
values  for  thrust  were  used  in  the  analysis. 

The  effects  of  turbulence  and  rain  on  aircraft  performance  were 
not  programmed  into  the  models.  Studies  by  Wlngrove  and  Bach 
[Ref:  2]  analytically  determined  that  the  rain  had  negligible  effect 
upon  DAL  191.  The  NTSB  Report  (National  Transportation  Safety 
Board)  came  to  the  same  conclusion.  What  must  be  considered,  is 
the  effect  of  turbulence  upon  a  pilot's  ability  to  control  the  aircraft. 
The  effect  of  turbulence  on  the  ability  of  a  pilot  to  execute  a 
particular  escape  maneuver  was  kept  in  mind  during  the  writing  of 
the  final  conclusions. 

All  the  aircraft  models  are  limited  to  one  axis  of  freedom.  The 
effects  of  coupled  lateral  and  directional  modes  upon  the  longitudinal 
response  were  not  modeled.  From  past  analysis  and  personal 
aviation  experience,  this  limitation  should  not  significantly  Impact 
aircraft  performance.  Again,  aircraft  control  may  become  difficult  if 
an  inherent  mode  is  excited  (such  as  the  phugoid  at  low  airspeed  emd 
high  angle  of  attack).  Escape  maneuvers  that  may  produce  such 
dynamic  modes  were  noted. 
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The  effect  of  changing  dynamic  pressure  on  stick  forces  was 
evaluated.  This  is  an  important  consideration  from  the  pilot’s 
perspective.  A  loss  of  20-30  knots  airspeed  can  induce  a  significant 
nose-down  pitch  force. 

C.  DATA  REDUCTION  AND  ANALYSIS 

Aircraft  performance  for  three  critical  flight  phases  was  evaluated 
when  exposed  to  a  DAL  191  type  microburst  windshear.  Approach  to 
landing  and  takeoff  flight  phases  were  scrutinized  for  the  best  escape 
maneuver  that  could  be  applied.  The  on-station  flight  phase  was 
examined  for  the  performance  of  a  P-3  while  encountering  the 
microburst  windshear  and  attempting  to  maintain  altitude. 
Calculated  data  were  recorded  and  converted  to  applicable  units  (eg., 
ft/s  to  knots).  The  data  were  then  presented  in  a  graphical  format 
for  analytical  and  subjective  comparison. 

For  all  cases,  tabulated  data  were  converted  to  graphical  form. 
Strip  graphs  for  key  dependent  variables  were  produced  for  each 
aircraft  performing  a  specific  escape  maneuver.  They  are  listed  in 
Appendices  F  and  G.  Combination  graphs  of  flight  path,  airspeed,  and 
specific  energy  compare  escape  maneuvers  for  each  aircraft  in  a 
given  microburst  encounter.  The  abscissa  axis  for  all  the  graphs 
represents  the  distance  from  the  microburst  center  in  feet.  A 


50 


three-step  approach  was  carried  out  to  deduce  the  “best"  escape 
manuever  for  a  particular  encounter  using  these  graphs. 

First,  the  inertial  reference  frame  of  flight  path  performance  was 
analyzed.  This  involved  using  the  combination  flight  path  graph.  A 
list  is  generated  containing  the  highest  to  lowest  altitude  obtained  at 
a  specific  point  from  the  microburst  center.  For  the  approach  to 
landing  scenario,  only  the  maneuvers  that  resulted  in  flight  paths 
staying  above  the  descent  path  (landing  glideslope)  were  chosen. 
Altitude  values  were  compared  for  each  valid  escape  maneuver  at 
x=5000  feet.  For  the  takeoff  scenario,  only  the  escape  maneuvers 
that  did  not  lead  to  ground  impact  were  listed.  Here,  altitude  values 
were  compared  at  x^^OOO  feet. 

The  second  step  was  to  cross  reference  each  selected  maneuver's 
airspeed  from  the  combination  airspeed  graph.  This  supplied  an 
insight  of  the  airmass  reference  performance.  The  airspeed  values  at 
the  point  of  interest,  and  throughout  the  microburst  encounter, 
allowed  a  subjective  analysis  of  the  validity  of  the  results.  Airspeeds 
that  fell  deep  within  the  power-on  stall  region  were  scrutinized  with 
the  strip  graphs  comparing  calculated  AO  A  values.  Any  maneuver 
that  resulted  in  a  very  low  airspeed,  high  AOA  condition  was  rejected 
from  the  list. 

The  third  step  evolved  using  the  specific  energy  combination 
graph  for  a  final  resolution.  Specific  energy  incorporates  both  the 
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Inertial  and  the  alrmass  reference.  The  eligible  maneuvers  that  had 
the  highest  specific  energy  value  at  the  point  of  interest  analytically 
are  the  "best"  escape  maneuvers  for  that  aircraft  configuration, 
given  that  type  of  microburst  encounter. 

i.  Approach  to  Landing  Microburst  Encounter  Analysis 
Figure  12,  Figure  13,  and  Figure  14  depict  the  results  obtained 
for  different  escape  maneuvers  executed  by  a  P-3  at  89,500  pounds 
gross  weight.  Table  6  compares  the  results  observed  at  x=5000  feet. 
The  constant  airspeed  and  12  unit  AOA  escape  maneuvers  were 
rejected  initially  due  to  resulting  flight  trajectories  below  the  descent 
path  (3o  landing  approach  glideslope)  as  depicted  in  Figure  12.  The  20 
unit  AOA  escape  maneuver  was  subsequently  rejected  because  of  low 
airspeed,  high  AOA  observed  in  Figure  13.  Figure  14  shows  that  the 
5o  theta  escape  maneuver  results  In  the  highest  specific  energy  value 
of  the  remaining  list. 


TABLE  6.  RELATIVE  VALUE  RANKING  FOR  A  P-3  AT  89,5000LBS. 


maneuver 

altitude 

airsoeed 

soecific  cncrav 

20  unit  AOA 

1 

rejected 

— 

const,  altitude 

2 

5 

5 

15o  theta 

3 

4 

4 

15  unit  AOA 

4 

3 

3 

lOo  theta 

5 

2 

2 

5o  theta 

6 

1 

1 

Figure  15,  Figure  16,  and  Figure  17  depict  the  results  obtained 
for  different  escape  maneuvers  executed  by  a  P-3  at  114,000  pounds 
gross  weight.  Table  7  compares  the  results  observed  at  x=5000  feet. 
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Figure  13.  Airspeed  data  comparison  of  a  P3  at  89,5001bs 
performing  different  escape  maneuvers  during  an 
approach  to  landing  microburst  encounter. 


54 


'5  .c 

“  5;  ^ 
-*•  c  ^ 
M3* 

oN 

O  ' 


in 


-  Q  -■5 

O  •*'  i 

*  <  *  ^ 

■c  -c  s; 

Ok  X  0*  . 

*  c  * 

•D  3  -O  M 

O  iO  10  O 

^  ^  O 


< 

O 

< 


o 

(M 


Figure  14.  Specific  energy  data  comparison  of  a  P3  at 
89,5001bs  performing  different  escape  maneuvers  diiring 
an  approach  to  landing  microbtarst  encounter. 
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The  constant  airspeed,  12  unit  AOA,  and  5o  theta  escape  maneuvers 
were  rejected  initially  due  resulting  flight  trajectories  below  the 
descent  path  as  depicted  in  Figure  15.  The  20  unit  AOA  and  constant 
altitude  escape  maneuvers  were  subsequently  rejected  because  of  low 
airspeeds,  high  AO  As  observed  in  Figure  16.  Figure  17  shows  that  the 
15  unit  AOA  escape  maneuver  results  in  the  highest  specific  energy 
value  of  the  remaining  list. 


TABLE  7.  RELATIVE  VALUE  RANKING  FOR  A  P-3  AT  114,000LBS. 


manguygr 

altitudg 

airspeed 

specific  energy 

20  unit  AOA 

1 

rejected 

— 

const,  altitude 

2 

rejected 

-  -  - 

15®  theta 

3 

3 

3 

10®  theta 

4 

2 

2 

15  unit  AOA 

5 

1 

1 

Figure  18,  Figure  .  19,  and  Figure  20  depict  the  results  obtained 
for  different  escape  maneuvers  executed  by  a  3-englne  heavy  airline 
transport.  Table  8  compares  the  results  observed  at  x=5000  feet. 
The  constant  airspeed,  5o  theta,  euid  iO®  theta  escape  maneuvers 


TABLE  8.  RELATIVE  VALUE  RANKING  FOR  A  3-EN6INE  AIRUNE 


TRANSPORT. 

maneuver 

altitude 

airspeed 

specific  energy 

const,  altitude 

1 

2 

2 

15®  theta 

2 

1 

1 

were  rejected  initially  due  resulting  flight  trajectories  below  the 
descent  path  as  depicted  in  Figure  18.  Figure  20  shows  that  the  15® 
theta  escape  maneuver  results  in  the  highest  specific  energy  value. 
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Figure  15.  Fl^ht  path  data  comparison  of  a  P3  at 
114,0001bs  performing  different  escape  memeuvers  during 
an  approach  to  landing  microhurst  encounter. 
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Figlire  16.  Airspeed  data  comparison  of  a  P3  at  114,0001bs 
performing  different  escape  maneuvers  during  an 
approach  to  landing  microhurst  encounter. 
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Figure  18.  Altitude  data  comparison  of  a  3-engine  heavy 
adrline  transport  performing  different  escape  maneuvers 
during  an  approach  to  landing  microhurst  encounter. 
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Figxare  19.  Airspeed  data  comparison  of  a  3-engine  heavy 
airline  transport  performing  different  escape  maneuvers 
during  an  approach  to  landing  microburst  encounter. 
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Figure  20.  Specific  energy  data  comparison  of  a  3-engine 
heavy  airline  transport  performing  different  escape 
maneuvers  during  an  approach  to  landing  microburst 
encounter. 


62 


Figure  21,  Figure  22,  and  Figure  23  depict  the  results  obtained 
for  different  escape  maneuvers  executed  by  a  T-44.  Table  9 
compares  the  results  observed  at  x=5000  feet.  The  constant 
airspeed,  10  unit  AOA,  and  16  unit  AOA  escape  maneuvers  were 
rejected  initially  due  to  resulting  flight  trajectories  below  the  descent 
path  as  depicted  in  Figure  21.  The  25  unit  AOA  and  constant  altitude 
escape  maneuvers  was  subsequently  rejected  because  of  low  airspeeds 
and  high  AOAs  observed  in  Figure  22.  Figure  23  shows  that  the  5o 
theta  escape  maneuver  results  in  the  highest  specific  energy  value  of 
the  remaining  list. 


TABLE  9.  RELATIVE  VALUE  RANKING  FOR  A  T-44. 


maneuver 

altitude 

airspeed 

specific  enercv 

const,  altitude 

- 1 

rejected 

— 

15°  theta 

2 

3 

3 

25  unit  AOA 

3 

rejected 

— 

lOo  theta 

4 

2 

2 

5o  theta 

5 

1 

1 

Figure  24,  Figure  25,  and  Figure  26  depict  the  results  obtained 
for  different  escape  maneuvers  performed  by  a  P-3  at  89,500  pounds 
gross  weight  operating  on  3  engines.  Table  10  compares  the  results 
observed  at  x*5000  feet.  The  constant  airspeed,  12  unit  AOA,  and  5o 
theta  escape  maneuvers  were  rejected  initially  due  to  resulting  flight 
trajectories  below  the  descent  path  as  depicted  in  Figure  24.  The  20 
unit  AOA  and  constant  altitude  escape  maneuvers  were  subsequently 
rejected  because  of  low  airspeeds  and  high  AOAs  observed  in 
Figure  25.  The  150  theta  escape  maneuver  was  also  rejected  for  the 
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Figure  21.  Altitude  data  compeurison  of  a 
performing  different  escape  meuieuvers  during 
approach  to  landing  microburst  encounter. 


T-44 

an 
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Figxare  22.  Airspeed  data  competrison  of  a 
performing  different  escape  maneuvers  during 
approach  to  landing  microburst  encounter. 


T-44 

an 
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Figure  23.  Specific  energy  data  comparison  of  a  T-44 
performing  different  escape  maneuvers  during  an 
approach  to  landing  microburst  encoxjinter. 
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Figxire  25.  Airspeed  data  comparison  of  a  P-3  at 
89^5001bs  operating  on  3  engines  performing  different 
escape  maneuvers  during  an  approach  to  landing 
microhurst  encounter. 
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Figlire  26.  Specific  energy  data  compeurison  of  a  P-3  at 
89,5001135  operating  on  3  engines  performing  different 
escape  meuieuvers  during  an  approach  to  leuiding 
microburst  encounter. 
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same  reasons.  Although  the  airspeed  is  only  within  the  flow 
separation  region  (steady  state  stall  buffet),  analysis  of  the  AOA  (see 
Figure  F34)  shows  that  critical  alpha  was  sustained  on  the  latter 
parts  of  the  maneuver.  Figure  26  shows  that  the  15  unit  AOA  escape 
maneuver  results  in  the  highest  specific  energy  value  of  the 
remaining  list. 


TABLE  10.  RELATIVE  VALUE  RANKING  FOR  A  P-3  AT  89,5000LBS 
OPERATING  ON  3  ENGINES. 


maneuver 

altitude 

airsoeed 

soecific  enerav 

const,  altitude 

1 

rejected 

— 

15o  theta 

2 

rejected 

— 

20  unit  AOA 

3 

rejected 

— 

lOo  theta 

4 

2 

2 

15  unit  AOA 

5 

1 

1 

2.  Takeoff  Microburst  Encounter  Analysis 


Figure  27,  Figure  28,  and  Figure  29  depict  the  results  obtained 
for  different  escape  maneuvers  performed  by  a  P-3  at  90,0001bs  gross 
weight.  Altitude,  airspeed,  and  specific  energy  values  are  compared 
at  x=4000ft.  All  six  escape  maneuvers  resulted  in  ground  plane 
clearance  as  depicted  in  Figure  27.  The  15  unit  AOA  escape  maneuver 
was  subsequently  rejected  because  of  low  airspeed,  high  AOA  as 
observed  in  Figure  28.  Figure  29  shows  that  all  escape  maneuvers, 
except  the  15  unit  AOA  escape  maneuver,  resulted  in  a  grouped  single 
value.  This  shows  that  all  relevant  escape  maneuvers  provide  the 
same  final  specific  energy. 
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Figure  28 
90,0001bs 
during  a 


Figure  29.  Specific  energy  compeu'ison  of  a  P-3  at 
90,0001bs  performing  different  escape  maneuvers 
during  a  takeoff  microburst  encounter. 
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Figure  30,  Figure  31,  and  Figure  32  depict  the  results  obtained 
for  different  escape  maneuvers  performed  by  a  P-3  at  120,000  pounds 
gross  weight.  Value  comparison  is  made  at  x=40C0  feet.  The 
constant  airspeed,  constant  altitude,  12  unit  AOA,  15  unit  AO  A,  and  5o 
theta  escape  maneuvers  were  rejected  initially  due  ground  impact  as 
depicted  in  Figure  30.  The  15o  theta  escape  maneuver  was 
subsequently  rejected  because  of  low  airspeed  and  high  AOA  as 
observed  in  Figure  31.  The  only  viable,  however  marginal, 
performance  observed  was  for  the  lOo  theta  escape  maneuver. 
Figure  31  shows  that  the  10°  theta  maneuver  resulted  in  an  airspeed 
drop  of  30  Knots  below  V2  at  one  point  followed  by  a  sustained 

airspeed  of  20  Knots  below  V2. 

Included  with  the  120,000  pound  P-3  data  is  the  resulting  flight 
performance  when  Vr  was  Increased.  A  large  gain  in  altitude  and 

specific  energy  over  other  performance  profiles  is  observed. 
Increasing  the  rotate  speed  by  18  Knots  coupled  with  flying  a  constant 
lOo  theta  provides  a  2  fold  Increase  In  altitude  and  a  21%  increase  in 
specific  energy.  This  Increase  in  specific  energy  is  important  when 
compared  to  the  spread  of  specific  energy  values  of  different  escape 
maneuvers. 
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Figure  31.  Airspeed  comparison  of  a  P-3  at 
120,0001bs  performing  different  escape  maneuvers 
during  a  takeoff  microburst  encounter. 
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Figure  33,  Figure  34,  and  Figure  35  depict  the  results  obtained 
for  different  escape  maneuvers  performed  by  a  P-3  at  135,0001bs 
gross  weight.  It  can  be  graphically  observed  that  a  P-3  loaded  to 
maximum  takeoff  weight  is  at  the  mercy  of  a  strong  microburst 
encountered  at  takeoff.  Although  the  15o  theta  maneuver  misses  the 
ground,  the  low  airspeed  and  high  AOA  experienced  disqualifies  it  as 
viable. 

Figure  36,  Figure  37,  and  Figure  38  depict  the  results  obtained 
for  different  escape  maneuvers  performed  by  a  T-44.  Altitude, 
airspeed,  and  specific  energy  values  are  compared  at  x=4000  feet.  All 
five  escape  maneuvers  resulted  in  ground  plane  clearance  as  depicted 
in  Figure  36.  All  five  escape  maneuvers  provided  adequate  margin  of 
airspeed  and  AOA  as  depicted  in  Figure  37.  Note  that  power-on  stall 
speed  for  the  T-44  is  less  than  90  knots.  Figure  38  shows  that  all 
escape  maneuvers  resulted  In  a  grouped  single  value  of  specific 
energy.  As  seen  above,  all  escape  maneuvers  for  a  given  aircraft 
provide  the  same  final  specific  energy  during  a  takeoff  microburst 
encounter.  The  only  exception  is  when  Vr  is  increased. 

The  results  of  one  variant  to  the  takeoff  encounter  are 
depicted  in  Figure  39,  Figure  40,  and  Figure  41.  Here,  the  microburst 
center  is  moved  an  added  800  feet  from  the  liftoff  point.  This  results 
in  a  delayed  encounter  with  the  severest  part  of  the  wlndshear.  A 
P-3  at  120,00  pounds  gross  weight  was  used  to  perform  the  different 
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Figure  34.  Airspeed  comparison  of  a  P-3  at 
135,0001bs  performing  different  escape  maneuvers 
during  a  takeoff  microburst  encounter. 
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Figure  41.  Specific  energy  comparison  of  a  P-3  at 
120,0001bs  performing  different  escape  maneuvers 
during  an  after  takeoff  microburst  encounter. 
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escape  maneuvers  including  the  Increased  Vr  technique.  Table  11 


compares  the  results  observed  at  x=4000  feet.  The  12  unit  AOA 
escape  maneuver  was  rejected  initially  due  to  ground  impact  as 
depicted  in  Figure  39.  The  15  unit  AOA  and  15o  theta  escape 
maneuvers  were  subsequently  rejected  because  of  low  airspeeds  and 
high  AOAs  observed  in  Figure  40.  Figure  41  shows  that  the  close 
grouping  of  the  specific  energy  values  still  exists  for  the 
non-lncr eased  Vr  maneuvers.  Again,  the  Increased  Vr  resulted  in  a 

significant  Improvement  in  altitude  and  specific  energy. 


TABLE  11.  RELATIVE  VALUE  RANKING  FOR  A  P-3  AT  120,000LBS 
WITH  AN  AFTER  TAKEOFF  ENCOUNTER. 


maneuver 

altitude 

airspeed 

specific  energy 

increased  Vr 

1 

2 

1 

15o  theta 

2 

rejected 

— 

lOo  theta 

3 

4 

2 

const,  altitude 

4 

3 

3 

15  unit  AOA 

5 

rejected 

_  _  _ 

const,  airspeed 

6 

1 

4 

3.  On-Station  Encounter  Analysis 

The  analysis  showed  that  a  P-3  at  120,0001bs  gross  weight 
successfully  navigated  a  strong  microburst  during  loiter  operations. 
Figure  42  graphically  depicts  the  results.  Altitude  deviation  was  no 
greater  than  ±  20ft.  Theta  input  never  exceeded  10  o.  AOA  remained 
below  critical  angle.  Additional  thrust  was  not  applied  until  40  knots 
of  airspeed  was  lost  In  the  encounter  (a  drop  from  210  to  170  knots). 
Power  was  added  at  x=2700  feet  and  a  time  lapse  of  23sec  from  the 
Initial  point  of  the  encounter. 
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Figure  42.  P-3  at  120,000lbs  gross  weight  encountering  a 
microburst  windshear  while  on-station  loiter. 


D.  OBSERVATIONS 


The  results  obtained  from  the  analytical  analysis  for  each 
encounter  are  integrated  with  other  relative  information  to  answer 
the  three  posed  questions: 

1)  What  is  the  optimum  microburst  escape  procedure  given  the 
type  of  encounter  and  the  available  flight  performance 
information? 

2)  Is  the  optimum  escape  procedure  effected  by  change  of  gross 
weight  or  available  thrust? 

3)  What  flight  instrument  indications,  flight  control  feeling,  and 
dynamic  response  would  be  expected  during  the  optimum 
escape  procedure? 

The  answers  to  these  questions  are  directed  to  the  P-3.  However, 
comparison  is  made  to  the  T-44  and  3-englnc  heavy  airline  transport 
where  appropriate.  Other  Important  issues  must  also  be  addressed 
for  comparison  of  results  to  other  wlndshear  studies.  These  issues 
include  the  impact  of  wing  loading  (W/S),  thrust  to  weight  (T/W),  and 
early  liftoff  speeds  for  the  P-3. 

1.  Available  Aircraft  Flight  Instruments 

The  available  flight  information  and  type  of  presentation  are 
extremely  influential  on  the  choice  of  the  optimum  escape  procedure. 
Both  the  P-3  and  T-44  have  "conventional"  flight  instruments  [Ref:  8 
and  Ref:  9].  Conventional  implies  gyroscopic-stabilized  attitude  and 
heading  indication  in  combination  with  pitot-static  airspeed,  altitude, 
and  vertical  speed  indication  (VSI).  Angle  of  Attack  (AOA)  indication 
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is  installed  on  the  P-3  and  T-44.  However,  the  AOA  is  heavily 
damped  and  is  primarily  designed  for  steady-state  flight  conditions. 
Standard  radio  navigation  equipment  is  installed  on  both  aircraft 
encompassing  an  Integrated  Landing  System  (ILS),  VHF  Omni  Radio 
(VOR),  and  Tactical  Navigation  (TACAN).  The  P-3  in  addition  has 
Inertial  Navigation  Systems  (INS),  but  no  Information  Is  presented  at 
the  flight  station.  Both  aircraft  incorporate  autopilots  and  flight 
directors.  The  P-3  autopilot  is  controlled  by  Control  Wheel  Steering 
(CWS)  input  only.  The  T-44  autopilot  is  controlled  through  thumb 
wheels  and  can  couple  to  the  flight  director.  Flight  directors  on  both 
aircraft  can  provide  navigational  steering.  However,  the  primary 
design  of  the  P-3  flight  director  is  tactical  while  the  primary  design 
of  the  T-44  flight  director  is  navigational.  In  summary,  both  the  P-3 
and  T-44  have  Attitude  Indicators  (AI),  damped  AOA  indication,  and 
pitot-static  instruments  for  reference  while  executing  an  escape 
maneuver. 

2.  Optimal  Escape  Procedure 

Flying  a  constant  theta  with  reference  to  the  AI  seems  to  be 
the  optimal  escape  procedure  during  a  microburst  encounter.  This 
conclusion  Is  supported  analytically  and  qualitatively  for  the  P-3  and 
T-44.  Analytically,  superior  performance  was  obtained  In  a  few 
circumstances  from  constant  AOA  maneuvers  with  an  alpha  value 
between  approach  and  stall.  However,  qualitatively  the  AOA 
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indication  in  both  aircraft  does  not  lend  itself  as  a  viable  reference 
for  two  reasons.  First,  the  AO  A  indication  may  be  significantly  in 
error  because  of  instrument  lag  in  turbulent  conditions.  Second,  P-3 
and  T-44  pilots  seldom  use  AOA  indication  during  the  approach  phase. 
The  transition  from  attitude/airspeed  to  attitude/ AOA  reference 
during  an  intense  situation  as  a  microburst  encounter,  is  more  than 
can  be  expected  from  a  pilotage  standpoint.  The  only  viable 
alternative  to  flying  constant  AOA  was  a  constant  theta  maneuver. 
Note  that  quantitatively  in  all  circumstances,  a  constant  theta 
maneuver  was  superior  to  or  closely  matched  with  any  other 
attempted  escape  maneuver. 

Theta  values  ranging  from  5o  to  150  were  Identified  as  optimal 
for  an  approach  to  landing  encounter.  The  theta  value  of  10®,  if  not 
optimal,  provided  suitable  recovery  for  the  P-3  in  the  configurations 
analyzed.  The  light,  the  heavy,  and  the  3-engine  P-3  successfully 
navigated  the  microburst  encounter  utilizing  the  lOo  theta  escape. 
This  result  is  important  from  the  standpoint  that  lOo  is  easy  to 
remember  and  easy  to  read  on  the  AI.  Also,  lOo  theta  worked  well 
for  the  T-44  in  the  approach  configuration  (approach  flap  setting).  A 
theta  of  15o  was  optimum  for  the  T-44  In  the  takeoff  configuration 
(flaps  up).  This  is  important  because  future  P-3  pilots  receive  their 
first  training  in  T44s. 
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3.  The  Effect  Of  Specific  Energy  On  Survival  Probability 

Escaping  from  a  takeoff  microburst  encounter  has  a  lower 
probability  for  success  compared  to  the  approach  to  landing 
encounter.  The  type  of  escape  maneuver  selected  in  a  takeoff 
encounter  has  a  lesser  effect  on  the  outcome.  This  result  is  caused 
by  the  excess  airspeed  and  additional  altitude  available  during  the 
approach  phase.  (Note  that  the  P-3  community  flies  an  unusually 
fast  approach  speed  compared  to  their  civilian  counterparts.  The  P-3 
approach  Vrai  airspeed  is  determined  from  1.35V6+5  knots  [Ref.  20], 

which  equates  to  1.4Vs  for  a  gross  weight  of  114,000  pounds.  The 
industry  standard  for  Vref  is  I.3Vs+an  added  factor  as  dictated  by 

the  type  of  aircraft.)  Specific  energy  can  be  used  for  an  interesting 
comparison  between  an  approach  to  landing  and  a  takeoff  microburst 
encounter.  A  P-3  weighing  89,000  pounds  has  specific  energy  values 
of  1230  feet  and  1040  feet  at  x=500  feet  and  x=4000  feet  respectively 
during  an  approach  to  landing  encounter  (reference  Figure  14). 
During  a  takeoff  encounter,  a  P-3  at  90,000  pounds  has  a  specific 
energy  values  of  810  feet  and  980  feet  at  the  same  x  distances  from 
the  microburst  center  (reference  Figure  29).  Both  executed  a  lOo 
theta  escape  maneuver  and  both  successfully  navigated  the 
microburst.  Note  that  the  approach  to  landing  phase  lost  energy 
while  the  takeoff  phase  gained  energy. 
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An  important  aspect  of  specific  energy  comparison  is  seen  in 
the  increased  rotate  speed  data.  Figure  4i  graphically  depicts  a 
120,000-pound  P-3  starting  with  an  additional  specific  energy  of  200 
feet  compared  to  the  normal  rotate  speed  profiles.  The  exit  specific 
energy  was  175  feet  greater  than  the  comparable  lOo  theta 
maneuver.  The  175  feet  were  translated  almost  completely  to 
altitude  (reference  Figure  30  and  Figure  31).  Data  support  that  an 
increased  rotate  speed  is  far  more  beneficial  than  any  particular 
escape  maneuver  during  a  taKeoff  microburst  encounter. 

4.  Stick  Force  vs.  Off-Trim  Airspeed  For  The  P-3 

The  expected  change  in  flight  control  "feel"  for  the  P-3  is  light 
compared  to  large  aircraft.  Flight  tests  show  that  stick  force  versus 
change  from  trim  airspeed  results  in  a  shallow  gradient  for  this  size 
aircraft  [Ref.  I8:p.  37].  In  the  takeoff/approach  flap  configuration  (l8o), 
with  an  aft  C.G.  of  293  M.A.C.,  a  stick  force  of  an  11-pound  pull  is 
needed  with  a  30-knot  decrease  from  trim  airspeed.  Stick  force 
increases  to  a  16  pound  pull  with  the  most  forward  C.G.  of  16*  M.A.C. 
Flight  test  data  [Ref.  18]  also  show  that  the  stick  force  to  trim 
airspeed  gradient  is  not  appreciably  effected  by  landing  gear  position 
or  gross  weight.  The  same  flap  setting  (takeoff/approach)  is  used  for 
landing  approach  and  takeoff.  Thus,  the  elevator  "feel"  is  expected  to 
be  the  same  for  either  approach  to  landing  or  takeoff  microburst 
encounters. 
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The  P-3  Flight  Simulator  (Device  2F87F)  was  used  to  confirm 
the  flight  test  data  relating  to  off-speed  to  stick  force  and  landing 
gear  drag.  The  simulator  allows  one  or  more  flight  parameters  to  be 
frozen.  This  feature  is  used  by  freezing  the  altitude  to  500  feet  and 
trimming  the  aircraft  with  maximum  power  at  reference  airspeed. 
The  airspeed  is  then  reset  and  frozen  to  the  low  airspeed  value 
expected  during  a  microburst  encounter.  The  control  force  required 
to  maintain  a  pitch  attitude  was  then  directly  measured  and  printed 
out.  The  P-3  flight  simulator  data  (Appendix  H)  showed  that  the 
elevator  pull  needed  at  Vtrim  minus  30kts  decreased  by  5  pounds 

when  full  power  was  applied.  Level  acceleration  maneuvers  were 
performed  with  the  landing  gear  extended,  retracted,  and  during 
retraction.  Simulator  data  showed  that  the  landing  gear  retraction 
cycle  does  not  increase  drag. 

5.  Weight,  Wing  Loading,  and  Thrust  to  Weight  Effects 

It  Is  obvious  from  the  analysis  that  a  heavy  P-3  is  effected  by 
microburst  windshear  to  a  higher  degree  than  a  light  P-3.  The  T-44 
was  relatively  unaffected  by  the  windshear.  The  question  arises 
whether  weight,  wing  loading,  or  thrust  to  weight  has  the  largest 
effect  on  microburst  survival.  Table  12  compares  the  weight  (Wt), 
wing  loading  (W/S),  thrust  to  weight  (T/W),  and  specific  energy  values 
for  three  different  aircraft  during  an  approach  to  landing  microburst 
encounter.  The  specific  energy  loss  (Es  loss)  column  shows  the  energy 


95 


(In  feet)  lost  during  the  encounter.  The  specific  energy  spread  (E* 

spread)  column  shows  the  spread  (In  feet)  of  the  different  escape 
techniques  measured  after  the  encounter.  Apparently,  the  higher 
the  thrust  to  weight  ratio,  the  less  energy  is  lost.  Light  wing  loading 
equates  to  an  individual  escape  procedure  having  less  effect  upon  the 
result.  Weight,  as  a  separate  parameter,  is  irrelevant.  In  other 
words,  increased  thrust  to  weight  decreases  the  effect  of  the 
windshear  upon  the  aircraft  and  light  wing  loading  allows  increased 
efficiency  of  energy  transformation. 


TABLE  12.  WEIGHT,  WING  LOADING,  AND  THRUST  TO  WEIGHT 
EFFECTS  ON  Pg  VALUES. 


aircraft 

Wt  (lbs) 

W/S 

im 

Eg  less 

Ee  soread 

P-3 

89,500 

68.85 

0.37 

380 

145 

P-3 

114,000 

88.08 

0.29 

455 

155 

T-44 

8,280 

39.43 

0.37 

285 

50 

3-eng. 

362,000 

79.07 

0.35 

345 

100 

hvy. 

6.  Early  Liftoff  Speed  For  The  P-3 

The  analyses  did  not  look  at  the  effect  of  a  microburst 
windshear  encountered  during  a  takeoff  roll.  The  FAA  Windshear 
Training  Aid  [Ref.  14]  points  out  that  If  windshear  is  encountered  past 
Vi,  the  aircraft  should  be  committed  to  flight  no  later  than  with 

2000  feet  of  runway  remaining.  Although  the  airspeed  may  fall  and 
stay  below  Vr,  most  aircraft  can  go  airborne.  In  this  regard,  a 

theoretical  liftoff  speed  was  calculated  for  the  P-3.  Equations  38,  40, 
and  42  were  used  in  conjunction  to  solve  for  V.  The  Clo  data  were 


96 


extrapolated  from  flight  test  data  [Ref.  18].  Table  13  compares  the 
recommended  Vr  speed  [Ref.  20]  to  the  theoretical  liftoff  speed  for  a 

rotate  pitch  attitude  of  lOo.  Ground  effect  was  not  considered  in  the 
calculations. 


TABLE  13.  THEORETICAL  LIFTOFF  SPEEDS 


weight  (lbs) 
80,000 
90,000 
100,000 
110,000 
120,000 
130,000 
140,000 


Yr-iUsl 

theoretical  liftoff  (kts) 

115 

90 

115 

95 

115 

101 

117 

105 

123 

no 

129 

115 

134 

119 

8.  Summary 

The  original  three  questions  are  answered  for  the  P-3: 

1)  The  optimum  approach  to  landing  microburst  escape 
procedure  Is  to  set  and  maintain  maximum  power  while 
simultaneously  setting  a  lOo  pitch  attitude.  For  the  takeoff 
microburst  encounter,  increase  rotate  speed  to  140  knots 
then  pitch  to  and  maintain  lOo. 

2)  The  optimum  escape  procedures  remain  the  same  for  all 
gross  weights.  The  approach  to  landing  encounter  escape 
procedure  is  the  same  for  four  or  three  operating  engines. 

3)  The  airspeed  will  rapidly  decay  and  remain  abnormally  low 
during  a  microburst  penetration.  An  elevator  force  of  5  to 
10  pounds  can  be  expected  to  maintain  a  lOo  pitch  attitude 
with  full  power.  Flight  tests  show  that  the  P-3  has  no 
unusual  short  period,  phugoid,  nor  cross-couple  dynamics  to 
contend  with. 
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The  answers  for  the  T-44  should  be  similar  to  the  P-3.  The  lOo  pitch 
attitude  was  optimal  for  the  approach  to  landing  encounter.  A 
pitch  attitude  of  15o  is  optimal  for  takeoff,  the  difference  being  that 
the  flaps  are  up  in  the  takeoff  configuration.  Not  enough 
performance  data  were  available  to  predict  the  expected  stick  forces. 
However,  personal  flight  experience  has  shown  that  the  T-44  exhibits 
longitudinal  forces  and  responses  very  similar  to  the  P-3. 
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IV.  CONCLUSIONS 


The  aim  of  this  thesis  was  to  analytically  produce  microburst 
escape  flight  procedures.  Study  was  directed  at  two  turboprop 
aircraft  the  Navy  routinely  operates.  The  P-3  and  7-44  were  closely 
scrutinized  for  the  optimal  escape  procedure  if  confronted  with  a 
microburst  penetration.  The  microburst  modeled  was  patterned 
after  that  encountered  by  a  Delta  Airlines  L-lOll  at  Dallas/Ft.  Worth 
on  August  2,  1985.  The  conclusions  drawn  arc  presented  within  the 
format  of  the  FAA  Windshear  Training  Aid  [Ref.  14].  Viable  flight 
procedures  and  associated  precautions  are  presented  for  both  aircraft. 
The  precautions  should  be  adhered  to  if  an  encounter  with  a 
microburst  windshear  is  a  possible  expectation.  NATOPS  change 
recommendations  reflecting  these  conclusions  will  be  forwarded  to  the 
appropriate  aircraft  model  managers. 

A.  MICROBURST  ESCAPE  PROCEDURES  FOR  THE  P-3 

A  target  pitch  attitude  of  10  o  is  optimal  for  all  weights  and  with 
four  or  three  engines  operating.  This  target  pitch  attitude  is  used  for 
approach  to  landing  and  taKeoff  microburst  encounters  with  the  flaps 
in  the  Takeoff/ Approach  position.  The  airspeed  at  which  the  P-3  lifts 
off  is  the  critical  factor  in  the  takeoff  scenario. 
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1.  Approach  To  Landing  Microburst  Encounter 

Precautions; 

1)  Avoid  thunderstorm  conditions.  Delay  the  approach  if 
possible. 

2)  Use  TaKeoff/Approach  flap  setting.  Do  not  extend  flaps 
to  Land  position  until  the  runway  is  made. 

3)  Use  a  precision  approach  procedure. 

4)  The  approach  speed,  Vref,  should  be  equal  to  1.35V3  +  5 
knots.  Additional  airspeed  is  not  normally  warranted. 
This  Vrcf  will  provide  a  comfortable  margin  above  stall 
speed.  Note  that  additional  airspeed  increases  landing 
distance  (5  knots  faster,  lOf  farther). 

5)  Determine  expected  descent  rate. 

6)  Consider  using  the  autopilot  (ASW-31)  with  the  altitude 
hold  switch  in  the  off  position.  The  autopilot  will 
provide  wing  rock  damping  and  pitch  attitude 
augmentation  in  the  face  of  turbulence.  It  will  also 
provide  pitch  attitude  hold  if  a  microburst  escape 
maneuver  is  warranted.  Insure  autopilot  disconnect 
passing  200  feet  AGL. 

7)  Attain  a  stabilized  airspeed  approach  before  passing 
1000  feet  AGL.  Minimize  power  lever  movement 
beyond  this  point.  Maintain  glideslope  with  pitch.  The 
airspeed  indicator  will  serve  as  a  windshear  indicator. 

8)  Strong  consideration  should  be  given  to  executing  the 
microburst  escape  procedure  if  one  of  the  indications  is 
observed: 

a)  a  rapid  and  sustained  airspeed  loss  of  20  knots  below 

Vrafl 

b)  a  descent  rate  500  feet  per  minute  greater  than  the 
predetermined  value; 

c)  greater  than  1  dot  low  on  the  ILS  glideslope  associated 
with  airspeed  10  knots  below  Vref;  or 

d)  a  "well  below  glideslope"  call  on  a  PAR  associated 
with  airspeed  10  knots  below  Vref 
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Escape  procedure: 

1)  Apply  maximum  power  (power  levers  to  the  stops). 

2)  Set  and  maintain  a  pitch  attitude  of  IQo  on  the  AI.  Do 
not  attempt  to  recover  airspeed. 

3)  Once  a  positive  rate  of  climb  is  established,  select 
landing  gear  up. 

4)  Do  not  raise  the  flaps  until  the  airspeed  has  increased 
above  140  knots  indicating  exit  from  the  windshear. 

2.  Takeoff  Microburst  Encounter 

Precautions: 

1)  Delay  the  takeoff  if  able. 

2)  Select  the  longest  suitable  runway. 

3)  Perform  takeoff  planning  for  adverse  conditions  as 
prescribed  by  the  NATOPS  Flight  Manual. 

4)  Increase  rotate  airspeed  to  140  knots,  or  the  airspeed 
that  is  attained  with  2000  feet  of  runway  remaining. 
This  is  determined  by  using  the  Four  Engine 
Acceleration  Chart  in  NATOPS  and  following  this 
procedure: 

a)  Subtract  two  thousand  feet  from  the  available 
runway  length. 

b)  Enter  the  chart  with  the  adjusted  runway  length, 
pressure  altitude,  temperature  and  gross  weight.  Exit 
the  chart  with  the  corresponding  airspeed  value.  Use 
140  knots  if  the  chart  value  is  higher. 

c)  Corrections  for  runway  slope,  winds,  or  standing 
water  is  not  required  for  THIS  prediction. 

5)  Thoroughly  brief  the  takeoff  procedure,  voice  calls,  and 
windshear  indications  among  the  flight  crew. 
Penetration  of  a  windshear  during  takeoff  will  be 
indicated  by  a  loss  of  airspeed  or  no  airspeed 
acceleration. 
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Procedure: 


1)  Abort  the  takeoff  if  windshear  is  indicated  before 
reaching  refusal  airspeed.  Reaching  refusal  speed 
ground  roll  distance  with  airspeed  significantly  below 
OR  above  refusal  speed  indicates  windshear. 

2)  Continue  the  takeoff  if  windshear  is  experienced  after 
refusal  speed.  Rotate  the  nose  when; 

a)  predicted  increased  rotate  speed,  usually  140  knots,  is 
attained,  or 

b)  2000  feet  of  runway  remain. 

3)  Rotate  to  a  pitch  attitude  of  10 1>.  DO  NOT  delay  rotating 
the  nose  because  of  low  airspeed. 

4)  Increase  pitch  attitude  toward  15o  if  ground  Impact  is 
imminent. 

5)  Raise  the  landing  gear  when  above  100  feet  AGL  with  a 
positive  rate  of  climb. 

6)  Do  not  raise  the  flaps  until  normal  climb  airspeed  is 
regained,  indicating  clear  of  the  microburst. 

3.  On-Station  Loiter  Microburst  Encounter 

No  immediate  effect  on  the  aircraft  performance  was 

observed.  However,  certain  precautions  should  be  taken  and 

procedures  followed  If  flying  through  convective  activity  at  low 

altitude: 

1)  Use  the  autopilot  altitude  hold  In  dual  axis  mode. 
Note  that  if  the  autopilot  is  not  used,  the  pilot  can 
expect  abnormal  elevator  control  force  changes. 

2)  Minimize  power  lever  movement.  Do  not  pull  off 
power  for  high  airspeed  indication. 
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3)  A  penetration  of  a  microburst  will  be  Indicated  by  a 
significant  loss  of  airspeed.  The  first  and  foremost 
reaction  should  be  the  addition  of  full  power.  Monitor 
autopilot  input  and  altitude  hold. 

4)  If  a  climb  is  deemed  warranted,  set  lOo  pitch  attitude 
using  the  CWS  function  of  the  autopilot. 


B.  MICROBURST  ESCAPE  PROCEDURES  FOR  THE  T-44 

Target  pitch  attitude  escape  procedures  are  Just  as  effective  for  a 
light  turboprop  as  the  T-44.  In  its  case,  the  target  pitch  values  are 
different  between  takeoff  and  landing.  This  change  is  owed  to 
different  flap  settings.  Many  of  the  same  precautions  and 
procedures  are  the  same  between  the  T-44  and  P-3. 

1.  Approach  To  Landing  Microburst  Encounter 
Precautions: 

1)  Avoid  thunderstorm  conditions.  Delay  the  approach  if 
possible. 

2)  Use  Approach  flap  setting.  Do  not  extend  flaps  to  full 
down  position  unless  required  and  until  the  runway  is 
made. 

3)  Use  a  precision  approach  procedure. 

4)  The  approach  speed,  Vrvtf  should  be  equal  to  120  knots. 
Additional  airspeed  is  not  warranted.  This  Vrgf  will 
provide  a  comfortable  margin  above  stall  speed.  Note 
that  additional  airspeed  Increases  landing  distance  (5 
knots  faster,  10*  farther). 

5)  Determine  expected  descent  rate. 
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6)  Consider  using  the  autopilot  in  the  coupled  mode  if 
executing  an  ILS.  However,  close  monitoring  of  the 
autopilot  is  required  because  of  abnormal  elevator  pitch 
forces  expected  in  a  windshear.  Insure  autopilot 
disconnect  passing  200  feet  AGL. 

7)  Attain  a  stabilixed  airspeed  approach  before  passing 
1000  feet  AGL.  Minimize  power  lever  movement 
beyond  this  point.  Maintain  glldeslope  with  pitch.  The 
airspeed  indicator  will  serve  as  a  windshear  indicator. 

8)  Strong  consideration  should  be  given  to  executing  the 
microburst  escape  procedure  if  one  of  the  Indications  is 
observed; 

a)  a  rapid  and  sustained  airspeed  loss  of  15  knots  below 
Vrcf) 

b)  a  descent  rate  500  feet  per  minute  greater  than  the 
predetermined  value; 

c)  greater  than  1  dot  low  on  the  ILS  glldeslope  associated 
with  airspeed  10  knots  below  Vref;  or 

d)  a  "well  below  glldeslope"  call  on  a  PAR  associated 
with  airspeed  10  knots  below  Vref. 

Escape  procedure: 

1)  Apply  maximum  power  (power  levers  to  the  stops). 

2)  Set  and  maintain  a  pitch  attitude  of  lOo  on  the  AI.  DO 
NOT  attempt  to  recover  airspeed, 

3)  Once  a  positive  rate  of  climb  is  established,  select 
landing  gear  up. 

4)  Do  not  raise  the  flaps  until  the  airspeed  has  increased 
above  120  knots,  indicating  exit  from  the  windshear. 


104 


2.  Takeoff  Microburst  Encounter 


Precautions: 

1)  Delay  the  takeoff  if  able. 

2)  Select  the  longest  suitable  runway. 

3)  Perform  takeoff  planning  for  adverse  conditions  as 
prescribed  by  the  NATOPS  Flight  Manual. 

4)  Increase  rotate  airspeed  to  120  knots,  or  the  airspeed 
that  is  attained  with  2000  feet  of  runway  remaining. 

5)  Thoroughly  brief  the  takeoff  procedure,  voice  calls,  and 
windshear  indications  among  the  flight  crew. 
Penetration  of  a  windshear  during  takeoff  will  be 
Indicated  by  a  loss  of  airspeed  or  no  airspeed 
acceleration. 

Procedure; 

1)  Abort  the  takeoff  If  windshear  is  indicated  before 
reaching  91  knots  airspeed. 

2)  Continue  the  takeoff  if  windshear  is  experienced  after  91 
knots  airspeed.  Rotate  the  nose  when: 

a)  120  knots;  or 

b)  2000  feet  of  runway  remain. 

3)  Rotate  to  a  pitch  attitude  of  15o.  DO  NOT  delay  rotating 
the  nose  because  of  low  airspeed. 

4)  Increase  pitch  attitude  toward  20  «>  if  ground  impact  is 
imminent. 

5)  Raise  the  landing  gear  when  above  100  feet  AGL  with  a 
positive  rate  of  climb. 
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V.  RECOMMENDATIONS 


Recommendations  are  broken  into  two  sections.  The  first  set  are 
for  analysis  refinement  and  further  research.  The  second  set  are 
directed  to  the  P-3  operators. 


A.  ANALYSIS  REFINEMENT  AND  CONTINUED  RESEARCH 

There  were  some  weak  areas  of  this  analysis  that  were  pointed 
out  earlier.  Although  valid  results  were  obtained,  refinement  of  the 
windshear  model  and  the  aircraft  equations  of  motion  would  further 
explore  the  edges  of  the  operating  envelopes.  The  following 
recommendations  are  provided  for  the  microburst  windshear  model: 

1)  Vortex  wind  field  models  seem  to  provide  the  closest 
representation  of  a  microburst  windshear,  short  of  applying 
Navier-Stokes  principles.  Schultz's  multiple  vortex  model 
combined  with  source  flows  should  be  combined  and  fitted  to 
other  previously  recorded  microbursts. 

2)  A  more  sophisticated  parameter  estimation  scheme  should  be 
applied  to  enhance  and  accelerate  parameter  fits  of  recorded 
data  to  the  vortex/source  flow  model. 

3)  Realistic  boundary  conditions  to  the  windshear  model  should  be 
defined.  This  would  help  explore  the  entry/exit  aircraft 
response. 

The  following  recommendations  are  for  the  aircraft  performance 
models; 
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1)  The  lift  and  drag  equations  should  be  improved  from  second  to 
fourth  order  equations.  This  change  will  enhance  the 
representation  of  the  aircraft  performance  at  low  airspeed/high 
angle  of  attack. 

2)  An  available  thrust  algorithm  should  be  integrated  into  the 
aircraft  performance  model.  Comparison  of  turbojet,  turbofan, 
turboprop,  and  reciprocating  engine  propulsion  may  lead  to 
significant  effects  on  performance  through  windshear. 

3)  A  state-space  matrix  of  aircraft  dynamics  should  be  integrated 
into  the  performance  algorithm.  Some  aircraft  may  exhibit 
unexpected  dynamics  associated  with  microburst  windshear 
and  escape  maneuvers. 


B.  P-3  OPERATIONS  AND  FLIGHT  CREW  TRAINING 

An  engineering  analysis  provides  the  performance  specifications 
for  penetrating  a  microburst  and  escaping.  However,  certain  issues 
must  be  addressed  by  the  P-3  community  so  as  to  successfully 
employ  the  results: 

1)  Flight  crew  coordination  training  is  essential  for  successful 
employment  of  microburst  escape  procedures.  All  members  of 
the  flight  deck  must  understand  the  teamwork  required  to 
execute  these  procedures.  Therefore,  comprehensive  windshear 
training  should  be  developed.  The  FAA  Windshear  Training  Aid 
and  this  thesis  provides  an  Initial  foundation. 

2)  A  windshear  algorithm  should  be  developed  for  the  P-3  Flight 
Simulators.  There  is  no  effective  inflight  means  to  expose 
flight  crews  to  microburst  effects.  Nor  can  escape  procedure  be 
effectively  practiced.  The  only  proven  instructional  means  are 
with  a  flight  simulator. 

3)  A  brief  dissertation  on  microburst  and  gust  front  windshear 
should  be  provided  in  the  flight  station  NATOPS  manual. 
NATOPS  change  recommendations  reflecting  the  conclusions  of 
this  study  should  be  adopted. 
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APPENDIX  A 


WINDSHEAR  MODEL  ALGORITHMS 


The  following  are  MatLab®  functions  which  model  the  winds 
experienced  at  DFW  on  August  2,  1985  by  two  separate  aircraft. 
AA539  penetrated  the  microburst  at  2500  feet.  DAL191  penetrated 
the  same  microburst  while  attempting  a  landing. 

A.  AA539  WINDSHEAR  MODEL  ALGORITHM 


function  [uX|Uz]"iicrobursl(ze,x) 


X  miTIRL  URLUES  AND  INPUTS 
X  Earth  position 
Xze  •  altitude  (ft) 

Xx  horizontal  position  (ft) 

X  Uortex  Ring  Di eons  ions 
RL-8503.3j 
riL-2001.1; 
kL-431968.8; 

roeL-riL/0.371; 

yL-3350.4; 

LUoralt"3100.6; 

Rs-1701.7; 

ris-323.9; 

ks-57204.9; 

rois"rls/0.371; 

ys*830.9j 

sUoralt*2333.6; 


Xradius  of  ring  fi latent  (ft) 

Xcore  radius  (ft) 

Xuortex  strength  of  the  ring 
Xf  i  laeent(ft'"2/sec) 

Xrelated  core  radius 
XIarge  ring  y  lateral  position  (ft) 
X large  uortex  altitude  RGL  (ft) 
Xradius  of  ring  f I  latent  (ft) 

Xcore  radius  (ft) 

Xuortex  strength  of  the  ring 
Xf  I  latent  (ft''2/sec) 

Xrelated  core  radius 

Xstall  ring  y  lateral  position  (ft) 

Xstall  uortex  altitude  RGL  (ft) 


no 


%  PRinRRV  LRRGE  RING  CRLCULRTIONS 
z*ze-LUonQltj 
if  x»"0, 
if  yL--0, 
gxL*Oj 

uzL*-kL*RL'‘2/(2*(Rr2+z"2)"1.5)j 

zerock“Oj 

end 

eiee 

zerock»999; 

r*sqrt(x''2'^yL^2);  Iradiol  distance  froe  the  z  axis 

rl"sqrt(z^2‘^(r>RL)''2);  Xclosest  distance  to  the  point  of 

X interest  froi  the  ring  fi latent. 
r2»sqrt(z^2+(r+RL)'‘2);  Xfarthest  distance  to  the  point  of 

X interest  froo  the  ring  fiiatent. 
Iaibda"(r2-rl  )/(r2+r1 );  Xscaling  ten 

n1-(r+RL)/r2; 
n2«(r-RL)/rl ; 
sigl*z/r2+z/rl ; 
slg2«z/r2-z/rl ; 

t  of-0 . 75*sqrt (1-1 aibda"2 ) +0 . 25 i 
de  I  tat  •(kL*  I  atbda)/(p  I  Hot) ; 

de  1 1  a2-(0 . 2955«kL«  I  atbda^3)/(p  I  Hoi'^2«8qrt  (1-1  atbda^2) ) ; 
dS I dr«0 . 788*de I ta1 *(n1 -n2)-0 . 394* I atbda*de I t o1 *(nl +02) 

♦de  I  ta2*(  (n1  -n2)- 1  aibda*(nUn2)) ; 

Xderiuative  of  the  streae  function  ert  r 
dS i dz-de I ta1 *(0 . 788*8 1 g2-0 . 394* I atbda*8 i g1) 

♦de I ta2*(8 i g2- I aobda*8 I g1 ) ; 

Xderivatiue  of  the  streat  function  trt  z 
dipL-(1-cxp(-rr2/(0.1*roiL'‘2))); 

Xuelocity  dotping  factor 

uxL*(x/r)*(1/r)*dSidzj  Xhorizontal  velocity  in  x  direction 

uzL”-(1/r)*d5idr;  Xvertical  velocity 

end  Xif 
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X  PRinRRV  SnRLL  RIHG  CRLCULRTIONS 
Z"ze-8Uora 1 1  j 
If  x—0, 
if  y9««0, 
vx8“0; 

uzs»-ks*R8"2/(2*(Rs"2+z'‘2)"l  .5) ; 
end 
el8e 

r''8qrt((x-50)''2+y8''2);  Xradial  di8tance  fpoi  the  z  axle 

rl"8qrt(z''2+(r-R8)^2);  Xcloeeet  dletonce  to  the  point  of 

Xintereet  froo  the  ring  filooent. 
r2*8qpt(z'‘2+(r+R8)'‘2);  Xfoptheet  dietonce  to  the  point  of 

Xintereet  froo  the  ring  filnent. 
laibdo-(p2-pl)/(p2+pl);  Xecoling  ten 

n1-(p+R8)/p2; 
n2-(p-R8)/pl j 
8igl»z/p2+z/pl ; 

8lg2*z/p2-z/p1 ; 

t  oi*0 . 75*8qpt (1-1 aibdo^2)+0 . 25 ; 
de I ta1 "(ke* I aebda)/(p i *toi) ; 

de I ta2-(0 . 2955«k8« I aibdo^3)/(p i «toe^2*8qrt (1-1 aebda^2) ) ; 
dS i dr-O . 788*de I t a1 «(n1 -n2)-0 . 391» I aibda«de I t a1 
•(n1+n2)+delto2*((ni-n2)-latbda*(n1+n2)); 

Xderivotive  of  the  etreao  function  ort  r 
dS i dz-de I ta1 *(0 . 788«8 i g2-0 . 394« I aebda*8 i g1 ) 

‘^^de  i  ta2*(8  i  g2- 1  aibda*8  i  g1 ) ; 

Xderiuative  of  the  streoe  function  ort  z 

if  (pr2/(0.i*poi8*2))>411, 
dip8”1 ; 

else  Xuelocity  daeping  factor 

dips- 1  -exp  ( -r  P2/  ( 0 . 1  •roa8''2 ) ) ; 
end 

ux8*(x/r)*(1/r)*dSidz;  Xhorizontol  uelocity  x  direction 

uz8*-(1/r)*dSldr;  Xuertical  uelocity 

end  Xif 
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iriRGE  LRR6E  RiHG  CRLCULRTIONS 


% 

x»x; 
z*-ze-LUoralt; 
if  x“»0, 
if  yL— 0, 
yxLi"0; 

uzLi-kL*Rr2/(2*(Rr2+z"2)"l  .5); 
end 
elee 

r*eqrt(x'‘2-*-yL^2);  Iradial  distance  froe  the  z  axis 

rl*sqpt(z''2+(r-RL)'‘2);  Sciosest  distance  to  the  point  of 

I interest  froe  the  ring  fi latent. 
r2"sqrt(z'‘2+(r+RL)‘^2);  Ifarthest  distance  to  the  point  of 

X interest  froe  the  ring  fi latent. 
Iaebda*(r2-r1 )/(r2+r1 );  Xscaling  tere 

nl"(r+RL)/r2; 
n2*(r-RL)/rl; 
slgl*z/r2+z/rl ; 
slg2*z/r2-z/r1 j 

t  Oi*0 . 75*sqrt  <1-1  atbda''2  )^0 . 25 ; 
delta1*(l(L*laebda)/(pi*toe); 

de  1 1  a2-(0 . 2955*kL*  I  atbda"3)/(p  I  *tOi'‘2*sqrt  (1-1  aebda''2) ) ; 
dSidr-0.786*delta1*(n1-n2)-0.391*laBbda«delta1 
•(n1 ♦n2)+de I ta2*( (n1 -n2)- I aebda*(n1 +n2) ) ; 

Xderivatiue  of  the  streae  function  ert  r 
dS I dz-de I t a 1 «(0 . 788«s i g2-0 . 391« I atbda«s i g1 ) 

‘•’de  I  ta2*(s  I  g2- 1  aBbda*s  I  g1 ) ; 

Xderivatiue  of  the  streae  function  ert  z 
depL  i  -(I  -exp(-r  r2/(0 . 1  •roiL'‘2) ) ) ; 

Xveiocity  daeping  foctor 

uxLi*(x/r)*(1/r)*dSidz;  Xhorizontai  ueiocity  in  x  direction 

uzLi«-(1/r)*dSidr;  Xuerticai  velocity 

end  Xif 
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X  mRGE  SnRLL  RING  CRLCULRTIONS 
z*-ze“sUoralt j 
if  x»*0, 
if  y8"0, 
vxsi*0; 

yzs i --k8»Rs'‘2/(2*(R8''2+z^2)"l . 5) ; 
end 
el8e 

r*8qrt((x+50)''2+y8'‘2)j  Xradlol  distance  fro®  the  z  axis 

r1»8qrt(z‘'2+(p-R8)''2)j  Xclosest  distance  to  the  point  of 

Xinterest  fro®  the  ring  fila®ent. 
r2*8qrt(z‘^2+(r+R8)'‘2)j  Xfarthest  distonce  to  the  point  of 

Xinterest  fro®  the  ring  fiia®ent. 
Ia®bda“(r2-rl )/(r2+rl );  Xscaling  ter® 

n1“(r+R8)/r2; 
n2-(r-R8)/rl ; 

8lgl-z/r2+z/rl ; 

8ig2*z/r2-z/rl j 

t oi-0 . 75*8qrt (1-1 a®bda"2)+0 . 25 ; 
delta1*>(ks*!aibda)/(pi*toe}; 

de  I  ta2*(0 . 2955*ks*  I  a®bda''3)/(p  I  *tOi'‘2*8qrt  (1-1  a®bda'^2) ) ; 
dS i dr-0 . 788*de I ta1 *(n1 -n2)-0 . 394* I oibda*de I ta1 
•(n1  ♦n2)+de  I  ta2*(  (n1  -n2)- 1  a®bda*(n1  ’•’n2) ) ; 

XderivQtive  of  the  strea®  function  ert  r 
dS i dz*de 1 1  a 1 « ( 0 . 788*8 i g2-0 . 391* I a®bda*s i g 1 ) 

-^de  1 1  a2*  (s  i  g2- 1  a®bda*8  i  g1 ) ; 

Xderivotive  of  the  strea®  function  ert  z 

If  (rr2/(0.1*rois'‘2))>111, 
dap8i«1 ; 

else  Xvelocity  doaping  factor 

daps  1  ■  1  -exp(-r  P2/(0 . 1  *roe8‘'2) ) ; 
end 

vx8i“(x/r)*(1/r)*dSidz;  Xhorizontal  velocity  x  direction 

vz8l“-(1/r)*dSidr;  Xuertical  velocity 

end  Xif 
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X  sunnRTioN 

If  2erock”999, 
dip>dipL*dtps*dipL i *dtp3 i ; 
ux»-(dip*(uxL+uxs+uxL i +yx8 i ) ) j 
y2*dip*(uzL+uz8-yzL i -yz8 i ) ; 
el  88 

yx*-(yxL+yx8+yxL i ♦yx8 i ) ; 
yz«yzL+yz8-yzL i -yz8 i ; 
end 

B.  DAL191  WINDSHEAR  MODEL  ALGORITHM 
function  [yx,yz]*eicrobur8t(ze,x) 

X  INITIRL  URLUES  RNO  INPUTS 


X  Earth  poeition 

Xze  ■  altitude  (ft) 

Xx  horizontal  position  (ft) 

X  Uortex  Ring  Oieensione 

RL-7000; 

Xradius  of  ring  fi latent  (ft) 

riL-2004.1; 

Xcore  radius  (ft) 

kL«13 1968. 8*1; 

Xuortex  strength  of  the  ring 

Xfi  latent  (ft''2/88c) 

roeL*riL/0.371 j 

Xrelated  core  radius 

yL-300; 

XIargs  ring  y  lateral  poeition  (ft) 

LUoralt”3400; 

Xiarge  uortex  altitude  R6L  (ft) 

Lxdi8p*2500; 

Xiarge  ring  displacetent  in  x  (ft) 

R8-f300; 

Xradius  of  ring  fi latent  (ft) 

ri8"323.9; 

Xcore  radius  (ft) 

k8-57204.9*2.3j 

Xuortex  strength  of  the  ring 

Xfi latent  (ft^2/sec) 

roi8»r 18/0.371 ; 

Xrelated  core  radius 

Xstall  ring  y  lateral  position  (ft) 

8Uoralt*‘800; 

Xstall  uortex  altitude  RGL  (ft) 

8Xdisp«300; 

Xstall  ring  displacetent  in  x  (ft) 
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t  PRiriRRV  LRR6E  RING  CRLCULRTIONS 

z«2e-LUorQ I t ; 

If  X--0, 
if  yL«0, 
uxL*0j 

uzL— kL*RL'‘2/(2*(RL"2+z'‘2)''1 .5) ; 
zerock‘0; 
end 
else 

zerock“999; 

r*8qrt(x^2‘>-yL^2)j  Xradial  distance  fnoe  the  z  axis 

r1»sqrt(z'‘2+(r-RL)^2);  Iclosest  distance  to  the  point  of 

X interest  froe  the  ring  fi latent. 
r2*8qrt(z'‘2+(r+RL)'“2);  Xfarthest  distance  to  the  point  of 

Xinterest  froi  the  ring  fi latent. 
Iatbda>(r2-rt)/(r2‘*^r1);  Xscaiing  lert 

nl-(r+RL)/r2; 
n2*(r-RL)/rl; 
slgl*z/r2+z/rl ; 

8lg2*z/r2-z/rl ; 

tot*0 . 75*sqrt  (1-1  atbda''2)-^0 . 25; 
delta1«(kL*latbda)/(pi*toi); 

de  1 1  a2-(0 . 2955«kL«  i  atbda^3)/(p  i  *tor2*8qrt  ( 1  - 1  atbda^2) ) ; 
dSidr-0.788*delta1«(n1-n2)-0.394*latbda«delta1«(n1+n2) 

■^de  I  ta2*(  (n1  -n2)- 1  atbda*(n1  '•^n2) ) ; 

Xderivative  of  the  street  function  ert  r 
dS i dz-de I ta1 «(0 . 788«8 i g2-0 . 394« latbda«8 i g1) 

-^de  I  ta2*(8  i  g2- 1  atbda*s  i  g1 ) ; 

Xderivatiue  of  the  street  function  trt  z 
dtpL«  ( 1  -exp  ( -r  1  ''2/  (0.1  •rotL'‘2 ) ) ) ; 

Xuelocity  datping  factor 

uxL*(x/r)*(l/r)*dSidz;  Xhorizontal  velocity  in  x  direction 

uzL"-(l/r)*dSidr;  Xuertical  velocity 

end  Xif 
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X  PRinRRV  SHRLL  RING  CRLCULRTIONS 
z*ze-8Uoralt; 

if  X--0, 

if  ys—O, 
ux8"0; 

vzs--k3*R8"2/  ( 2*  (  R8'‘2+z*2  ) "  1 . 5 )  ; 
end 
else 

r“8qrt((x-50)'‘2+y8'‘2);  Xrodiol  dl8tance  froi  the  z  oxie 

r1*8qpt(z^2+(r-R8)"'2);  Xcioeeet  dietonce  to  the  point  of 

Xintereet  froe  the  ring  fiioeent. 
r2"sqrt(z''2+(r*Rs)''2)i  Xfarthest  dietonce  to  the  point  of 

Xintereet  froe  the  ring  filoaent. 
laBbda"(r2-r1)/(r2‘»'r1);  Xecoiing  tero 

nJ"(r+R8)/r2; 
n2*(r-R8)/rl; 

8lgl"z/r2+z/rl ; 

8lg2*z/r2-z/r1 j 

toi*0 . 75*8qrt  (1-1  aibda''2)+0 . 25; 
deltaf*(k8*(atbda)/(pi*toe); 

de I t a2-(0 . 2955«k8« I aebda^3)/(p i »to8^2«8qrt (1-1 aebda^2) ) ; 
dSidr-0.788«delta1«(n1-n2)-0.394«lanbda«deltQl 
*(nl  ♦n2)+de  I  tQ2*(  (nl  -n2)- 1  aibda*(n1  <^02) ) ; 

Xderivotive  of  the  etreon  function  ort  r 
dS i dz-de I tol «(0 . 788«8 i g2-0 . 394* I oibda«8 i g1 ) 

<*’de  I  to2*(8  i  g2- 1  oibdo’s  i  g1 ) ; 

XderiMOtiue  of  the  etreon  function  nrt  z 

if  (rr2/(0.1*rone'‘2))>M1, 
dnpe*1; 

elee  Xvelocity  donping  factor 

dipe-l  -exp(-r  1  ''2/(0 . 1  *ron8'‘2) ) ; 
end 

ux8”(x/r)*(1/r)*dSidz;  Xhorlzontol  uelocity  x  direction 

uze*-(1/r)*dSidr;  Xuerticol  uelocity 

end  Xif 
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inflGE  LRRGE  RING  CRLCULrTIONS 


t 

X"Xj 

z“-ze-LUoralt; 
if  x"0, 
if  yL--0, 
uxLi»0; 

vzLi-kL*RL"2/(2*(Rr2+z"2)"l  .5); 
end 
else 

r*8qrt(x"'2-^yL''2);  Iradial  distance  froe  the  z  axis 

r1"sqrt(z'^2+(r-RL)''2)i  Iciosest  distance  to  the  point  of 

Kinterest  froe  the  ring  filaeent. 
r2>sqrt(z^2‘>’(r‘»^RL)'^2);  ^farthest  distance  to  the  point  of 

linterest  froi  the  ring  fi latent. 
Iaebda"(r2-r1)/(r2+r1);  tscaling  ten 

nl"(r+flL)/r2; 
n2*(r-RL)/rl; 

8lgi«z/r2+z/rl; 

slg2*z/r2-z/rl; 

t  oi«0 . 75*8qrt ( I - 1 aibda^2)+0 . 25  j 
de I  tat ■(kL* I aebda)/(p I *toe) ; 

de  I  ta2-(0 . 2955*kL*  I  aibda"3)/(p  I  »toe'‘2*8qrt  (1-1  aebdo'‘2) ) ; 
dSidr-0.788«delta1«(n1-n2)-0.394»loebda«delta1 
•(nl  ♦n2)+de 1 1 a2*( (n1 -n2)- 1 aebda*(n1 +02) ) ; 

Ideri votive  of  the  streoo  function  trt  r 
dS i dz-de I  tat  *(0 . 788*8 1 g2-0 . 394* I oibdo«s i g1 ) 

<^de  I  ta2*(8 1  g2- 1  aibda*s  I  g1 ) ; 

Iderivotive  of  the  streoo  function  ort  z 
dopL  I  ■  ( 1  -exp  ( -r  1  '‘2/  (0.1  •rofL''2 ) ) ) ; 

Ivelocity  dooping  foctor 

vxLi«(x/r)*(l/r)*dSidz;  Xhorizontol  velocity  in  x  direction 

vzLi*-(l/r)*dSidr;  Xverticoi  velocity 

end  Xif 


118 


X  mflGE  SHALL  RING  CALCULATIONS 

z«-3e-8Uoralt; 

If  x-0, 

If  y8«0, 
ux8i*0; 

uzs I «-k8*R8'‘2/(2*(Rs"2+z''2)'‘l  .5) ; 
end 
el  86 

r*8qrt((x‘<-50)''2+y8''2);  Xradial  dietonce  froe  the  z  axle 

rl*8qpt(z'‘2+(r-R8)^2)j  Xcloeeet  dietance  to  the  point  of 

Xintereet  froe  the  ring  filoient. 
r2*8qrt(z'‘2+(r+R8)*2);  Xfortheet  dietance  to  the  point  of 

Xintereet  froi  the  ring  filonent. 
Iaibda"(r2-rl)/(r2+rl);  Xecollng  ten 

nl"(r+R8)/r2; 
n2*(r-R8)/rl; 

8lgl*z/r2+z/rl ; 

8lg2*z/r2-z/r1; 

toi*0 . 75*8qrt (1-1 oibda^2)+0 . 25; 
deltQl«(k8*iatbda)/(p(*toi); 

de I t a2-(0 . 2955«k8« I aibda^3)/(p I «toe^2»8qrt (1-1 anbda^2) ) ; 
dS I dr-0 . 788«de I t o1 «(n1 -n2)-0 . 391« I anbdQ«de I ta1 
•(n1  ♦n2)+de  I  ta2*(  (n1  -n2)- 1  anbda*(n1  <^n2) ) ; 

Xderi votive  of  the  etreon  function  ort  r 
dS  i  dz-de  1 1 0 1  *  ( 0 . 788*8  i  g2-0 . 394*  I  aibda*8  i  g  1 ) 

-^de  I  to2*(8  i  g2- 1  anbda*8  i  g1 ) ; 

Xderi votive  of  the  etreoe  function  ert  z 

if  (rr2/(0.1*roe8"2))>141, 
dep8i*1 ; 

else  Xvelocity  doeping  factor 

dope  I  ■!  -exp(-r P2/(0 . 1  *roe8''2) ) ; 
end 

vxei"(x/r)*(1/r)*dSidz;  Xhorizontol  velocity  x  direction 

vz8i*-(1/r)*dSidr;  Xverticoi  veiocity 

end  Xif 
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X  SOURCE  PROFILE  FOR  UZ  CORRECTION 

Salt**10000;  Xaource  altitude 

gagaz*2772000;  Xsource  strength/2T 

gagax" 1355396;  Xsource  8trength/2T  far 

uzsource^O.O; 

uxsource*0.0; 

if  ze>800, 

vzeource*-!  5*(LUora  I  t-ze)/(  (LUora  I  t-ze)'"2+x*2)‘'0 .5; 

else 

Radlu8l-(x'‘2+(Salt-ze)"2)"0.5; 

Xradiue  fros  prisary  source  to  the  flight  pt 
Red  I  u82-(x"2+(Se  1 1  ♦ze)"2)''0 . 5; 

Xradius  fros  loage  source  to  the  flight  pt 
uz80urce*gagaz*(  (Sa  I  t**-ze)/Rad  i  u82^2-(Sa  I  t-ze)/Rad  I U8P2) ; 
end 

ux80urce"gagax*x*  ( 1  /Rad  I  us  1  ''2*  1  /Rad  i  us2''2 ) ; 

X  sunnflTioN 
If  zerock—999, 
dtp«dspL*dsps*dspLi*dipsi ; 
if  x<-2000, 

ux«(dop*(uxL'*'Uxs<*’UxL  I  '•’uxs  i ) ) ; 

else 

vx*vxsource; 

end 

uz*-(dap*( vzL+vzs-vzL i -uzs i )-uzsource) ; 
else 

ux»-( wxL+uxs+uxL I ♦uxs I ) j 
uz*uzL^uzs-uzL I -uzs I ; 
end 
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APPENDIX  B 


PARAMETER  SENSITIVITY  ALGORITHM 


The  following  are  MatLab®  m-file  and  function  algorithms  used  to 
determine  the  parameters  for  the  windshear  model  that 
approximates  the  DAL191  recorded  winds. 


A.  M-FILE  DRIVER 

XParaieter  sensitiuity  prograi  for  the  DRL191  eind  eodel. 

clear 

cic 

{Recorded  flight  data  froe  DRL191. 

X[ref  nuiber,  altitude(ft  RGL),  diet  fre  RB  center(ft),  x  find 
(ft/e),  y  find  (ft/e)]. 


DRL-Il 

987  -8642 

2  10; 

2 

975 

-8361  -1 

11 

3 

962 

-8080  -3 

11 

4 

948 

-7799  -4 

11 

5 

935 

-7518  -6 

14 

6 

922 

-7237  -9 

15 

7 

909 

-6956  -12 

12 

8 

894 

-6676  -12 

3 

9 

875 

-6396  -9 

6 

10 

856 

-6117  -11 

9 

11 

842 

-5838  -16 

4 

12 

623 

-5561  -25 

-18 

13 

607 

-5287  -36 

-14 

M 

790 

-5016  -39 

-13 

15 

777 

-4749  -29 

-20 

16 

764 

-4486  -26 

-17 

17 

755 

-4229  -22 

-17 
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18 

717 

-3977 

-29 

-17; 

19 

712 

-3730 

-20 

-11; 

20 

738 

-3187 

-17 

-23j 

21 

730 

-3219 

-12 

-30j 

22 

711 

-3013 

-8 

-36; 

23 

698 

-2778 

0 

-31; 

21 

680 

-2512 

3 

-32; 

25 

662 

-2306 

0 

-31; 

26 

617 

-2069 

-8 

-31; 

27 

635 

-1831 

-18 

-21; 

28 

627 

-1593 

-20 

-26; 

29 

620 

-1355 

-10 

-36; 

30 

608 

-1117 

-8 

-10; 

31 

592 

-880 

18 

-1; 

32 

585 

-639 

32 

19; 

33 

582 

-391 

19 

16; 

31 

580 

-139 

18 

21; 

35 

573 

120 

19 

-3; 

36 

518 

385 

38 

-21; 

37 

511 

656 

51 

-19; 

38 

162 

936 

65 

9; 

39 

111 

1228 

66 

9; 

10 

351 

1531 

72 

-13; 

11 

277 

1815 

71 

-12; 

12 

193 

2169 

72 

-33; 

13 

119 

2501 

76 

-2; 

11 

70 

2851 

78 

15; 

15 

12 

3211 

61 

17; 

16 

20 

3579 

81 

1; 

17 

0 

3951 

85 

2] 

XParateter  aatrix: 


[2300 

2500 

2800; 

XLorge  ring  x  disploceient 

1250 

1300 

1350; 

XSbqII  ring  core  radius 

57201.9*2 

57201.9*2.3 

57201.9*2.5; 

XSeaii  ring  uortex  strength 

-100 

-300 

-200; 

XLorge  ring  y  posit 

300 

500 

700; 

XSiaii  ring  X  dispioceient 

-150 

1 

150; 

XSeaii  ring  y  dispiocseent 
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1386000*1.8  1386000*2  1386000*2.2;  ISource  z  strength 

193628*6.5  193628*7  193628*7.53;  ISource  x  strength 


uxerr*0.0;wzerr"0.0;tn*0;bestris»999;  Unit  lol  Izot  ion 

i-1; 

z-1; 

for  k*1:3, 
hooe 

I 

I 

z 

k 

best res 

for  l■1:3,  for  ■■1:3,  for  n*1:3, 
for  0*1 ;3,  for  p*1 :3, 
uxerr*0 . 0 ; uzerr"0 . 0 ; 
for  c"l!47, 

[Ux, Uh]*po i croburst (OflKc, 2) « DflLCc, 3) , Ps(  1 « i ) , Ps(2, z) , 
Ps(3,k),Ps(4J),Po(5,i),Ps(6,n),Ps(7,o),Ps(8,p)); 
vxerr*uxerr*(0flL(c,1)-Ux)^2; 
uzerr-uzerr+ (OflL  (c « 5  )-Uh }  '"2} 
end 

rt8vx*(uxerr/40)''0 .5; 
risuz*  ( uzerr /40 )  '‘0 . 5 ; 
net  ot  *  ( ( wxerr  ♦  vzerr  )/80)''0.5; 

■n*in<^1 ; 

Pn(in,;HP8(M),Ps(2,z),Ps(3,k),P«(4,l), 

Ps(5,i),P8(6,n),Ps(7,o),Ps(8,p),rtsux,ri8vz,ri8tot]; 
if  bestris>ri8toti 
BESTP-Pn(in,:); 
bestris^restot; 
save  BEST 
end 

end,  end,  end,  end,  end,  end, 
eoue  porooeterell 

cleor  ox  eh  reeux  rietot  en  Pfl  BESTP 
vxerr*0.0;uzerr*0.0;en*0;  linitiol izot ion 
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i-lj 

z-2; 

for  k*l :3, 

hole 

« 

I 

z 

k 

beet res 

for  l■l:3,  for  ■■1:3,  for  n“l:3, 
for  o“l :3,  for  p»l :3, 
uxerr“0 . 0 ; uzerr-0 . 0 ; 
for  c*l :47, 

[Ux,Uh]"paicroburst(DfiL(c,2),DRL(c,3),Pe(t ,  i),Pe(2,z), 
Pe(3,k),Pe(1,l),Pe(5,i),Pe(6,n),Pe(?,o),Pe(6,p)); 
uxerr -uxerr  ♦  ( DflL  ( c ,  4 )  -Ux )  ''2 ; 
uzerr-uzerr  ♦  ( DRL  ( c ,  5 )  -Uh )  ^2 ; 
end 

rtevx"  ( uxerr/40 )  '‘0.5; 
reeuz*  ( uzerr /40 )  '‘0.5; 
net  ot  ■  ( ( wxerr+yzerr )  /OO )  '‘0 . 5 ; 
antin'*’ t; 

Pn(in,!)-[Pe(1,l),Pe(2,z),Pe(3,k),Pe(4J), 

Pe(5, i) ,Pe(6, n) , Pe(7,o) ,Pe(8,p) ,rieyx, raeuz, netot ] ; 
if  beetrie>rietot, 

BESTP-Pn(an,:); 
beetne^netot; 
eove  BEST 
end 

end,  end,  end,  end,  end,  end, 
eoue  paraaetere12 

clear  ex  ih  raeux  rietot  in  Pfl  BE5TP 

yxerr*0.0;yzerr*0.0;an«0;  Unit  iol  izat  ion 

i-1; 

z-3; 

for  k*l:3, 
hole 

i 
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z 

k 

bestras 

for  I-IjS,  for  t'hS,  for  n*l:3, 
for  0*1 ;3,  for  p*1 :3, 
uxerr*0 . 0 ; yzerr“0 . 0  j 
for  c*1 :17, 

[Ux^Uh]*poicrobur8t(DRL(c,2},0RL(c^3),P8(1 , i ),P8(2,z), 
P8(3,k),P8(4J),P8(5,t),P8(6,n),P8(7,o),P8(8,p)); 
uxerr“uxerr+  ( DRL  ( c ,  4 )  -Ux )  '‘2 ; 
wzerr"uzerr+ (0RL(c,5)-Uh)^2; 
end 

rt8ux*  ( oxer  r/40 )  ''0 . 5  j 
ri8UZ”  ( uzerr /10 )  '‘0 . 5 ; 
ri8t  ot  ■  ( ( oxer  r+uzer  r )  /80 )  '‘0 . 5 ; 

■n«en+1 ; 

PfKin, !)*[P8(1, l),p8(2,z),P8(3,k),P8(1, I), 

p8(5«e)4p8(6in)|P8(7,o),P8(8,p),re8vx«re8uz,re8tot]; 
if  be8tro8>re8tot> 

BESTP-PnCen,:); 
be8tre8«ro8tot; 
save  BEST 
end 

end,  end,  end,  end,  end,  end, 
save  paraeetersIB 

clear  ex  eh  rievx  restot  in  PH  BESTP 

vxerr*0.0;vzerr*0.0;en"0;  f initial ization 
1-2; 
z-1; 

for  k*1 :3, 
hoee 
I 

z 

k 

best res 

for  1*1:3,  for  i*l:3,  for  n*l;3, 
for  0*1 :3,  for  p*l ;3, 
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uxerr»0 . 0 ; uzerr“0 . 0 ; 
for  c"1 ;17, 

[UX|Uh]*piicroburst(DRL(cj2)|DRL(c«3)>P9(1 , i )^P8(2,z}^ 
P8(3,k),P8(4, l),P8(5,i),P9(6,n),P8(7,o),P8(8,p)); 
wxerr-uxerr  +  ( DRL  ( c ,  1)  -Ux )  ^^2 ; 
vzerr“yzern+(DRL(cj5)-Uh)'‘2; 
end 

ni8vx» ( vxerr/10 ) ^0 . 5  j 
Pi8uz»(vzerr/10)''0.5j 
ni8t  ot  ■  ( ( yxerr+yzerr )  /80 )  ''0 . 5 ; 

■n“tn+l I 

Pn(in,:)-[P8(l,i),P8(2,z),P8(3,k),P8(1J), 

P8(5,i),P8(6,n),P8(7,o),P8(8,p),ri8yx,ri8yz,r88tot]; 
if  be8tri8>rB8tot| 

BESTP-Pn(in,:); 
be8tri8"ra8tot ; 

8aye  BEST 
end 

endj  end,  end,  end,  end,  end, 
eaye  paraaeter821 

clear  ex  eh  rasux  raetoi  an  Ptl  6ESTP 

yxerr-0.0;yzerr*0.0;an«0;  Unit  iol  izat  ion 
1-2; 
z-2; 

for  k*l ;3, 
hoae 
I 

z 

k 

be8tra8 

for  l■1;3,  for  a"l:3,  for  n*l!3, 
for  o"l :3,  for  p"1 ;3, 
yxerr-0 . 0 ; yzerr"0 . 0 ; 
for  c"l 547, 

[ Ux , Uh ] "pa i cr obur 8t ( DRL ( c , 2 ) , DRL ( c , 3 ) , P8 ( 1 , i ) , P8 ( 2 , z ) , 
P8(3,k),P8(4,l),P8(5,a),P8(6,n),P8(7,o),P8(B,p)); 
vxerr -y  xerr  ♦  ( DRL  ( c ,  4 )  -Ux )  ^*2  ; 
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uzerr"uzepp+  ( DflL  ( c ,  5 )  -Uh )  '*2 ; 
end 

piewx- ( yxepp/40 ) ^0 . 5  j 
peeyz*(yzepp/40)'‘0 .5; 

Piet  ot  ■  ( ( wxepp+uzepp )  /80 )  '‘0 . 5 ; 

■n"in+1 ; 

Pn(in,  :)-[P8(l, l),P8(2,z),P8(3,k),P8(1, 1), 

P8(5,i),P8(6,n),P8(7,0),P8(8,p),Pi8UX,Pi8UZ,Pe8t0t]j 
if  be8tpe8>pe8tot, 

BESTP-Pn(in,!); 

be8tpe8*pe8tot; 

8aue  BEST 
end 

end,  end,  end,  end,  end,  end, 
eaue  papaeeteP822 

cleap  fx  eh  peeux  peetot  en  Pft  BESTP 

uxepp"0.0;uzepp*0.0;en«0j  finitiai izotion 
1-2; 
z-3; 

fop  k*l 53, 
hoee 

I 

z 

k 

best pee 

fop  l■l:3,  fop  i"1:3,  fop  n*l;3, 
fop  o"1;3,  fop  p*1;3, 

Mxepp*0 . 0 ; vzepp*0 . 0 ; 
fop  c*1:47, 

[Ux,Uh]-peicpobuP8t(DRL(c,2),DflL(c,3),P8(1,i),P8(2,z), 

P8(3,k),P8(1,  l),P8(5,e),P8U,n),P8(7,o),P8(8,p))j 
vxepp-Mxepp+(DflL(c,  4)-Ux)''2j 
wzepp-uzepp*  ( DHL  ( c ,  5 )  -Uh )  ^'2 ; 
end 

pieux" ( vxepp/40 ) ^0 . 5 ; 
pe8uz*(uzepp/40)''0 . 5 ; 
pest  ot  ■  ( ( vxepp+vzepp )  /BO )  '*0 . 5 ; 
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■n*in+l ; 

Pn(tn,:)-[P8(l,i),P8(2,z),P8(3,k),p8(1J), 

P8(5,t),P8(6,n),P8(7,o),P8(8,p),ri8yx,ri8yz,ri8tot]; 
if  be8tri8>ni8tot, 

BESTP-Pn(in,:)j 

be8tri8-ri8tot; 

8aye  BEST 
end 

end,  end,  end,  end,  end,  end, 
eoue  paraeeter823 

clear  ex  eh  reeux  reetot  in  PH  6ESTP 

yxerr*0.0;Mzerr"0.0jen>0;  Xinit ioi izat ion 

i-3; 

z-1; 

for  k*l :3, 
hoee 

I 

z 

k 

beet r 18 

for  l■l!3,  for  ■■1:3,  for  n«1:3, 
for  0*1 :3,  for  p«l :3, 
wxerr*0 . 0 ; vzerr«0 .0; 
for  c»l ;17, 

[ Ux , Uh ] -pi I croburet ( DflL ( c , 2 ) , DflL ( c , 3 ) , P8 ( 1 , I) , Pe ( 2 , z ) , 
P8(3,k),P8(1,l),P8(5,«),P8(6,n),P8(7,o),P8(B,p)); 
vxerr-uxerr ♦ ( DHL ( c , 1 ) -Ux ) ^2 ; 
yzerr"wzerr+  ( DflL  ( c ,  5 )  -Uh )  ''2 ; 
end 

ri8VX"  ( vxerr /40 )  ''0 . 5 ; 
ri8vz*  ( uzerr /10 )  ''0 . 5 ; 
r  Bet  ot  ■  ( ( uxerr+yzerr )  /BO )  ''0 . 5 ; 

■n-in-^l ; 

Pn(tn,;)-[P8(1,i),p8(2,z),P8(3,k),P8(1,l), 

P8(5, i) , P8(6, n) , P8(7,o) ,P8(B,p) ,rB8vx, reeuz, rietot ] ; 
if  be8trB8>ri8tot, 

BESTP-Pn(in,:); 
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bestri8«ri8tot; 

8aue  BEST 
end 

end,  end,  end,  end,  end,  end, 
eaue  paronetereSI 

clear  ex  eh  rteux  reetot  in  Pfl  BESTP 

uxerr*0.0;vzerr»0.0jin“0j  Xinit iai izat ion 

1-3; 

z-2; 

for  k-1 :3, 
hoee 

I 

z 

k 

beetree 

for  1-1 ;3,  for  ■-1;3,  for  n-l;3, 
for  0-1 !3,  for  p-1 s3, 
uxerr-0 . 0 ; uzerr-0 . 0 ; 
forc-l:17, 

[Ux,Uh]-peicrobur8t(0fiL(c,2),DflL(c,3),P8(l,l),P8(2,z), 

P8(3,k),P8(1,l),P8(5,i),P8(6,n),p8(7,o),P8(8,p)); 
yxerr-wxerr+ ( DflL  ( c ,  4 )  -Ux )  ^^2 ; 
uzerr-wzerr+ ( DflL  ( c ,  5 )  -Uh )  '‘2 ; 
end 

ri8yx-( vxerr/40)^0 .5; 
ri8vz-(  vzerr/10)''0 . 5 ; 
net  ot  -  ( ( uxerr+yzerr )  /BO )  ''0 . 5 ; 

■n-en+l ; 

Pn(in,:)-[P8(l,l),P8(2,z),P8(3,k),P8(4,l), 

PeCS, ■) , PeCB, n) , P8(7, o) , PeCB, p) , rieox, reeuz, reetot ] ; 
if  be8tre8>re8tot, 

BESTP-Pn(en,!); 
beetree-reetot; 
eaue  BEST 
end 

end,  end,  end,  end,  end,  end, 
eaue  paraeetere32 
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clear  ex  eh  resux  restot  en  PH  BESTP 


uxerr*0.0;uzerr"0.0;en*0;  Xinit iai izat ion 

i-3; 

z-3; 

for  k«l :3, 
hooe 
I 

z 

k 

best res 

for  l-hS,  for  ■■1:3,  for  n*1:3, 
for  o"1 :3,  for  p*l :3, 
yxerr-O.Ojwzerr-O.Oj 
for  c*1 :17, 

[Ux,Uh]"piicroburst(DflL(c,2),DfiL(c,3),Ps(1 , i),Ps(2,z), 
Ps(3,k),p8(4,l),Ps(5,i),Ps(6,n),Ps(7,o),P8(8,p)); 
wxerr  "uxer  r+  ( DflL  ( c ,  1 )  -Ux )  '‘2 ; 

Mzer  r -yzerr*  ( ORL  ( c ,  5 )  -Uh )  '^2 ; 
end 

ns  vx*  ( vxerr /40 )  ''0 . 5 ; 
risuz*  ( vzerr/40 )  ''0 . 5 ; 
rest  ot  ■  ( ( uxerr+vzerr )  /80 )  ''0 . 5 ; 

■n^en-t-l ; 

Pn(in,:)-[P8(l,l),P8(2,z),P8(3,k),P8(1,l), 

P8(5,i),P8(6,n),P8(7,o),Pe(8,p),ri8vx,ri8yz, restot]; 
if  be8tri8>ri8tot, 

8ESTP-Pn(tn,:); 
be8tri8*ri8tot; 
save  BEST 
end 

end,  end,  end,  end,  end,  end, 

eaue  paratetersSS 

clear 

load  BEST 

uxerr*0 . 0 ; uzerr»0 . 0 ; 
for  l^hlT, 
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J 


hole 

I 

X-0flL(l,3); 

[Ux,Uh]"piicroburst(OflL( i ,2),DflL( i ,3),BESTP(1)>6ESTP(2), 

BESTP(3),BESTP(1),BESTP(5),BESTP(6),BESTP(7),BESTP(8)); 
■odeKi,  :)>[i,UX|Uh]; 
gxerr»yxerr+(DflL( i ,1)-Ux)^2j 
uzerr-vzerr+CDflK I ,5)-Uh)^2; 
end 

nievx*(yxepn/40)''0 .5; 
ri8yz“(wzerp/40)'‘0.5; 

Pist  ot  ■  ( ( yxepp+yzepp )  /BO )  ^'0 . 5 ; 
ansiep^Oj 

•hile  oneiep'^OOO, 
diepCThe  hopizontal  RHS  is:  ') 
disp(pieyx) 

diepCThe  yeptical  RfIS  ie:') 
diep(pieyz) 

diepCThe  totai  eyetei  RHS  ie:') 
diep(pietot) 

diepCThe  foiioiing  piote  ape  auai table: ’) 
diepC  1*hopizontal  linde') 
diepC  2"yeptical  linde') 
diepC  999-end') 
aneiep-inputCpiat  nuibep  '); 
if  aneiop”1, 
piot(DRL(:,3),DRL(:,1)) 
hold 

plot(DRL(:,3),iodel(:,2)/-') 

pauee 

hold 

e  lee  if  aneiep— 2j 
piot(DRL(:,3),DRL(:,5)) 
hoid 

piot(DRL(:,3),iodei(:,3)C-') 
pauee 
hoid 
end I  end 
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B.  WINDSHEAR  FUNCTION 


function  [vXjVzl-poicrobupotfze^xAxdiop^Ro^ko^yLjSxdisp^RLAi-) 
X  INITIRL  URLUES  RND  INPUTS 


X  Earth  position 

Xze  ■  altitude  (ft) 

Xx  horizontal  position  (ft) 

X  Uortsx  Ring  Di sens  ions 

XRL 

Xrodius  of  ring  fiiaoent  (ft) 

rlL-2004.1; 

Xcore  rod i us  (ft) 

XkL 

Xuortex  strength  of  the  ring 
Xfilossnt  (^""2/800) 

rosL>riL/0.371 ; 

Xrslatsd  core  radius 

XyL 

XIargs  ring  y  lateriai  position  (ft) 

LUoraltOlOO; 

X large  uortsx  altitude  RGL  (ft) 

XLxd i sp 

XIarge  ring  displacesent  in  x  (ft) 

XRs 

Xradlus  of  ring  filasent  (ft) 

rls-323.9; 

Xcore  radius  (ft) 

Xks 

Xuortex  strength  of  the  ring 
Xflloient  (^^^2/888) 

ross«rls/0.371; 

Xrslatsd  core  radius 

ys-830. 9/830.9; 

Xssall  ring  y  lotsriol  position  (ft) 

sUoralt*800; 

XsmII  uortsx  altitude  RGL  (ft) 

Xsxdisp 

Xsull  ring  displacesent  in  x  (ft) 

X  PRINRRV  LRR6E  RING  CRLCULRTIONS 

z*ze>LUora I t ; 

zerock-999j 

r"8qrt  ( (x-Lxd  I  »p)''2+yL'‘2) ; 

r1-oqrt(z^2*(r-RL)^2); 

r2-sqrt(z^2+(r+RL)^2); 

Iaibda”(r2-p1 )/(r2^r1 ); 

n1"(r+RL)/r2; 

n2-(r-RL)/r1j 
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8lgl*z/r2+z/rl; 

8lg2*z/r2-z/r1 ; 

t  Of  *0 . 75*8qr  t  (1-1  afbda'‘2 ) +0 . 25 ; 
de I  to t ‘(kL* I afbda)/(p i *tof) ; 

de  I  ta2-(0 . 2955*kL*  I  afbda''3)/(p  I  *lof''2*8qrt  (1-1  afbda"2) ) ; 
dS i dr»0 . 788*de 1 1  a 1 ♦ ( n 1 -n2 ) -0 . 391* I a»bdo*de 1 1  a 1 
*(nl+n2)+della2*((nl-n2)-laibda*(nl+n2))j 
dS i dz*de I tol *(0 . 788*8 i g2-0 . 394* I afbda*8 I gl ) 

♦de 1 1  a2* ( 8 1 g2- 1 atbda*8 i g 1 ) ; 
dfpL*(  1  -exp(-rP2/(0 . 1  *rofL^2) ) ) ; 
uxL*(x/r)*(l/r)*dSldz; 
yzL*-(l/r)*dSldr; 

X  PRItlRRV  SnRLL  RING  CRLCULRTIONS 
z"ze-8Uora I t ; 

r"8qrl  ( (x-8xd  1 8p)'‘2+y8'‘2) ; 
r  1  ■8qrt  ( z'‘2+ ( r-R8 ) '‘2 ) ; 
r2-8qrt (z^2+(r+R8)^2) ; 

Iafbda*(r2-rl )/(r2+rl ); 

nl*(r+R8)/r2; 

n2«(r-R8)/rl; 

8lgl*z/r2+z/rl; 

8lg2"z/r2-z/rli 

tof *0 . 75*8qrt ( I - i afbda^2)^0 . 25; 
de I tal ■(k8* I afbda)/(p i *tof ) ; 

de  1 1 a2-(0 . 2955*k8*  I  qibda''3)/(p  I  •toi"2*8qrl  (1-1  afbda"2) ) ; 
dS 1 dr-0 . 78e*de I t a 1 *(nl -n2)-0 . 394* I afbda*de I to 1 
*(nl+n2)+delta2*((nl-n2)-  Iafbda*(nl+n2)); 
dS i dz-de I tol «(0 . 786*8 i g2-0 . 394* I afbda*8 i gl ) 

♦de I ta2*(8 i g2- I afbda*8 i gl ) ; 

If  (rr2/(0.1*rof8"2))>444, 
dfp8*l ; 
elee 

dope*  1  -exp  ( -r  1  ^2/  (0.1 5*rof8'*2 ) ) ; 
end 

ux8«(x/r)*( 1 /r)*dS I dz; 
uz8*-(l/r)*dSldr; 
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X  inflGE  LRR6E  RING  CRLCULRTIOHS 
z"-ze-LUoraIt; 

r-8qrt  ( (x-Lxd  i  8p)''2+yL'‘2) ; 
r1-8qrt(z"2+(r-RL)'‘2); 
r2*8qrt  ( z'‘2+ ( r+RL )  ^^2 )  j 
Iaibda»(r2-rl )/(r2+rl ); 
n1»(r+RL)/r2j 
n2»(r-RL)/r1 ; 

8jgl*z/r2+z/rl ; 

8ig2"z/r2-z/r1 ; 

t  ot«0 . 75*8qrt  (1-1  atbda'‘2 ) +0 . 25 ; 
delta1"(kL*la8bda)/(pi*to8}; 

de  1 1  a2*(0 . 2955*kL*  I  a8bda''3)/(p  i  *tot'‘2*8qrt  (1-1  aibda^2) ) ; 

dSidp-0.788*deltal«(nl-n2)-0.391*laBbda*deltal*(nl-^n2)-i'delta2*((nl-n 

2)-la8bda*(nl+n2))j 

dS i dz-de 1 1  a 1  * ( 0 . 788*8 i g2-0 . 394* I aBbda*8 i g 1 ) 

-^de  I  ta2*(8 1  g2- 1  aibda*8  i  gl ) ; 
dipL I -( 1 -exp(-r r2/(0 . 1 *rofL"2) ) ) ; 
uxL i ■(x/r)*( 1 /r)*dS I dz; 
vzLI»-(l/r)*dSldr; 

X  INRGE  SRRLL  RING  CRLCULRTIOHS 

z--ze-8Uoralt; 

r»8qrt  ( (x-8xd  1 8p)'‘2+y8''2) ; 

rl  "sgrt  (z''2+(r-R8)''2)  j 

r2-8qpt  (z''2+(p+R8)'‘2)  ; 

Ia8bda*(p2-pl )/(p2+p1 )j 
nl«(p+R8)/p2; 
n2*(p-R8)/pl ; 

8igl«z/p2+z/pl; 

8lg2*z/p2-z/pl  j 

toi*0 . 75*8qpl  (1-1  a8bdQ^2)’»0 . 25 ; 
de I tal ■(k8* I aibda)/(p I *toi) ; 

de  I  ta2«(0 . 2955*k8*  I  aibda^3)/(p  i  *toi'‘2*8qpl  (1-1  a8bda''2) ) ; 
dS i dp-0 . 788*de I tal *(nl -n2)-0 . 394* I aabda*de I tal 
*(nl+n2)+delta2*((nl-n2)-  Iatbda*(nl't’n2)); 
dSldz-deltal*(0.788*8lg2-0.394*latbda*8lgl) 
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-^de  I  tQ2*(s  i  g2- 1  aibda*8  i  g1 ) ; 

If  (rr2/(0.1*rot8"2))>M4, 
dipsi*1 j 
else 

dips  i  ■!  -exp(-r  P2/(0 . 1 5*roi8'‘2) ) ; 
end 

ux8i"(x/r)*(l/r)*dSidz; 

uz8l»-(l/r)*dSidr; 

X  SOURCE  PROFILE  FOR  U2  CORRECTION 

Salt"10000;  Xsource  altitude 

Xgagaz  Xsource  8trength/2f 

Xgagax  Xsource  8trength/2i  for 

uz8ource>0.0; 

ux8ource"0.0; 

If  ze>800, 

uz80urce--1 5*(LUora  I  t-ze)/(  (LUora  I  t-ze)''2+x''2)'‘0 . 5 ; 
else 

Radlus1-(x'‘2+(Salt-zer2)"0.5; 

Rad  I  u82-(x^2+(Sa  I  t+ze)^2)'‘0 .  s 

uzsource-gagaz* ( (Sa 1 1  +ze ) /Rad  I u82^2-(Sa 1 1 -ze ) /Rad  I  us  r2 ) ; 
end 

wxsource«gagax*x*(  1  /Rad  I  us  1  ^2+1  /Rad  I  us2'‘2) ; 

X  sunnRTioN 
dip*dipL*dips*dipL I *dap8 I ; 

If  x<-2000, 

ux"(dap*(uxL'^uxs't’UxL  I  ^uxs  I ) ) ; 
else 

ux”ux80urce; 

end 

uz*-  ( dip*  ( uzL-^uzs- vzL  I  - vzs  I )  -uzsource ) ; 
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APPENDIX  C 


LANDING  APPROACH  ALGORITHMS 


Contained  in  this  appendix  are  the  MatLab®  programs  utilized  in 
calculating  the  aircraft  response  given  a  windshear  and  a  particular 
escape  maneuver.  They  are  listed  in  three  sections.  The  first  section 
contains  the  driver  routine  and  the  functions  required  to  solve  the 
differential  equations.  The  second  and  third  sections  contain  the 
specific  aircraft  parameter  functions  and  escape  maneuver  functions 
respectively.  The  windshear  function  is  listed  in  Appendix  A. 

A.  DRIVER  AND  DIFFERENTIAL  EQUATION  SOLVER 

|««»»«*««»«»»»»»««»«»««»»«»*«« 

%  Descent  Profile  Rode  I  * 

I***************************** 

%  This  prograe  till  calculate  the  inertial  flight  path  angle  for 
%  a  given  aircraft.  It  utilizes  the  function  '■icroburst'  ehich 
X  node  Is  the  oind  shear  encountered  during  the  DFU  Del  to  Reel  dent. 

X  This  prograe  uses  energy  height  theory  to  calculate  inertial  space 
position. 

clear 

cic 

foreat  bank 
foreat  coepact 

I - 

XINITIRLIZRTIDH  RND  EHU I RONRENTRL  CONSTRNTS 
X  The  fol losing  initialize  constants  and  variables  for  the  eodel. 
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XRIRCRRFT  CONSTRNTS 

X  The  fol losing  constants  are  for  the  particular  aircraft  being 
studied. 

d i ep( '  LRND I NG  RPPRORCH  RODEL ' ) 

dispC  ’) 

aneeer^O.O; 

dispCType  aircraft: ') 
dispC  1.  P3  at  89,500ib') 

dispC  2.  P3  at  114, 5001b') 

diepC  3.  P3  at  89,5001b  3  engine') 
dispC  4.  Lion') 

dispC  5.  T44') 

ansoer*inputC  Enter  desired  aircraft  by  nusber:  '); 

If  anseer—l, 

[S,K,C0o,Tiax,Ut,CLo,CLalpha,alphaoax,q,U1,T,theta]>P3LtRpp; 

XP3  •  89, 500 lbs,  app  config. 

el  seif  ansoer”2, 

[  S ,  1C ,  COo ,  Tiax ,  Ut ,  CLo ,  CLa  I  pha ,  a  I  phaoax ,  q ,  U 1 ,  T ,  t  het  a  ]  ■P3HuRpp ; 

XP3  •  11 4, 500 lbs,  app  config. 

el  seif  ansoer”3, 

[S,)C,C0o,Toax,Ut,CLo,CLalpha,aiphaoax,q,U1,T,theta]*P3Eout; 

XP3  oith  one  engine  out. 

el  seif  ansoer"4, 

[S,IC,CDo,Tiax,Ut,CLo,CLalpha,alphaMx,q,U1,T,theta]*L1011Rpp;XL1011 
else  if  onsier**5, 

[ S , K , COo , Toax , Ut , CLo , CLa I pha , a I phaeox , q , U 1 , T , t het a ] ■T44Rpp ; XT44 
end 

XS«plane  fort  surface  area  in  ^^"2 
XK^coeff icient  for  CD  calculation 
XCDo*coeff icient  of  drag 
XTiax>tax  thrust  auailabie  in  ft>lbs 
XUt ■aircraft  eeight  in  lbs 
n*Ut/32.174;  X  mss  in  slugs 
XCLo*zero  ROR  lift  coefficient 
XCLalpho"! ift  curve  slope 
XalphaMx*stol  I  buffet  alpha  •/  approach  flaps 
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Xq^noiinal  pitch  rate  of  Sdeg/sec 

XU1 ■approach  equiualent  airepeed  in  ft/s 

XT"initiai  thrust  in  ibf 

I - 


{sssppogpQied  initial  values**************************************** 
iters*30;  Xnuaber  of  iterations 

deitat«1.0;  Xtiee  increeent  in  sec 

FIBcenter*  10300;  Xdi stance  110  center  is 

froe  runeay 

xstart>-500;  Xstarting  distance  froe 

no  center 

I******************************************************************* 


hstart^(n0center-xstart)*tan(0.052);  Xstarts  on  a  3deg 

gl ideslope 

dispC  Initial  Input  Section') 
dispC  ') 

dispCThe  fo I  losing  uaiues  are  preprograssed : ' ) 
dispC  0.  The  approach  is  based  upon  a  3deg  g I Ideslope. ' ) 
dispC  b.  The  starting  distance  (distance  froe  nO  center)  is:') 
dispCxstart) 

dispC  c.  The  nO  center  is  this  froe  the  end  of  the  runeay:') 
dispfnOcenter) 

diepC  d.  The  storting  altitude  is:') 
dlsp(hstart) 

dispC  e.  The  tiee  step  in  sec:') 
disp(deltat) 

dispC  f.  The  nunber  of  iterations  is:') 
disp( iters) 
dispC  ') 


TU*Tsax/Ut;  Xthrust/oeight  ratio 

UL*Ut/S;  Xeing  loading 

dispCThe  foi  toeing  uaiues  eere  suppi  ied/catcutatsd  fros  aircraft 
constants') 

dispC  a.  Uref  in  ft/s  (1 .35Uso'*^5kts  for  P3  and  1.3Uso  for  other 

aircraft:  ’) 

disp(UI) 

dispC  b.  Thrust  to  eeight  is:') 
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disp(TU) 

dispC  c.  Uing  loading  isi') 

disp(UL) 

dIopC  •) 

>aneuMer"0.0; 

diop( 'Select  the  escape  naneuMer:') 


dispC  1. 
f indshear') 

Stay  on 

gl  ide 

siope^ 

go 

■ i seed  at 

200ft 

RGL;  NO 

dispC  2. 
e i ndshear ' ) 

Stay  on 

gl  ide 

slope. 

go 

■ i seed  at 

200ft 

RGL;UITH 

diepC  3,  Escape  oith  constant  airspeed') 
dispC  4.  Escape  elth  constant  altitude') 
dispC  5.  Escape  elth  constant  theta') 
dispC  6.  Escape  elth  constant  alpha') 
dispC  7.  Escape  elth  tax  alpha') 

dispC  8.  Eaulate  081191  final  fee  seconds  using  recorded  theta') 
dispC  9.  Eaulate  DflL191  final  fee  seconds  using  recorded  theta  and 
finds') 

■aneuuer- 1 nput ( '  Ent  er  ae I ect i on  nuiber :  ' ) ; 

if  oaneuuer^-Si 

dispCUhat  value  for  escape  theta?') 
thetaaia*inputC  Enter  in  radians:  '); 
el  seif  aaneuver^O, 
dispCUhat  value  for  escape  alpha?') 
alphaaia-inputC  Enter  in  radians:  '); 
el seif  aaneuver**?, 
alphaaia*aiphaaax; 
elseif  ■aneuver>*8, 
hstart*hstart+50; 
end 

x1-xstart; 
hi*hstart; 

Ug-Ui; 

gaaai*-.052; 
rad  x-3deg) 


linitiai ization 
linitial ization 

Iground  speed  initialization  in  ft/e 
Sinitial  inertial  flight  path  angle  in 
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rof ■0.002377; 
trlgger*0.0; 
check 1*0.0; 
check2"0.0; 


lair  density  in  siuge/ft^'S  •  S.L. 
lutility  check  nueber 
futility  check  nueber 
lutility  check  nueber 


RLGORITHn 

The  actual  algorithe  is  broken  into  seueral  sealler  diuisions  and 
unct i ons : 

starting  uatues  -  as  it  ieplies 

loop  -  beginning  of  Iterating  process 

oindshear  inputs  -  uses  the  function  'eicroburst'  to  giue 

point  X  ond  z  axis  eind  speeds 

predictor  -  first  guess  at  the  solution  of  the  gouerning  OEs 
corrector  -  corrected  solutions  to  the  gouerning  DEs  and 
final  point  ualues. 

output  file  eanageeent  -  as  it  Ieplies 

aircraft  control  -  uses  eultiple  control  functions  for 

control  inputs  to  theta,  thrust,  and  ROR 
final  output  -  as  it  ieplies 


letarting  ualues: 
pt-1; 

Q-0.5*roi»Ur2; 

CL-Ut/(g*S); 

CD-CDo+K*CL"2; 

a  I pha* ( CL-CLo ) /CLa I pha ; 

gaeaa"theta-alpha; 

Es-Ur2/(2*32.174)+hl; 

Esdot"-12.32; 


Icounter 

Idynaeic  pressure 

Icoefficient  of  lift 

Icoefficient  of  drag 

IROR  required  for  1g  flight 

laireass  flight  path  angle 

Ispecific  energy  in  ft 

It  lee  rate  change  in  specific  energy 


y(l)-xl 
y(2)-hl 
y(5)-U1 
y(6)-gaeaa; 


ly(t)*x  position  in  ft 
ly(2)«aircraft  altitude  in  ft 
ly(5)*aircraft  airspeed  in  ft/s 
ly(6)«aircraft  aireaes  fit.  path  (rad) 


Rircraftd,  :)*[x1,  hi,  U1, -12. 36,  theta,  alpha,  1,T]; 

I nert i a  I ( 1 , : )-[Ug, 0, gaeaa, gaea 1 , 0, 0, 1 , Es, Esdot ] ; 

I  starting  output  files 

X - 


MO 


Xioop  control 
Spresent  pt  in  space 


X I  OOP : 
cic 

•hile  theta*"999, 
pt*pt+l; 

% - 

Xiind  shear  inputs: 
x2*Ug*deltat’*-x1 ;  Xopproxlnote  horizontal  disploceoent 

h2"tan(gaBal)*(x2-xi)*h1 }  Xopproxltote  vertical  disptoceient 

if  laneuaer**!, 

Ux 1 1 -0 ; Ux2 1 -0 ; Ux 1 2-0 ; Ux22-0 ; Uh 1 1 -0 ; Uh2 1 -0 ; Uh 1 2-0 ; Uh22-0 ; 
el  seif  soneuaer— 9j 

[Ux1]«Uh1t]-recsinds(pt);  Xrecorded  finds  rros  DRL191 

[  Ux22 ,  Uh22  ] -rece  i  nds  ( pt 1 )  i 

Ux12-(Ux1HUx22)/2j 

Ux21-Ux12j 

Uht2-(Uh11«Uh22)/2: 

Uh21-Uh]2; 

else 

[Uxlt,Uht1]-iicroburst(h1,x1);  XThese  four  colls  to  the  function 
CUx2f ,Uh2l]-oicroburst(h2,xl);  Xnlcroburst  finds  the  point  wind 
[Ux12,Uhl2]-iicroburst(h1|X2);  Xsheor  for  deteroination  of  the 
[Ux22,Uh22]-oicroburst(h2,x2);  Xoueroge  block  find  shear. 

XThe  aboue  is  used  in  conjunction 
Xother  calls  ohen  no  tind  shear  is 
Xdesired. 
end 

UindX(pt)-Ux1t;  XUind  shear  eatrices 

UindH(pt)-Uh11j 

dUxdx- ((Ux12-Ux11)*( Ux22~Ux2 1 ) ) / ( 2* ( x2'X 1));XThese  f oreu I  as 
dUxdh-((Ux22-Ux12)^(Ux21-Ux11))/(2*(h2-h1));Xcalculate  the 
dUhdx-((Uh12-Uh1 1 )>(Uh22-Uh21 ))/(2»(x2-x1 ));Xdi f ferent iai  change 
dUhdh-((Uh22-Uh12)^(Uh21-Uh11))/(2»(h2>h1));Xln  elnd  shear. 

y(3)«Ux11;  Xy(3)-eind  shear  in  the  x  dir.  ft/s 

y(1)-Uhlij  Xy(1)-eind  shear  in  the  z  dir.  ft/s 


HI 


Xpredicton  and  intenediate  ualues 
ydot -ypr i tes ( y , dUxdx , dUxdh , dUhdx , dUhdh , t  het  a , T , S « H , Q , CL , CD , 
alpha^alphonax); 

Xcall  to  the  OE  function 

for  n»l :6, 

y1 (n)*y(n)+deltat*ydot(n);  XEuler  first  step  fed  predictor 
yhalf(n)«y(n)-^deltat/2*ydot(n);  XEuler  half  step  fed  predictor 
end 

alpha”theta-‘y(6);  XROR 

CL>CLo‘*-aipha*CLaipha;  Xcoefficient  of  lift 

CD"CDo+IC*CL''2;  Xcoefficient  of  drag 

Q*0.5*roe*yhaif(5)''2;  Xdynaeic  pressure 


Xcorrector  and  nee  uaiue  foreulation 
ydot -ypr i eesfyho I f ^  dUxdx^  dUxdh, dUhdx, dUhdh^  thet a>  T>  5, n, Q, CL^  CD^ 
aipha^aiphaeax); 

for  n«l :6, 

y2(n)"yha I f (n)-i’de 1 1 at /2*ydot (n) ; XR i chardson  extrapo lotion 
y(n)"2*y2(n)-y1(n);  Xcorrector  scheoe. 

end 

xl«y(I);  Xnee  x  posit  in  ft 

h1"y(2);  Xnee  altitude  in  ft  RGL 

Ug*y(5)*cos(y(6))<^y(3);  Xground  speed  in  ft/sec 

R0C-y(5)*8in(y(6))‘*’y(1);  XRate  of  Clieb  in  ft/sec 

gaiai>atan(ROC/Ug);  Xnee  a i mass  FP  angle  in  rad 

E8"y(5)'^2/6i.3id*y(2);  Xnee  specific  energy 

E8dot«y(5)*ydot(5)/32.174‘^ydot(2);Xnee  trc  of  specific  energy 

aipha"theta-y(6);  Xnee  ROR 

CL"CLo-«>alpha*CLalpha;  Xnee  coefficient  of  lift 

CD*CDo'^K*CL''2;  Xnee  coefficient  of  drag 

Q*0.5*roe*y(5)^2;  Xnee  dynaeic  pressure 
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lout put  file  eanogeoent: 

fllrcroftCpt, ;)«[y(l), y(2),y(5),ydot (5), theta, alpha, pt,T]; 

I nert i a  I (pt , : )*[Ug, ROC, y(6) , gaea i , y(3) , y(1) , pt , Es, Esdot ] ; 

hooe 

dispC  pt  X  h  U 

airspeed  accel  ') 

opl -[pt , R i rcraf t (pt , 1 ) , R i rcraft (pt , 2) , R I rcraf t (pt , 3) , 

Rircpaft(pt,1)]; 

disp(opl) 

dispC  pt  theta  alpha  gaeaa  gaoai') 

op2"[pt, R i rcraf t(pt, 5), Ri rcraft (pt, 6), Inert iai (pt,3), 

InertiaKpt, 

di8p(op2) 

diepC  pt  Ug  ROC  Thrust') 

op3"[pt , I nert I  a  I (pt , 1) , I nert I  a  I (pt , 2) , R I rcraf t (pt , 8) ] ; 
dl8p(op3) 

dispC  pt  Es  Esdot') 

op5>[pt,  InertiaKpt,  8),  InertiaKpt, 9)]; 
di8p(op5) 

dispC  pt  Ux  Uh') 

op1«[pt,  InertiaKpt,  5),  InertiaKpt, 6)1; 
dl8p(op4) 

I - 


{aircraft  control: 

{functions  called  for  theta,  thrust 
land  alpha  control. 

if  oaneuuer  <•  2, 

If  y(2)<200, 
trlgger-1; 

[theta,T]"nRcontroKq,deltat,Teax,theta,T,alpha,alphaiax,y); 
elself  trigger  1, 


M3 


[thetaJ]*GPcontroi(qideltat|ROC  Jiax,theta,T,alpha,alphaiax^U1  > 

h8larl,x8lart,y); 

el self  trigger  ■■  I, 

[theta J]-nRcontrol(q,deltat,Teax,theta J;alpha,alphaeax^y); 
end 

el  seif  ■aneuuer<8> 
if  y(5)  <  (U1-33.33), 
trigger*!; 
end 
else 

If  y(2)<-1004, 

[ t het a , T , check 1 ] "DAL t het a ( Teax , t het a , T , U 1 , y , check 1 ) ; 
else 

[theta|T]-GPcontrol(q,deltat,R0CJeax,thetaJ,alpha,aiphaaax,U1, 

hstartjxstartjy); 

end 

end 

if  aoneuaer  >*3, 

If  trigger  ■■  1, 
if  aaneuuer  ■*  3| 

[ t het a , T I  check 1 ] -CRSesc ( q , de 1 1 at Jaax , t het a J , a  I pha , a  I phaaax , 

Ul,y, check!) 

el  seif  aaneuuer  ■■  1, 

[ t het a , T , check ! ] "CRLTesc ( q , de 1 1 at , Taax , t het a , T , a  I pha I  a  I phaaax , 

y, check l,ydot) 

el  seif  aaneuver  ■■  5, 

[thetaj T]*CTHesc(q, del  tat, Taax, theta, T, alpha, alphaaax,thetaaia) 
el  seif  aaneuver  ■■  6, 

[theta, T]*CROResc(q, de I  tat, Taax, theta, T, alpha, alphaaax,alphaa la) 
el  seif  aaneuver  ■■  7, 

[theta, T]*nRXosc(q, del  tat, Taax, theta, T, alpha, aiphaa la) 
end 
end 
end 

if  pt» iters,  Xnuaber  of  iterations 

theta-999; 


\44 


l9top  at  ground  inpact 


elseif  y(2)<0, 
theta-999; 
else 
end 

end  lohile  algorithn  loop 

% - 

Xoutput  files: 

aneeer- i nput ( ' Do  you  eant  a  flight  path  plot?  1-VES  '); 

If  anseer--!, 

plot(fllrcraft(:,1),fllrcraft(;,2)) 

pause 

anseer-O.O; 

end 

anseer- I nput ( 'Do  you  eant  Xt  h,  U,  and  airspeed  accel?  1-VES  '); 

If  anseer— 1, 

dispC  pt  X  h  U  airspeed  accel  ') 

opl-[fllrcraft(:,7),fllrcraft(:J),filrcraft(:,2),fllrcraft(;  ,3), 
fllrcraft(:,1)]; 
dispCopI) 

else 

end 

aneeer- I nput ('Do  you  eant  theta,  alpha,  gaeaa,  and  gaeoi?  I-VE5  '); 
if  aneeer— I, 

diepC  pt  theta  alpha  ganoa  gaeai') 

op2-[Rircraft(:,7),Rircraft(:,5),Bircraft(:,6), Inert ial(:,3), 
lnertial(:,1)]; 
disp(op2) 
else 
end 

oneeer- input ('Do  you  eant  Ug,  ROC,  ond  Thrust?  1-VES  '); 

If  aneeer— I, 

diepC  pt  Ug  ROC  Thrust') 

op3-[Rircraft(:,7), lnertial(:,l),lnertial(:,2),Rircraft(:,8)]; 
disp(op3) 
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else 

end 

anseer* i nput ( ' Do  you  eont  Es  and  Esdot?  1>VES  '); 
if  onseer**!, 

dispC  pt  Es  Esdot') 

op5”[R i rcraft ( : , 7) , I nert i a l( : , 8) J  nert i a  I ( : « 9) ] ; 

dispCopS) 

else 

end 

anseer” I nput ( ' Do  you  eant  Ux  and  Uh?  1"VES  '); 
if  anseer*"lj 

dispC  pt  Ux  Uh'} 

op1*[Rircraft(:  j7)Jnert  ial(:i5),  Inert  ial(:,6)]; 

disp(opl) 

else 

end 
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function  [ydot]*ypriie3(y,dUxdx,dUxdh,dUhdX|dUhdh, theta, 

T,S,n,Q,CL,CD,alpha,aiphaoax) 

X  This  function  calculates  the  y'  solutions  to  the 

X  coupled  systeo  equations: 

X  y(1 ) '*xdot«U*cos(gaiaa)-«-Ux 

X  y(2) '=hdot=U*sin(gaiaa)+Uh 

X  y(3)  '>‘Uxdot*’dUx/dx*xdot*dUx/dh*hdot 

X  y(1) •»Uhdot=dUh/dx*xdot*dUh/dh*hdot 

X  y (5) ‘ »Udot *(T*cos(a I pho) )/H-D/l1-g*s i n(gataa)“Uxdot *cos(gaiaa) 

X  -Uhdot*sin(gaeaa) 

X  y(6) '■T*sln(alpha)/(l1*U)+L/(n*U)-g*cos(gaiaa)/U 

X  -^Uxdot*s  i  n(gaeaa)/U-Uhdot*cos(gaoaa)/U 


ydot(l)“y(5)*cos(y(6))+y(3); 
ydot(2)«y(5)*sln(y(6))+y(1); 
ydot (3)*dUxdx*ydot ( 1 )+dUxdh*ydot (2) ; 
ydot (1)*dUhdx*ydot ( 1 )+dUhdh*ydot (2) ; 
if  alpha<alphaiax, 

ydot  (5)-T/f1*cos(thet  o-y  (6)  )-Q*S*CD/n-32 . 1 74*s  I  n(y(6) ) 
-ydot(3)*cos(y(6))-ydot(4)*sin(y(6)); 
ydot (6)-T/(n*y (5) )*s i n(t het o-y (6) )+Q»S*CL/(n*y(5) ) 

-32.171/y(5)*cos(y(6))+ydot(3)*sln(y(6))/y(5)- 

ydot (1)/y(5)*cos(y(6) ) j 
else 

ydot(5)"T/n*cos(alphaoax) -Q*S»CD/n-32 . I 74*s I n ( y ( 6 ) ) 

-ydot (3)*cos(y(6) )-ydot (1)*  sin(y(6) ) ; 
ydot(6)-T/(n*y(5))*sin(alphaMx)*0*S»CL/(n*y(5))-32.l74/y(5) 

*cos(y(6))<»ydot(3)*  sin(y(6))/y(5)-ydot(4)/y(5)*cos(y(6)); 
end 
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B.  AIRCRAFT  PERFORMANCE  PARAMETER  FUNCTIONS 


fund  ion[SjlC,CDo,TiaX|Ut  jCLo, 
“P3Ltflpp() 

XP3  at  89, 500 1 bs^  app  flaps, 
S-1300; 

K-O. 05041; 

C0o"0.0567; 

T»ax«33400; 

Ut-89500; 

CLo*0.600; 

CLaipha*5.73j 

alphaiax"0.244; 

q*0.0873; 

U1-236; 

T-5000; 

theta--0.013; 


CLalpha,alphaiax,q,U1 ,T, theta] 
gear  doen 

Xplane  fore  surface  area  in  ft ‘^2 
Xcoefficient  for  CO  calculation 
Xcoefficient  of  drag 
Xeax  thrust  auailabie  In  ft-lbs 
Xaircraft  eeight  in  lbs 
Xzero  OOP  lift  coefficient 
Xlift  curve  slope 

Xstall  buffet  alpha  e/  approach  flaps 
Xnoeinal  pitch  rate  of  5deg/sec 
Xapproach  equivalent  airepeed  in  ft/s 
Xinitial  thrust  in  Ibf 
Xinitial  eater  line  deck  ongle  in  rad 


f unct i on [ S , K , COo , Teax , Ut , CLo , CLa I pha , 0 I phoeox , q , 01 , T , t het a] 
■P3HvRppO 

XP3  at  11 4, 000 1  be,  app  flaps,  gear  doen 


S-1300; 

K-O. 05041; 

COo-0.0567; 

Teax-33400; 

Ut-1 14500; 
CLo-0.800; 
CLalpha-5.73; 
alphaBax-0.244; 
q-0.0873; 
01-262; 

T-8000; 

theta--0.013; 


Xpione  fore  surface  area  in  ft‘'2 
Xcoefficient  for  CD  caiculation 
Xcoefficient  of  drag 
Xeax  thruet  avaiiable  in  ft- I be 
Xaircraft  eeight  in  I be 
Xzero  ROR  lift  coefficient 
Xlift  curve  elope 

Xetoll  buffet  olpha  e/  approach  flaps 
Xnoeinal  pitch  rate  of  5deg/eec 
Xapproach  equivalent  airepeed  in  ft/e 
Xinitial  thruet  in  ibf 
Xinitial  eater  line  deck  angle  in  rad 
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fund  ion[S,ICjCDo,Tiax, 
-LlOllflppO 
XL1011  at  362>000lb8, 
S-1578; 

K-0.059; 

CDo”0.108; 

Tiox* 126000; 

Ut-362000; 

CLo-0.532; 

CLaipha*4.96; 

alphatax*0.3I4; 

q*0.0873; 

U1-227; 

T-41500; 

theta"0.118; 


Ut^CLOjCLalpha^alphaiaXiq^UI theta] 

Rpp  flaps,  gear  dotn 

Xplane  fort  surface  area  in  ft ^2 
Xcoefficient  for  CD  calculation 
Xcoefficisnt  of  drag 
Xiax  thrust  auailable  in  ft-lbs 
Xaircraft  oeight  in  lbs 
Xzero  ROfl  lift  coefficient 
Xlift  curve  slope 

Xstali  earning  alpha  •/  approach  flaps 

XClOdeg) 

Xnooinai  pitch  rate  of  Sdeg/sec 
^approach  equivalent  airspeed  in  ft/s 

X(136kt8) 

Xinitiol  thrust  in  Ibf 

XInitial  eater  line  deck  angle  in  rad 


function[S,K,CDo, 
■TlIRppO 
XT11  at  8280 lbs, 
S-210; 

K-0.0503; 

CDo-0.120; 

Taax>3023; 

Ut-8280; 

CLo"0.587; 

CLaipha*6.24; 

alphaeax-0.241; 

q*0.0873; 

UI-203; 

T-1250; 

theta>-0.02; 


Teax,Ut,CLo,CLalpha,alphaeax,q,U1,T, theta] 

app  flaps,  gear  doen 

Xplane  fore  surface  area  in  ft'‘2 
Xcoefficient  for  CD  calculation 
Xcoefficient  of  drog 
Xeax  thrust  avaiiabie  in  ft-lbs 
Xaircraft  eeight  in  ibs 
Xzero  ROfl  lift  coefficient 
Xlift  curve  slope 

Xstali  buffet  alpho  e/  approach  flaps 
Xnoeinol  pitch  rate  of  Sdeg/sec 
Xapprooch  equivolent  airspeed  in  ft/s 
Xinitiol  thrust  in  Ibf 
Xinitiol  eater  line  deck  angle  in  rad 
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fund  ion[S,K|CDO|Tiax 
•P3Eoul() 

XP3  d  89, 500 lbs I  app 
S-1300; 

K-0. 05041; 

C0o*0.0630; 

Tiax-33400*(3/4); 

Ut-89500; 

CLo"0.600; 

CLalpha”5.73; 

alphanax*0.244; 

q*0.0873; 

Ul-236; 

T-5000; 

theta“-0.013; 


Ut|CLo,CLalpha,aiphaiax,q,U1,T,theta] 

flaps,  gear  doen,  one  engine  out. 

Xplane  fore  surface  area  in  ft ^2 
Xcoefficient  for  CO  calculation 
Xcoefficlent  of  drag  (inc  for  eng.  out) 
Xeax  thrust  auaiiabie  in  ft-lbs 
Xaircraft  eeight  in  lbs 
Xzero  flOfl  lift  coefficient 
Xlift  curue  slope 

Xstall  buffet  alpha  e/  approach  flaps 
Xnoeinal  pitch  rate  of  5deg/sec 
Xapproach  equivalent  airspeed  in  ft/s 
Xinitial  thrust  in  Ibf 
Xinitial  eoter  line  deck  angle  in  rad 
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C.  FLI6HT  PATH  CONTROL  AND  ESCAPE  MANEUVER  FUNCTIONS 


f unci i on [ t het Qout , Tout ] ® GPcont ro I ( q , de I t at , ROC , Tiox , t het a , T , 

a  I pha , a  I phoaax , U 1 , hot  art , xot  art  ^  y } 
XThio  control  function  tries  to  oa into  in  o  3deg  glide  path. 
Tout*!; 

thetaout*theta; 

•antedh*helort-(abe(xelart-y( 1 ) )*0 . 052H) ; 


hdi ff*y(2)-eantedh; 

Udlff«y(5)-Ul; 

X  theta  input: 
if  alpha<0| 
t  het  aout  *1 hot  a-a I pha ; 
el  self  alpha<alphaeaxj 
if  ROOOj 

t het aout *1 het a-q*de I  tat ; 
e lee  if  R0C<-25.0j 

t  het  aout  *1  het  a‘*'q*de  1 1  at ; 
eleeif  hdiff>100, 

t  het  aout  *t  het  a*q/2*de 1 1  at ; 
eleeif  hdiff<-100, 

t  het  aout  *1  het  a‘^q/2*de  1 1  at ; 
end 

if  t het aout >0.521 I 
t het aout *0.521; 
eleeif  thetaout<-0.175, 
theloout*-0.175; 
end 
else 

t  het  aout  *1 het  a- ( a  I pho-a I phaeax ) ; 
end 

XThruet  input: 

If  Udlff>0.5, 

Tout  *1-0 . 20*T ■ax*de 1 1 at ; 

If  Tout<0, 

Toul*0; 

else 


Idee i red  altitude  for  3deg  6P 
Xalt  difference  (neg  ehen  ioe) 
Xairepeed  diff  (neg  ehen  eloe) 

X(1) check  for  negatiue  ROR 

X(2)check  for  tax  auailabie  ROR 
X(3)check  for  ROC  >  Ofpe 

Xif  so  pitch  den  q*tiBe 
X(1)check  for  ROC  <  1500fpB 
Xif  eo  pitch  up  q*tiBe 
X(5)check  for  alt  >100  diff 

Xif  eo  pitch  den  1/2*q*tiBe 
X(6)check  for  alt  <100  diff 

Xif  eo  pitch  up  1/2*q*tiBe 

Xnoe  check  for  theta  lieit  of 
X  -lOdeg  to  +30  deg 


Xcheck  speed  diff 
Xif  diff  >  5kts  above  app  speed 
Xreduce  thruet  by  20X  of  Teax 
Xbut  not  ieee  than  0. 
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end 

eteeif  Udiff<-6.5, 

Toiit”T+0.25*Teax*deltat;  Xif  diff  <  Skte  beloe  app  speed 

if  Tout>Teax,  Xinc  thrust  by  25X  of  Tsax 

Tout«Teax;  Xbut  not  eore  than  Teax. 

else 
end 


funct ionCthetaout^ Tout ]*nRcontroi(q, del  tat Jeox, theta, T, alpha, 

alphaeax,y) 

XThis  control  function  initiates  a  norsal  eissed  approach. 

X  Hax  poeer  is  added  and  theta  increased  to  lOdeg.  hax  flOR 
Xis  applied  ehen  required. 

Tout"T; 

thetaout "theta; 

X  theta  input: 

if  alpha<alphaeax,  X(1)check  for  tax  auailable  ROR 

If  thetaout>«0.175,  X(2)check  for  lOdeg  deck  angle 

thetaout"0.175; 

else 

thetaout"theta-^q/2*deltat;  X(3)lf  theta  not  lOdeg,  pitch  to 

lOdeg 

end 

else 

thetaout "theta-Ca I pha-a I phaeax) ; 
end 

XThrust  Input: 
if  Tout<Teax, 

Tout"T+Tiax/2*deltat ; 
if  Tout>TMx, 

Tout"Teax; 

end 

else 

end 
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f  unct  i  on  [  t  het  Qout  I  Tout ,  check  1  ] -CRSedc  ( q ,  de  1 1  at  Jux ,  t  het  a  I T , 

0 1 pho , 0 i phonox , U 1 , y « check 1 ) 

Xflircroft  eecope  noneuver  consisting  of  constant  airspeed, 
thetaout "theta; 
if  check 
if  y(5)<U1, 
thetaout'O; 
eise 

thetaout"q*deitat; 
check1"1 ; 
end 
eise 

if  y(5)>-(Ul+8.1), 
t  het  aout  "t  het  a‘^q*de  i  t  at ; 
elseif  y(5)<-(U1-8.4) 
t  het  aout "t  het  a-q*de i t  at ; 
end 
end 

if  aipha>aiphaeaX| 
t het aout "theta- (a i pho-a I phaeax) ; 
end 

if  T<Tmx, 

T  out "T +1 ■ax*0 . 2*de i t  at ; 
else 

Tout"Teax; 

end 
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f unct i on [ t het Qout , Tout , check 1 ] >CflLTesc ( q , de I t at « Tiax j t het a , T , 

0 1 pho 1 0 1 phoeox I y ^ check 1 « ydot ) 

Iflircroft  control  ueing  constant  altitude  escape  eaneuaer. 
thetoout "theta; 
if  check 1»"0; 

check1"y(2); 

end 

if  y(2)<check1-20| 
if  ydot(2)<"0| 
thetaout"theta‘*^q*deitat ; 
end 

el  seif  y(2)>checkl‘*-20| 
if  ydot(2)>-0, 
thetaout"theta-q*de I  tot ; 
end 
end 

if  aipha>aiphaeaX| 
t  het  oout "t  het  0- ( 0 1 pho-o I phoeox ) ; 
end 

if  T<Teaxj 

Tout«T+T»ax*0 . 2*de I  tot ; 
else 

Tout"Tiax; 

end 
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fund i onEthetaout , Tout l"CTHesc(q,de I tot , Toox, theta, T, 

a  I pha, a  I phaeox , t hetaa i ■) 
XConetont  theta  escape  eaneuuer. 
if  theta<thetaaia, 
t  het  aout  "t  het  a-f^q’de  1 1  at ; 
else 

thetaout"thetaaie; 

end 

if  alpha>alphaeax, 
t  het  aout -t  het  a- ( a  I pha-a I phaeax ) ; 
end 

If  KTmx, 

T  out "T +T eax’0 . 2*de 1 1  at ; 
else 

Tout*Teax; 

end 


funct lonEthet aout, Tout ]"CROflesc(q,deltat,Tiax, theta, T, 

a I pha , a I phaeax , a  I phaa I n ) 

XConstant  alpha  escape  eaneuuer. 
if  alpha<alphaaie, 
t  hetaout  ■theta’i’Ca  I  phaa  I  e-a  I  pha) ; 
if  thetaout>(theta'^q*deltat), 
t  het  aout  "t  het  a'''q*de  1 1  at ; 
end 

else 

t het aout d het a- (a  I pha-a I phaa I e) ; 
if  thetaout<(theta-q*deltat), 
t het aout d het a-q*de I  tot ; 
end,  end 

if  aipha>alphaeox, 
t  het  aout "t  het  a- ( a  I pho-o i phaeax ) ; 
end 

if  T<Teax, 

Tout*T+Teax*0 . 2*de I  tat ; 
else,  Tout"Teax;  end 
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f  unct  i  on [  t  het  Qout « Tout  ] "nflXesc  ( q ,  de  1 1  at ,  Tiax ,  t  het  Q ,  T I  a  I  pho  j  a  I  phaa  i  ■ ) 
Xnaxiaui  alpha  escape  aaneuuer. 
if  alpha<alphaaia« 
t  het  aout  "t  het  a-^  ( a  I  phaa  i  a-a  I  pha } ; 

If  thetaout>(theta-*-q*deltat}| 
t  het  aout "t  het  a+q*de 1 1  at ; 
end 

else 

thetaout*theta-(a I pho-a I phaa is) ; 
if  thetaout<(theta-q*deltat)> 
t  het  aout "t  het  a-q*de 1 1  at ; 
end 
end 

If  T<Taax, 

Tout-T+Taax*0 . 2*de I  tat ; 
else 

Tout"Taax; 

end 
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funct  ion[thetaout4Tout|Check1out]**DflLtheta(Tiax,theta  J,U1  check! ) 
I  This  control  function  tries  to  eoulote  DRL191 
Xlost  18  sec  of  flight, 
if  check Xstorting  point 
check! ■32; 
end 

X  theta  input : (theta, Xpooer) 

DRL!9!-[0.0698!3!7007!; 

0.06263!8530  7!;  0.051! 0520687!; 

0.017! 238898  72;  0.017! 23889873;  0.038397213573; 

0.03!1!59265  7!;  0.02268928027!;  0.038397213569; 

0.051!052068  69;  0.051! 05206869;  0.051! 05206870; 

0.0785398!63  68;  0. ! 239! 8376867;  0. 132615023! 67; 

0.1623156201  67;  0.185001900767;  0.1 9! 9862! 7767; 

0.2007128639  67;  0.200712863969;  0.200712863971; 

0.2231021112  72;  0.253072711585;  0.267035375589; 

0.2710166925  92;  0.271016692592;  0.267035375595; 

0.2600510585  92;  0.237361778292;  0.237361778292; 

0.2530727115  87;  0.253072711587;  0.223102111292; 

0.1553313031  93;  0.062631853098;  0.031115926596; 

0.0226892802  99;  -0.005235987100;  -0.068067810100; 

-0.111862327  100;  -0.130899693100;  -0.020913951100; 

0.0698131700  100;  0.0925021503100;  0.0311159265100; 

-0.005235987  100;  0.0511052068100;  0.0311159265100]; 
thetaout  *081191  (check!  *^1 , 1 ); 

Tout-DRL191(checkl«l,2)/100*Tiax; 
check  lout  ■check  I**-!; 
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APPENDIX  D 


TAKEOFF  ALGORITHMS 


Contained  In  this  appendix  are  the  MatLab®  programs  utilized  In 
calculating  the  aircraft  response  given  a  windshear  and  a  particular 
escape  maneuver.  They  are  listed  in  three  sections.  The  first  section 
contains  the  driver  routine.  The  second  and  third  sections  contain 
the  specific  aircraft  parameter  functions  and  escape  maneuver 
functions  respectively.  The  windshear  function  is  listed  in  Appendix 
A,  The  differential  equation  solving  routine  is  listed  in  Appendix  C. 

A.  DRIVER 

|»»«»»»»»«»«»»*»««»«««»»»«»»»« 

X  Takeoff  Profile  nodel  * 

I***************************** 

X  Thie  progroe  ell  I  calculate  the  inertial  flight  path  angle  for 
X  a  given  oircroft.  It  utilizee  the  function  'eicroburet'  ehich 
X  node  I e  the  find  ehear  encountered  during  the  DFU  Delta  Reel  dent. 

X  ThIe  prograe  ueee  energy  height  theory  to  calculate  Inertial  epace 
poeltlon. 

clear 

cic 

foreat  bank 
foreat  coepact 

X - 

XINITIRLI2RTI0N  RND  ENUIRONnENTRL  CONSTRNTS 
X  The  folloilng  Initialize  conetante  and  varlablee  for  the  nodel. 
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XRiRCRflFT  CONSTRNTS 

t  The  folloeing  constants  are  for  the  particular  aircraft  being 
studied. 

dispC  TRKEOFF  nODEL') 

dispC  ’) 

anseer^O.Oj 

dispCType  aircraft: ') 
dispC  1.  P3  at  90,0001b') 
dispC  2.  P3  at  120,000lb‘) 

dispC  3.  P3  at  135,0001b') 

dispC  4.  LlOlf) 

dispC  5.  T44') 

dispC  6.  P3  at  120,0001b,  Ur- HO') 

ansier- i nput C  Enter  desired  aircraft  by  nuiber:  '); 

If  ansoer— 1, 

[S,K,CDo,Ut,CLo,CLalpha,alphaiax,q,U1,U2,T,thetacliib]-P3LtToff; 
theta-0.087;  Xinitial ization 

el  seif  anseer— 2, 

[S,K,CDo,Ut,CLo,CLalpha,alphaiax,q,U1,U2,T,thetocliib)-P3HuToff; 
theta-0.130;  Xinitialization 

el  seif  ansoer— 3, 

[S,IC,CDo,Ut,CLo,CLalpha,alphaoax,q,Ul,U2,T,thetocliib]-P3HvHvToff; 
theta-0.133;  Xinitialization 

el  seif  ansoer— 4, 

[S,K,CDo,Ut,CLo,CLalpha,alphaoax,q,U1,U2,T,thetacliob]-L1011Toff;XL1 

Oil 

theta-0.2;  Xinitialization 

el  seif  ansoer— 5, 

[S,IC,CDo,Ut,CLo,CLalpha,aiphooax,q,U1,U2,T,thetaciiob]-T44Toff; 

XT44 

theta-0.200;  Xinitialization 

el  seif  ansoer— 6, 

[5,K,CDo,Ut,CLo,CLalpha,alphaoax,q,U1,U2,T,thetacliob]-oodP3HgTo; 
theta-0.130;  Xinitialization 

end 
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IS*plane  fori  surface  area  in  ft^2 
XK*coeff icient  for  CD  calculation 
XCDo*coeff icient  of  drag 
XUt*aircraft  eeight  In  lbs 
n>Ut/32.171;  X  lass  in  slugs 
XCLo*zero  ROD  lift  coefficient 
XCLa I pha* 1 1  ft  curve  s I  ope 
Xalphaiax*stal I  buffet  alpha  •/  approach  flaps 
Xq*noiinal  pitch  rate  of  Sdeg/sec 
XU1 ■rotate  speed  in  ft/s 
XU2"takeoff  safety  airspeed  in  ft/s 
XT*thrust  in  ibf 


X 


X***prograied  initial  gdues*************************************** 
iters*65;  Xnuiber  of  iterations 

del  tat *0.5;  Xtiie  increeent  in  sec 

xstart*-1200;  Xstarting  distance  froi 

ns  center 

hstart«2.0;  Xstarts  at  rotate 


dispC  Initio)  input  Section*) 
dIspC  •) 

dispCThe  fol  losing  values  are  preprogroiied: ') 

dispC  a.  The  flight  path  is  based  on  sox  poser.') 

dispC  b.  The  starting  distance  (dietonce  fros  HD  center)  is:') 

disp(xstart) 

dispC  c.  The  HD  center  is  this  fros  the  end  of  the  runsay:') 

disp(abs(xstart)) 

dispC  d.  The  tise  step  in  sec:') 

disp(deltat) 

dispC  e.  The  nuiber  of  iterotlons  is:*) 
disp( iters) 
dispC  ') 


TU*T/Ut;  Xthrust/seight  ratio 

UL*Ut/5i  Xsing  loading 

dispCThe  fol  losing  values  sere  suppt  ied/caiculated  fros  aircraft 
constants') 


k 
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dispC  a.  Takeoff  Safety  Speed:*) 
disp(U2) 

diepC  b.  Thruet  to  eeight  ie:') 
diepCTU) 

diepC  c.  Uing  loading  is:') 

dispdJL) 

dIspC  •) 

oaneuuer>0.0; 

dispC Select  the  escape  eaneuuep:*) 
dispC  1.  Takeoff  flight  path;  NO  eindshear') 
dispC  2.  Takeoff  flight  path;  UiTH  eindshear') 
dispC  3.  Escape  sith  constant  airspeed') 
dispC  1.  Escape  eith  constant  altitude') 
dispC  5.  Escape  eith  constant  theto') 
dispC  6.  Escape  eith  constant  alpha*) 
oaneuuer«inputC  Enter  selection  nueber:  '); 

if  eaneuuer^^S, 

dispCUhat  ualue  for  escape  theto?*) 
thetaaii«inputC  Enter  in  radians:  *); 
el  seif  eaneuuer-^b, 
diepCUhat  value  for  escape  aipho?’) 
alphaoie-inputC  Enter  in  radions:  '); 


end 

x1>xstart; 

llnitial ization 

h1«hetart; 

linitiai ization 

Ug-U1; 

Iground  speed  initialization  in  ft/s 

gaeai«0.001; 

Xinitial  inertial  flight  path  angle 

rad  (tdeg) 

roe”0. 002377; 

lair  density  in  eiugs/ft''3  •  S.L. 

trigger"0.0; 

futility  check  nueber 

check1*0.0; 

futility  check  nueber 

check2*0.0; 

« 

futility  check  nueber 

m -  " 

IflLGORITHN 
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'  The  actual  algorithe  is  broken  into  seueral  eealier  diuisione  and 
unctions: 

starting  ualues  -  as  it  ieplies 

loop  -  beginning  of  iterating  process 

•indshear  inputs  -  uses  the  function  'eicroburst'  to  giue 

point  X  and  z  axis  eind  speeds 

predictor  -  first  guess  at  the  solution  of  the  gouerning  DEs 
correcter  -  corrected  eolutions  to  the  gouerning  DEs  and 
final  point  ualues. 

output  file  eanagenent  -  as  it  ieplies 

aircraft  control  -  uses  eultiple  control  functions  for 

control  inputs  to  thetOi  thrust,  and  ROD 
final  output  -  os  it  ieplies 


1st anting  ualues: 

pt- 1; 

Xcounter 

Q-0.5*roe*Ur2; 

Xdynaeic  pressure 

CL-Ut/(Q*S); 

Xcoefficient  of  lift 

CD-CDo+K*CL"2j 

Xcoefficient  of  drog 

a  1  pha«(CL-CLo)/Cl.a  1  pha; 

XROR  required  for  1g  flight 

gaeaa«theta-alpha; 

Xalreoss  flight  path  angle 

Es-Ur2/(2*32.171)+hlj 

Xspecific  energy  in  ft 

Esdot«-12.32; 

Xtiee  rate  change  in  specific  energy 

y(l)*x1; 

Xy(1)«x  position  in  ft 

y(2)-hlj 

Xy(2)«oircraft  oltitude  in  ft 

y(5)-ui) 

Xy(5}«aircraft  oirspeed  in  ft/s 

y(6)-gaeaa; 

rad 

Xy(6)"aircraft  oireass  flight  path  in 

RircraftCl, :)*[x1,h1,U1, 

-t2.36,theta,alphaJ,T]; 

1 nert i a  1 ( 1 , : ) -[Ug , 0 , gaeaa , gaea i , 0, 0, 1 , Es , Esdot ] ; 

1  starting  output  files 

XIoop: 

cic 

•hile  theta**999, 

Xioop  control 

pt*pt+l; 

1 - 

Xpresent  pt  in  space 

162 


Xiind  shear  inputs: 

x2-Ug*deitat-<’x1;  Xapproxitate  horizontal  displaceoent 

h2"tan(gaeai)*(x2-x1  )‘*’h1 ;  Xapproxitate  uertical  displaceoent 

if  oaneuMer'"!, 

[Uxl 1, Uhl )]«oicroburst (hi, xt);  XThese  four  calls  to  the  function 
[Ux2t|Uh21]*oicroburst(h2|Xt);  Xoicroburst  finds  the  point  oind 
[Ux12,Uhl2]*oicroburst(hl,x2);  Xshear  for  deteroinotion  of  the 
[Ux22,Uh22]*oicroburst(h2|X2);  Xauerage  block  find  shear, 
else 

Ux11«0;Ux21-0;Ux12-0;Ux22-0;Uh11«0;Uh21-0;Uh12-0;Uh22-0; 

XThe  aboue  is  used  in  conjunction 
Xother  calls  then  no  find  shear  is 
Ides i red. 
end 

UindX(pt)"Ux11;  XUind  shear  oatrices 

UindH(pt)-Uh11; 

dUxdx-((Ux12-Ux1 1 )+(Ux22-Ux21))/(2«(x2-x1)); XThese  forouias 
dUxdh-((Ux22-Ux12)^(Ux21>-Ux11))/(2»(h2>h1));Xcalcutate  the 
dUhdx-((Uhf 2-Uht  t )«(Uh22-Uh2f }}/(2*(x2>xf ));Xdl f ferent iai  change 
dUhdh«((Uh22-Uht2)^(Uh21-Uh11))/(2«(h2-h1));Xin  find  shear. 

y(3)*Ux11;  Xy(3)*find  shear  in  the  x  direction 

ft/s 

y(1)*Uh11;  Xy(1)*find  shear  in  the  z  direction 

ft/s 

I - 

Xpredictor  and  interoediate  values 

ydot "ypr i oes ( y , dUxdx , dUxdh , dUhdx , dUhdh , t het a , T , 5 , n , Q , CL , CD , a  I pha , a  I p 
haoax); 

Xcaii  to  the  DE  function 


for  n»l ;6, 

y1(n)”y(n)‘i’deltat*ydot(n);  XEuler  first  step  fed  predictor 
yhalf(n)”y(n)-*-deltat/2*ydot(n);  XEuler  half  step  fed  predictor 
end 
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alpha*theto-y(6); 
CL'CLo’i’a  I  pha*CLa  I  pha ; 
CD-CDo+K*CL'‘2; 
Q«0.5*roi*yhalf(5)"2; 

% - 


XflOR 

Xcoefficient  of  lift 
Xcoefficient  of  drag 
Xdynatic  pressure 


Xcorrector  and  nee  ualue  foreulation 
ydot "ypr i ees ( yha I f I dUxdx , dUxdh , dUhdx , dUhdh , t he t a I T j S I n I Q , CL I  CO j a i pha 
lalphaeax); 


for  n«l :6, 

y2 ( n ) >yha I f ( n ) +de 1 1 at /2*ydot ( n ) ; XR i chordson  ext rapo lotion 
y(n)«2*y2(n)-y1(n);  Xcorrecter  scheee. 

end 


x1-y(l); 

h1*y(2)j 

Ug-y(5)*cos(y(6))+y(3); 

R0C-y(5)*sin(y(6))+y(1); 

gaiai«atan(ROC/Ug)j 

Es-y(5)^2/64.3«+y(2)i 


Xnee  x  posit  in  ft 
Xnee  oltitude  in  ft  RGL 

Xground  speed  in  ft /sec 
XRate  of  Clieb  in  ft/sec 
Xnee  airoass  FP  angle  In  rad 
Xnee  specific  energy 


Esdot*y(5)*ydot(5)/32.174^ydot(2);  Xnee  trc  of  specific  energy 


aipha*theta-y(6); 
CKLo'^o  I  pha*CLa  1  pha; 
CD-CDo+IC*CL''2; 
Q«0.5*roe*y(5)"2; 

* - 


Xnee  ROR 

Xnee  coefficient  of  lift 
Xnee  coefficient  of  drag 
Xnee  dynaeic  pressure 


Xoutput  file  eanageeent: 

Rircraft  (pt  I  :)-[y(1  ),y(2),y(5), ydotCS),  theta, alpha, ptj]; 

I nert I a  I (pt , : )-[Ug, ROC, y (6) , gaea I , y(3) , y(4) , pt , Es, Esdot ] ; 

hoee 

dispC  pt  X  h  U 

olrspeed  occel  ') 

op1*'[pt,Rircraft(pt,1),Rlrcraft(pt,2),Rircraft(pt,3),Rircraft(pt,4)] 

I 
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dispCopI) 

dispC  pt  theta  alpha  gaeaa 

gaeai ') 

op2*[pt , fl I rcraf t (pt , 5) , fl I rcraf t (pt , 6) , I nert i a  I (pt  ^  3) J  nert I  a  I (pt , 

1)]; 

di3p(op2} 

diepC  pt  Ug  ROC  Thrust'} 

op3*[pt , I nert i a  I (pt , 1) J  nert i a  I (pt « 2) , R i rcraf t (pt , 8) ] ; 
disp(op3) 

diepC  pt  Ee  Eedot') 

op5"[pt 1 1 nert i a  I (pt 1 8) J nert i a  I (pt , 9) ] ; 
disp(op5) 

diepC  pt  Ux  Uh') 

op4"[pt|  inertial(pt«5)  Jnertial(pt«6)]; 
di8p(op4) 

% - 


laircraft  control: 

Ifunctione  called  for  theta,  thrust 
if  y(2)>50,  land  alpha  control, 

if  ROC  <-0, 
trlgger-1; 
eleeif  y(5)<U2, 
trlgger-1; 
end 

eleeif  y(5)<U1, 
trigger-1; 
end 

If  trigger— 0, 
thet  o-t  heta‘^q*0 . 8*de  I  tot ; 
if  thetaHhetacI ieb, 
theta-thetaci ieb; 
end 
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else 

if  eaneuuer^'S, 

[theta«check1]"ToffCRS(q,deltat,theta,alpha,alphaeax,U2,yjCheck1 ) 
el  self  eaneuver^M, 

[t hetOi checki ]-Tof  fCRLT (q, de I  tat , theta, a  I pha, a  I phaeax, y, check 1 , ydot ) 
el  seif  eaneuaer^S, 

[theta]*ToffCTH(qideltat,theta,aipha,aiphaeax,thetaaie) 
el  seif  eaneaueM^d, 

[theta]"Tof f CR0fl(q, de I  tat , theta, a  I pha, a  I phaeax, a  I phaa I ■) 
end 
end 

If  pt*"lter8, 
theta-999; 
eleeif  pt>5. 

If  y(2)<0.0, 
theta-999; 
end 
end 

end 

g - 

loutput  files: 

anseer- input ('Do  you  eant  a  flight  path  plot?  1-VES  '); 
if  anseer— 1, 

plot(Rircraft(:,1),Rircraft(:,2)) 

pause 

ansser-O.O; 

end 

aneeer- input ('Do  you  oant  x,  h,  U,  ond  airspeed  accel?  1-VES  '); 

If  anseer— J, 

diepC  pt  X  h  U 

airspeed  accel  ') 

op1-[Rlrcraft(:,7),Rlrcraft(:, l),fllrcraft(:,2),Rlrcraft(;,3),Rlrcraf 
t(:,4)]; 


Inuiber  of  iterations 


lehiie  olgorithe  loop 
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disp(opl) 

else 

end 

aneeer" input ('Do  you  sont  thetOi  alpha,  gaoaa,  and  gaeai?  t«VES  '); 
if  aneeer”!, 

diepC  pt  theta  alpha  gaaaa 

gaeai ') 

op2*[Rircraft(:,7),flircraft(:,5),flircraft(:,6), Inert iai ( ; ,3), Inertia 

l(:,1)]; 

diep(op2) 

else 

end 

anseer* I nput ( ' Do  you  eant  Ug,  ROC,  and  Thrust?  1*VE5  '); 
if  anseer««1, 

dispC  pt  Ug  ROC  Thrust') 

op3*[Rircraft (:,?), Inert lai(:,1), Inert lal(:,2),Rircraft(;,6)]j 

disp(op3) 

else 

end 

anseer" input ('Do  you  eant  Es  and  Esdot?  1*VES  '); 
if  ansier**!, 

dispC  pt  Es  Esdot’) 

op5*[R I rcraft (:,?), Inertiai(:,8), Inert ial(;,9)]; 

disp(op5) 

else 

end 

onseer* input ('Do  you  eant  Ux  and  Uh?  I^VES  '); 
if  anseer**!, 

dispC  pt  Ux  Uh') 

op4«[Rircraft(:,7), lnertial(:,5),inertial(:,6)]; 

disp(opl) 

else 

end 
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B.  AIRCRAFT  PERFORMANCE  PARAMETER  FUNCTIONS 


function[S|iC>C0o,Ut,CLo,CLalpha,alphaiaXjq,U1  jU2Jjthetacl  iib] 


-P3UToff() 

XP3  at  QOjOOOIbSj  app  fiapSj 
S-1300; 

K-0. 05041; 

CDo-0.0551; 

Ut -90000; 

CLo-t.OOO; 

CLaipha>5.73; 

alphaiax-0.244; 

q-0.0873; 

Ul-204; 

U2-220; 

T-33400; 

thetoci iab-0. 175; 


gear  up 

Xpiane  fort  eurface  area  in  ft''2 
Xcoefficient  for  CD  calculation 
Xcoefficient  of  drag 
Xaircraft  eelght  in  Ibe 
Xzero  flOR  lift  coefficient 
XI I  ft  curve  elope 

Xetall  buffet  olpha  ■/  approach  flape 
Xnotinal  pitch  rate  of  5deg/8ec 
XLiftoff  epeed  (I21kt8) 

Xtakeoff  eafety  airepeed  in  ft/e 
Xthruet  in  Ibf 

Xeater  line  deck  angle  in  rad 


function[S/KjC0o,Ut 
-P3HvToff() 
XP3  at  I20,000lbe; 
S-1300; 

K-0. 05041; 
CDo-0.0551 ; 

Ut-I 20000; 

CLo-1 .000; 
CLalpha-5.73; 
alphaoax-0.244; 
q-0.0873; 

UI-214; 

U2-227; 

T-33400; 

thetaci  iib-O.nS; 


, CLoj  CLa I pha, 0 1 phaiax, q, Ul ^  U2 thetac I i ob] 

app  flaps,  gear  up 

Xpiane  fore  surface  orsa  In  ft^2 
Xcoefficient  for  CD  calculation 
Xcoefficient  of  drag 
Xaircraft  eeight  in  Ibe 
Xzero  flOR  lift  coefficient 
Xlift  curve  elope 

Xetoil  buffet  alpha  •/  approach  flaps 
Xnotinal  pitch  rate  of  Sdeg/eec 
XLiftoff  epeed  (I27kte) 

Xtokeoff  eafety  airepeed  in  ft/e 
Xthruet  in  Ibf 

Xiater  line  deck  angle  in  rad 
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f unct I  on [ S , K , CDo , Ut I CLo 
-P3HuHvToff() 
XP3  at  135, 000 lbs,  opp 
S-1300; 

K-O. 05011; 

CDo*0.0551 ; 

Ut-135000; 

CLo-1.000; 

CLalpha>5.73; 

alphaiax«0.211; 

q*0.0673; 

01-229;. 

02-239; 

T-33100; 

thetocl itb-0. 1715; 


I  CLo I pho , a  I phonox , q , 0 1 , 02 , T , t  hot  oc I i tb ] 
flaps,  gear  up 

Xplane  fore  surface  area  in  ft ''2 
Xcoefficient  for  CD  calculation 
Xcoefficient  of  drag 
Xaircraft  eelght  in  lbs 
Xzero  ROR  lift  coefficient 
XI ift  curue  slope 

Xstall  buffet  alpho  «/  approach  flaps 
Xnoeinai  pitch  rate  of  5deg/sec 
XLiftoff  sqeed  (J36kts). 

Xtakeoff  safety  airspeed  in  ft/s 
Xthrust  in  Ibf 

Xeater  line  deck  angle  in  rad 


f unct i on[ S , K , CDo , Ut , CLo 
-■odP3HuTo() 

XP3  at  120, 000 lbs,  app 
S-1300; 

K-O. 05011; 

CDo-0.0551; 

Ut-1 20000; 

CLo-1.000; 

CLalpho-5.73; 

alphaiax-0.211; 

q-0.0673; 

01-212; 

02-251; 

T-33100; 

thetaci l■b-0.175; 


,  CLa  I  pha,  a  iphaMX,q,  01,02,7,  thetocl  ieb] 

flaps,  gear  up.  Increased  rotate  speed. 
Xplane  fore  surface  area  in  ft^2 
Xcoefficient  for  CO  calculation 
Xcoefficient  of  drag 
Xaircraft  oeight  in  lbs 
Xzero  ROR  lift  coefficient 
XI ift  curue  elope 

Xetali  buffet  olpho  e/  approach  flaps 
Xnoeinai  pitch  rote  of  5deg/eec 
XLiftoff  speed  (115kts) 

Xtakeoff  safety  airspeed  in  ft/s 
Xthrust  in  Ibf 

Xeater  line  deck  angle  in  rad 
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fund ion[S,IC«CDo«Ut>CLo,CLalpha,alphaiaXiqiU1,U2J,thetc:cl  iab] 
•LIOIIToffO 

IL1011  at  462|000lb8|  Rpp  flaps,  gear  up 


S-4578J 

IC-0.059; 

CDo*0.098; 

Ut -162000; 
CLo-0.532; 
CLalpha-1.96; 
alphaiax-0.3Hi 

q-0.087j; 

U 1-238; 

U2-255; 

T-126000; 
thetaci i■b-0.209; 


Ip  I  one  fore  surface  area  in  ft ^2 
Icoefficient  for  CD  calculation 
Icoefficient  of  drag 
laircraft  eeight  in  lbs 
Izero  flOfl  lift  coefficient 
I.  ift  curue  slope 

Istall  earning  aipha  •/  approach  flaps 

l(18deg) 

Inoeinal  pitch  rate  of  Sdeg/sec 
Xrotate  speed  In  ft/sec  (lllkts) 
Xtakeoff  safety  airspeed  ft/s  (tSIkts) 
X  thrust  in  Ibf 

Xinltial  deck  angle  rad(12deg) 


funct lon[S|K,CDO|Ut,CLO|CLalpha,alphaeax,q,U1,U2|T,thetacl leb] 
-TIIToffO 

XT11  at  761 7 lbs,  app  flaps,  geor  doen 


S-2I0; 

K-O.OlO; 

CDo-0.010; 

«t-7817; 

CLo-0.0521; 

CLalpha-6.21; 

alphaeax-0.227; 

q-0.0873; 

U1-t52; 

U2-180; 

-3023; 

thetaci l■b-0.262; 


Xplane  fore  surface  area  In  ft "^2 
Icoefficient  for  CD  calculation 
Icoefficient  of  drag 
Xaircroft  eeight  in  lbs 
Izero  ROR  lift  coefficient 
II ift  curve  slope 
Xstoll  buffet  alpha  e/flaps  up 
Inoeinal  pitch  rote  of  Sdeg/sec 
Xliftoff  airspeed  in  ft/s(90kts) 
Xtakeoff  sofety  speed  ft/s  (107kts) 
Xinltial  thrust  in  Ibf 
Xinitia!  deck  angle  in  rad  (15  deg) 
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C.  FLIGHT  PATH  CONTROL  AND  ESCAPE  MANEUVER  FUNCTIONS 


fund  lonCthetaoutjCheckI  ] 

■To f  f CRS ( q ^  de 1 1  at , t  het  a ^  a i pha , a  I phoiax >  U2 , y  ^  check 1 } 
XRIncraft  escape  eaneuuer  const! Ing  of  constant  airspeed, 
thetaout “theta; 

If  check 1““0, 

If  y(5)<U2, 

t  het  aout “t  het  a-q*de 1 1  at ; 
if  thetaout<0.0^ 
the! aout *0.0; 
end 


else 

thetaout*q*deltat ; 
check1"1 ; 
end 
else 

If  y(5)>-(U2+8.1), 
thetaout  ■!  het  a<^q*dettat ; 
elself  y(5)<-(U2-e.1) 
thetaout ■theta-q*de I  tat ; 
end 
end 

If  alpha>alphaeax, 
thetaout*theta-(aipha-alphaeax); 
end 
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fund  i  on[  t  het  aout ,  check  t  ] 

■To  f  f CRLT ( q , de 1 1  at , t  het  0 , 0 1 pho « a i phoeax , y , check 1 , ydot ) 
IRircroft  control  using  constant  altitude  escape  eaneuaer. 
t het aout "theta; 
if  check 1«“0; 

check l"y(2)+20; 
end 

if  y(2)<checkl-10, 
if  ydot(2)<»0, 
t  het  aout "t  het  a+q*de 1 1  at ; 
end 

elseif  y(2)>check1‘^10j 
If  ydot(2)>*0, 
t het aout "t het a-q*de I  tat ; 
end 
end 

if  thetaout<Oj 
thetaout*0; 
end 

if  alpha>aiphaiax, 
t  het  aout "t  het  a- ( a i pha-a i phaeax ) ; 
end 


f unct i on [ t het aout ] "To f f CTH ( q , de 1 1 ot , t het a , a  I pho , a  I phaeax , t het 00 i ■ ) 
XConstant  theta  escape  eaneuuer. 
if  theta<thetaaie, 
t  het  aout  "t  het  a‘'’q*de  I  tat ; 
else 

thetaout"thetaaio; 

end 

if  aipha>alphaeax, 
t  het  aout "t  het  a- ( a  I pha-a I phaeax ) ; 
end 
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f unct i on [ t het aout ] -To f f CflOfl ( q I de 1 1 at , t het a , a  I pha , a  I phaiax , Q I phoQ i ■ ) 

IConstant  alpha  escape  eaneuver. 
if  alpha<alphaaii| 
that  aout -t  hot  a* (a  I phaa i t-a I pha } ; 

If  thetaout>(theta+q*dcllat), 
t  hot  aout -t  hot  a+q*de 1 1  at ; 
end 

else 

thetaout-theta-(a I pha-a I phaa !■) ; 

If  thetaout<(theta-q*deltat), 
thetaout-theta-q*de I  tat  j 
end 
end 

If  alpha>alphaiaxi 
t  het  aout -t  he t  a- ( a  I pha-a I phaiax ) ; 
end 
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APPENDIX  E 


ON-STATION  ALGORITHMS 

Contained  in  this  appendix  are  the  MatLab®  programs  utilized  in 
calculating  the  response  of  a  P-3  encountering  a  windshear  while 
operating  close  to  the  surface.  They  are  listed  in  three  sections.  The 
first  section  contains  the  driver  routine.  The  second  and  third 
sections  contain  the  aircraft  parameter  function  and  flight  path 
control  function  respectively.  The  windshear  function  is  listed  in 
Appendix  A.  The  differential  equation  solving  routine  is  listed  in 
Appendix  C. 

A.  DRIVER 

|«»«««»«««««««««»»««»»«»»««»»« 

%  P3  Loiter  Profile  Hodel  « 
leeeeeeeeeeeeeeeeeeeeeeeeeeeee 

clear 

cic 

foreat  bank 
foreat  coepact 

% - 

fINITiflLIZflTION  RHD  ENUiflONnENTRL  CONSTRNTS 

%  The  folloeing  initialize  conetante  and  ooriablee  for  the  eodel. 
IRiRCRRFT  CONSTRNTS 

I  The  folloeing  conetante  are  for  the  particular  aircraft  being 
studied. 

diepC  P3  ON-STRTION') 
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dlspC  •) 
check 1 ”0.0; 

U-0.0; 

answer  "t.O; 

check1"input( '  Enter  on-station  loiter  altitude(ft):  ')j 
U*input('  Enter  loiter  airspeed  in  knots:  *); 

U"U*1.687;  Iconoert  knots  to  ft/sec 

onseer* input ('Enter  0  for  no  eindsheor:  '); 

[ S , K , CDo , Ut , CLo I CLo I pho , 0 1 phoeox I q , T , t het 0 ] ■P3onst 0 ( U ) ; 

Xinitiof izQtion 

X5*plane  fore  surface  area  in  ft^2 
Xk”coeff icient  for  CD  calculotion 
XC0o*coeff icient  of  drag 
XUt ‘aircraft  eeight  in  lbs 
n"Ut/32.174]  X  eass  in  slugs 
XCLo‘zero  ROR  lift  coefficient 
XCLalpha‘1 ift  curue  slope 
Xalphaeax‘3tal I  buffet  alpha  e/  approach  flaps 
Xq‘noeinal  pitch  rate  of  Sdeg/sec 
XUI-rotote  speed  in  ft/s 
XU2*takeoff  safety  airspeed  in  ft/s 
XT‘thrust  in  Ibf 


X***prograeed  initial  go  lues*************************************** 
iter8‘H5;  Xnueber  of  iterations 

deltat‘0.25;  Xtiee  increeent  in  sec 

X8tart‘-5000;  Xstarting  distance  froe 

ns  center 

dispC  Initial  Input  Section') 
dispC  ') 

dispCThe  fo I  losing  ualues  are  preprogroeeed: ') 

dispC  a.  The  flight  path  is  based  on  constant  poser.') 

dispC  b.  The  storting  distance  (distance  fron  116  center)  is:') 

di8p(x8tart) 

dispC  c.  The  tine  step  in  sec:') 


175 


dispCdeltat) 

dispC  d.  The  nueber  of  iterations  is:') 
dispC iters) 
dispC  ‘) 


TU"T/Utj  Ithrust/eeight  rotio 

UL«Ut/Si  Xoing  ioading 

dispCThe  follooing  uoiues  eere  suppi  ied/ca!cu toted  froe  aircraft 
constants') 

dispC  a.  Loiter  Rir  Speed:') 
disp(U) 

dispC  b.  Thrust  to  eeight  is:') 
disp(TU) 

dispC  c.  Uing  ioading  is:') 

disp(UL) 

dispC  ') 


xl"xstart; 
h 1 "check 1; 
Ul-U; 

Ug-Uj 

goto i "0.001; 
rad  (I deg) 
roi"0. 002377; 
trigger"0.0; 
check2"0.0; 


SinitiaMzation 
linitiol ization 
linitial ization 

Iground  speed  initialization  in  ft/s 
linitial  inertial  flight  path  angle  in 

fair  density  In  siugs/fC'3  •  S.L. 
futility  check  nueber 
futility  check  nueber 


RLGORITHn 

The  actual  algorithe  is  broken  into  seuerol  sealler  divisions  and 
unctions: 

storting  values  -  as  It  leplles 

loop  -  beginning  of  iterating  process 

eindshear  inputs  -  uses  the  function  'eicroburst'  to  give 

point  X  and  z  axis  oind  speeds 

predictor  •  first  guess  at  the  solution  of  the  governing  DEs 
corrector  -  corrected  solutions  to  the  governing  DEs  and 
final  point  values. 

output  file  eanageeent  -  as  it  iepiies 

aircraft  control  -  uses  eultiple  control  functions  for 

control  inputs  to  theta,  thrust,  and  ROR 
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final  output  -  as  it  ioplieo 


Xstarting  ualues: 

pt«l; 

Xcounter 

Q-0.5*roi*Ur2; 

Xdynaoic  pressure 

CL-Ut/(Q*S)j 

Xcoefficient  of  lift 

CD-CDo+K*CL"2; 

Xcoefficient  of  drag 

a  1 pha*(CL-CLo)/CLa i pha; 

XROR  required  for  1g  flight 

gaeaa*theta-alpha; 

Xai reass  flight  path  angle 

E8-Ur2/(2*32.171)+hlj 

Xspecific  energy  in  ft 

t8(lot»U; 

Xt'iee  rate  change  in  spec) Tie  energy 

y(1)«xl; 

Xy(1)*x  position  in  ft 

y(2)«hl; 

Xy(2)*aircraft  altitude  in  ft 

y(5)-Ul; 

Xy(5)*aircraft  airspeed  in  ft/s 

y(6)-gaeaaj 

rad 

Xy(6)"aircraft  airnass  flight  path  in 

Rircraftdi  :)"[x1 1  hi  theta, alpha, t,T]; 

i nert i a l( 1 , : )«[Ug, 0, gaeaa, gaea i ,0,0, 1 ,E8, Esdot ] ; 

X  starting  output  files 

1 - 

XIoop: 

cic 

ehile  theta"*999. 

pt-pt+1; 

Xpresent  pt  in  space 

Xoind  shear  inputs: 

x2*Ug*deltot+xl ; 

Xapproxleate  horizontal  displaceeent 

h2"tan(gaeai  )*(x2-x1  )‘^h1 ; 

Xapproxieate  uertical  displaceeent 

if  anster'^O, 

[Ux11,Uh11]*picrobur3t(h1,x1);  IThese  four  colls  to  the  function 
[Ux21,Uh21]«iicrobuP8t(h2|Xl);  fsicroburst  finds  the  point  oind 
[UxI2|UhI2]*eicroburst(hI|X2);  Xsheor  for  detereination  of  the 
[Ux22,Uh22]*iicrobur8t(h2,x2);  Xauerage  block  oind  shear, 
if  pt<*8,  XThese  next  if  logic  coeeands 

Ux11”Ux11*pt/8;Uh11"Uh11*pt/6;  Xcontrol  the  entry  and  exit 
Ux12”Ux12*pt/8;Uh12"Uh12*pt/8;  Xeindshear  effects. 


el  seif  y(1)>*abs(xstart), 

If  check2“0, 
check2*pt ; 
end 

If  (6-pt+check2)>"0, 

Ux21-Ux21*(6-(pt-check2))/8;Uh21-Uh21*(8-(pl-check2))/8; 

Ux22-Ux22*(8-(pt-check2))/8jUh22-Uh22*(8-(pt-check2))/8; 

else 

Uxt1-0;Ux21-0;Ux12-0;Ux22-0;Uh1t-0;Uh21-0;Uh12-0jUh22-0; 

end 

end 

else 

Ux11-0jUx21-0;Ux12-0jUx22-0;Uht1-0;Uh21-0;Uh12-0;Uh22-0; 

IThe  aboue  is  used  in  conjunction 
lot her  colie  ehen  no  oind  shear  is 
Ides i red. 

end  lonseer 

UindX(pt)-Ux1]j  lUind  shear  latrices 

UindH(pt)-Uhnj 

dUxdx-((Ux12-Ux1 1 )>^(Ux22-Ux21 ))/(2*(x2-xl ));ITheee  forouias 
dUxdh-((Ux22-Mx12)^(Ux21-Ux11))/(2«(h2-h1));lcolculate  the 
dUhdx-((Uhl2-Uht I )^(Uh22-Uh2l ))/(2*(x2-xt ));ldi fferent iol  change 
dUhdh-((Uh22-Uh12)«(Uh21-Uh11))/(2«(h2-h1));lin  oind  shear. 

y(3)-Ux11;  ly(3)*sind  sheor  in  the  x  direction 

ft/s 

y(4)>Uh11;  ly(4)”sind  shear  in  the  z  direction 

ft/s 

I - 

Ipredictor  and  interiediate  values 
y dot  "y pr  i  MS  ( y ,  dUxdx  I  dUxdh ,  dUhdx ,  dllhdh ,  t  hst  a ,  T ,  S ,  n ,  Q I  CL ,  CD , 
aipha^alphasax);  leal  I  to  the  DE  function 

for  n*l :6, 

y1(n)*y(n)'*'deltat*ydot(n);  lEuler  first  step  fed  predictor 
yhalf(n)*y(n)<*'deltat/2*ydot(n);  lEuler  half  step  fed  predictor 
end 
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alphQ*th8ta-y(6); 
CL^CLo-t-a  I  pha*CLa  I  pha; 
CD-CDo+»C*CL"2; 
0-0.5*rot*yhalf(5)"2; 

X - 


XROR 

Xcoefficient  of  lift 
Xcoefficient  of  drag 
Xdynaoic  pressure 


Xcorrector  and  nee  ualue  foreulation 
ydot  "ypr  i  ees  ( yha  I  f ,  dUxdx ,  dUxdh ,  dUhdx ,  dUhdh  1 1  het  a ,  T ,  S ,  11 ,  Q  i 
CL,CD,alphaialphaeax); 

for  n-1 ;6, 

y2(n)*yha I f (n)+de I tat/2*ydot (n) ;XRi chardson  extrapo lotion 
y(n)»2*y2(n)-yl(n);  Xcorrector  scheee. 

end 


xi-y(l)j 

hl-y(2); 

Ug-y(5)*co8(y(6))+y(3)j 

R0C-y(5)*sln(y(6))+y(1); 

gaaai-atan(ROC/Ug); 

Es-y(5)"2/61.348+y(2); 


Xnee  x  posit  in  ft 
Xneo  altitude  in  ft  R6L 
Xground  speed  in  ft /sec 
XRate  of  Clieb  in  ft/sec 
Xneo  aireass  FP  angle  in  rad 
Xnee  specific  energy 


Esdot*y(5)*ydot(5)/32.174+ydot(2);  Xnee  trc  of  specific  energy 


alpha«theto-y(6); 
CL^CLo'^a  i  pha*CLo  I  pha ; 
CD-C0o+K*CL"2; 
0-0.5*roe»y(5)"2; 

X - 


Xnee  ROR 

Xnee  coefficient  of  lift 
Xnee  coefficient  of  drag 
Xnee  dynoeic  pressure 


Xoutput  file  eanageeent: 

Rircraftfpt,  :)-[y(1  ),y(2),y(5), ydot  (5), theta, alpha, ptj]; 
Inert iol (pt, :)-[Ug,R0C,y(6),goeai,y(3),y(4),pt,E8,E8dot]; 
hoee 

dispC  pt  X  h 

airspeed  accei  ') 

op1"[pt,Rircraft(pt,1),Rircraft(pt,2),Rircraft(pt,3), 

Rircraft(pt,1)]j 

dispCopI) 

dIspC  pt  theta  alpha 

ganai ') 

op2*[pt , R i rcraf t (pt , 5) , R i rcraf t (pt , 6) , I nert i a i (pt , 3) , 


gaeaa 


U 
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InertiaKpt,  4)]; 
disp(op2) 

dispC  pt  Ug  ROC  Thrust') 

op3« [pt 1 1 nert i a  I (pt  J ) J  nert i a  I (pt « 2) « R i rcro f t ( pt « 8 } ] ; 
di8p(op3) 

dispC  pt  Es  Esdot') 

op5*[pt,  lnertial(pt,8)  Jnertial(pt,9)]; 
di3p(op5) 

dispC  pt  Ux  Uh') 

op4”[pt  Jnertial(pt«5)j  Inert  iaf  (pt,6)]; 
dl8p(op1) 

I - 

laircraft  control: 

[ t het 0 J ] -CRLT ( q , de 1 1 at , t het 0 , 0 1 pho 1 0 1 phosox , y , check 1 , 
ydotJ,U) 

If  pt**lter8,  Inuiber  of  iterations 

theto*999; 
el  seif  y(1)>5D00« 
theta*999; 
end 

end  lehile  algorithi  loop 

I - 

Soutput  files: 

ansoer- input  ('Do  you  oant  a  flight  path  plot?  1*>VES  '); 
if  ansier”1, 

plot(Rircraft(:,niRircraft(:,2)) 

pause 

ansoer*0.0; 

end 

anster" input ('Do  you  lant  x,  h,  U,  and  airspeed  accel?  1*VES  '); 

If  anseer—l, 

dispC  pt  X  h  U 

airspeed  accel  ') 

opl“[fllrcraft(:,7),Rlrcraft(;,l),Rlrcraft(:,2),Rlrcraft(:,3), 

Aircraft (;,1)]; 
di8p(op1 ) 
else 
end 
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ansier" input ('Do  you  oont  theta,  alpha,  gaeaa,  and  gaeai?  1*VES  '); 
if  anseer”!, 

diepC  pt  theta  alpha  gaeaa 

ganai ') 

op2"[flircraft(:,7),fllrcpaft(:,5),fllrcraft(:,6), Inert ial(; ,3), 

Inert ial( : ,1)]; 
dlsp(op2) 
else 
end 

aneeer" I nput (' Do  you  eont  Ug,  ROC,  and  Thrust?  1*VES  '); 

If  aneeer”!, 

dispC  pt  Ug  ROC  Thrust') 

op3*[Rircraft(:,7), inert ial(:, t), lnert{al(:,2),Rircraft(:,6}]; 

dl8p(op3) 

else 

end 

aneeer* I nput ( 'Do  you  eant  Ee  and  Esdot?  1*VES  '); 

If  aneeer”!, 

diepC  pt  Ee  Eedot') 

op5*[Rlrcraft(:,7),  inertial(.‘,8).  Inert  ial(:,9)]; 

diepCopS) 

else 

end 

aneeer* I nput ('Do  you  eant  Ux  and  Uh?  1*VE5  '); 

If  aneeer**!, 

diepC  pt  Ux  Uh') 

op1*[Rircraft(:,7), lnertial(:,5), lnertial(:,6)]; 

di8p(op1) 

else 

end 
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B.  AIRCRAFT  PERFORMANCE  PARAMETER  FUNCTION 


funct ion[S,ICjCDo,Ut,CLo,CLalpha,alphaMax,q,T^theta]"P3on8ta(U) 
IP3  at  120,000lb8, flaps  and  gear  up. 

S*1300;  Xpiane  fori  surface  area  in  ft^2 

i("0. 05041;  Icoefficient  for  CO  calculation 

C0o*0.0213;  ^coefficient  of  drag 

Ut*120000;  laircraft  leight  in  lbs 

CLo*0.350;  Xzero  AOfl  lift  coefficient 

CLalpha"6.25;  Xiift  curue  slope 

aiphaiax-0.209j  Xstail  buffet  alpha  ■/  flaps  up 

q"0.0873;  Xnoiinal  pitch  rate  of  Sdeg/sec 

T- .  5*0 . 002377*U''2*S*(C0o+K*(Ut /( .  5* .  002377*U"2*S)  )"2)  ; 

Xthrust  in  Ibf 


theta-(Ut/(.5*.002377*U"2*S)-CLo)/CLalpha; 


Xiater  line  deck 
angle  in  rad 
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C.  FLIGHT  PATH  CONTROL  FUNCTION 


functlon[thetaoulJout]-CRLT(q,deltat|thetaialphaialphaiax,y, 

check l^ydot^T^U) 

XRircraft  control  using  constant  altitude  escape  oaneuger. 
thetaout "theta; 

Tout«T; 

if  y(2)<check1-5, 

If  ydot(2)<-0, 

thetaout*theta+q/2 . 5*de I  tat ; 
end 

el  self  y(2)>checkH5« 

If  ydot(2)>"0, 

thetaout ■theta-q/2 . 5*de I  tat ; 
end 
end 

I f  abe(ydot(2)*deltat)>abe(y(2)>checkt)« 
thetaout'thet  aout -e I gn ( ydot ( 2 ) ) *q/ 1 *de 1 1  at ; 
end 

if  aipha>alphaiaX| 
t het aout -t het a- (a  I pha>a i phaeax) ; 
end 

If  y(5)<-(U-67.5), 

Tout-T+deltat/2*33400; 

end 

If  Tout >33400, 

Tout "33100; 
end 
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APPENDIX  F 


APPROACH  TO  LANDING  ENCOUNTER  GRAPHIC  DATA 

This  appendix  contains  the  calculated  performance  of  different 
aircraft  and  weight  combinations  upon  encountering  a  microburst 
windshear  during  an  approach  to  landing.  Each  figure  is  for  a 
particular  aircraft  performing  a  specified  escape  maneuver.  Each 
figure  contains  four  graphs  depicting  altitude,  theta,  alpha,  airspeed, 
specific  energy,  and  time.  The  abscissa  in  all  graphs  is  x  (distance 
from  microburst  center). 
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Figxire  F2.  P-3  at  89,5001bs  gross  weight  performing  a 
constant  altitude  escape 
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Figure  F3.  P-3  at  89,5001bs  gross  weight  performing  a 
consteuit  5"  theta  escape. 
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Figlire  F4.  P-3  at  89,5001bs  gross  weight  performing  a 
constant  80  theta  escape. 
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Figiire  F5.  P-3  at  89,5001bs  gross  weight  performing  a 
constant  10®  theta  escape. 
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Figure  F6.  P-3  at  89,5001bs  gross  weight  performing  a 
constant  15"  theta  escape. 


0 


o  cn  d 

tlincCs) 


>(«) 


-500  0  500  1000  1500  2000  2500  3000  3500  4000  4500  5000  5500 


x(«) 


Figure  F7.  P-3  at  89,5001bs  gross  weight  performing  a 
constant  12  unit  AOA  escape. 
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Figure  F8.  P-3  at  89,500lbs  gross  weight  performing  a 
constant  15  unit  AOA  escape. 
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Figure  FIO.  P-3  at  114,0001bs  gross  weight  performing  a 
constant  airspeed  escape. 


194 


timvCs) 


12 


-2  I  »  i  »  i  ■  i — ■ — |— ' — I — '-—I — ' — i — I — i — ' — i — >--1 — • — i — r 


-500  0  500  10)0  1500  2000  2500  3000  3500  4000  4500  5000  5500 


x(«) 


Figure  F12.  P-3  at  114,0001bs  gross  weight  performing  a 
constant  5«  theta  escape. 
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Figiire  F13.  P-3  at  114,0001bs  gross  weight  performing  a 
constant  8o  theta  escape. 
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Figure  F15.  P-3  at  114,0001bs  gross  weight  performing  a 
constsuit  15a  theta  escape. 
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Figure  F18.  P-3  at  114,0001bs  gross  weight  performing  a 
constant  20  unit  AOA  escape. 
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Figiire  F23.  Three-engine  heavy  eurline  transport 
perfomning  a  constant  15«>  theta  escape. 
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Figxjire  F24.  T-44  performing  a  constant  mrspeed  escape. 
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Figure  F25.  T-44  performing  a  constant  altitude  escape. 
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Figxire  F29.  T-44  performing  a  constant  10  unit  AOA 
escape. 
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Figure  F30.  T-44  performing  a  const€uit  16  unit  AOA 
escape. 
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Figure  F31.  T-44  performing  a  constent  25  unit  AOA 
escape. 
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Figxire  F32.  P-3  at  89^5001bs  gross  weight  with  three 
engines  operating  performing  a  constant  5o  theta  escape. 
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Figure  F33.  P-3  at  89,5001bs  gross  weight  with  three 
engines  operating  performing  a  constant  10«  theta  escape 
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Figui  e  F34.  P-3  at  89,5001bs  gross  weight  with  three 
engines  operating  performing  a  constant  15®  theta  escape. 
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APPENDIX  G 


TAKEOFF  ENCOUNTER  GRAPHIC  DATA 

This  appendix  contains  the  calculated  performance  of  different 
aircraft  and  weight  combinations  upon  encountering  a  microburst 
windshear  at  or  immediately  after  lift  off.  Each  figure  is  for  a 
particular  aircraft  performing  a  specified  escape  maneuver.  Each 
figure  contains  four  graphs  depicting  altitude,  theta,  alpha,  airspeed, 
specific  energy,  and  time.  The  abscissa  in  all  graphs  is  x  (distance 
from  microburst  center). 
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Figure  61.  P-3  at  90,0001bs  gross  weight  performing  a 
constant  airspeed  maneuver  during  takeoff. 
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Figure  G3.  P-3  at  90,0001bs  gross  weight  performing  a 
constant  12  unit  AOA  maneuver  during  takeoff. 
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APPENDIX  H 


P-3  FLIGHT  SIMULATOR  DATA 


\ 
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FUEL 

TANK 

FUEL 

LOAD 

ENGINE 

1 

2 

3 

M 

BRAKES 

1 

3937 

SHP 

^B3 

MBM 

^99 

597 

PILOT 

7 

3937 

TIT 

5M9 

5M9 

551 

564 

COPILOT 

w 

FUEL  FLOU 

8M9 

BH9 

855 

893 

FLIGHT 

3 

3937 

IAS  115 

U 

3935 

ALT  490 

5 

996 

RUNUAY-  100 

MOFFETT 

NAS  RUY 

32R 

HDG  265.9 

FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00 : 00 : 00 

MET  TIMER 

00:03:58 

CONF I GURAT I ON/COND I T 1 ONS 

GROSS  UEIGHT 

69960 

PRESSURE  ALTITUDE 

524.0 

C.G. 

20.00 

CALIBRATED  AIRSPD 

1 14.  7 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

114.72 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 

195.41 

MACH  NUMBER 

0.  18 

FLIGHT/ AERO 

PITCH  ANGLE 

7.3 

BANK  ANGLE 

-1.4 

ANGLE  OF  ATTACK 

6.4 

SIDESLIP 

-1.4 

HEADING  ANGLE 

283.2 

RATE  OF  CLIMB  (FPM) 

181 

PITCH  UELOCITY  (D/S) 

-0. 156 

PITCH  ACCELERATION 

0.0004 

ROLL  UELOCITY  (D/S) 

-0. 164 

ROLL  ACCELERATION 

-0.0024 

YAU  UELOCITY  (D/S) 

-0.273 

YAU  ACCELERATION 

-0.0023 

NORTH-SOUTH  UELOCITY 

111.71 

NORTH-SOUTH  ACCEL 

-0.373 

EAST-UEST  UELOCITY 

-160.32 

EAST-UEST  ACCELERATION 

-0. 138 

UERTICAL  UELOCITY 

-2.99 

UERTICAL  ACCELERATION 

0.579 

LONGITUDINAL  ACCEL 

0.  1216 

TOTAL  PITCHING  MOMENT 

336 

LATERAL  ACCEL 

0.0117 

TOTAL  ROLLING  MOMENT 

-1664 

UERTICAL  ACCEL  (G'S) 

-0.9739 

TOTAL  YAUING  MOMENT 

-3575 

CONTROL  LOADING 

ELEUATOR  POSITION 

-5.43 

ELEUATOR  TRIM  TAB 

12.69 

COLUMN  FORCE 

16.06 

COLUMN  POSITION 

3.24 

RUDDER  POSITION 

-4.34 

RUDDER  TRIM  TAB 

3.79 

PEDAL  FORCE 

-11.06 

PEDAL  POSITION 

-0.47 

AILERON  POSITION 

-0.30 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-2.20 

ENGINES 

UHEEL  POSITION 

-1.  16 

TOTAL  THRUST 

5024 

THRUST  COEFFICIENT 

0.09 

THROTTLE  ANGLE 

46.4 

LATERAL  T.C. 

0.06 

ENGINE  S.H.P. 

460 

ENGINE  T.l.T. 

549 

UEIGHT  AND  BALANCE 

IXX  INERTIA  (/  1024) 

725 

lYY  INERTIA  (/  1024) 

862 

IZZ  INERTIA  (/  1024) 

1549 

CROSS  PRODUCT/ INERTIA 

43240 

NOTE:  UALUES  INUALID 

DURING  ATG 

-  TO  USE  COL 

MARKER  SU  FOR  SNAPS  SET 

COLSNP  TRUE 
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FUEL 

FUEL 

ENGINE 

1 

2 

3 

4 

BRAKES 

TANK 

LOAD 

SHP 

485 

486 

496 

600 

PILOT 

1 

3929 

TIT 

550 

550 

551 

564 

COPILOT 

3929 

Z 

FUEL  FLOU 

848 

848 

853 

893 

FLIGHT 

3 

3929 

IAS 

108 

M 

3927 

• 

ALT 

490 

5 

996 

RUNUAY-  100 

MOFFETT 

NAS  RUY 

32R 

HDG 

256.8 

FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00:00:00 

MET  TIMER 

00: 04: 27 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

89928 

PRESSURE  ALTITUDE 

524.0 

C.G. 

20.00 

CALIBRATED  AIRSPD 

108.6 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

108.59 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 

185.00 

MACH  NUMBER 

0.  17 

FLIGHT/AERO 

PITCH  ANGLE 

11.9 

BANK  ANGLE 

-0.3 

ANGLE  OF  ATTACK 

7.8 

SIDESLIP 

-1.0 

HEADING  ANGLE 

273.8 

RATE  OF  CLIMB  (FPM) 

804 

PITCH  UELOCITY  <D/S) 

0.  195 

PITCH  ACCELERATION 

-0.0002 

ROLL  UELOCITY  (D/S) 

0.086 

ROLL  ACCELERATION 

0.0003 

YAU  UELOCITY  (D/S) 

-0.023 

YAU  ACCELERATION 

-0.0007 

NORTH-SOUTH  UELOCITY 

80.  19 

NORTH-SOUTH  ACCEL 

-0.056 

EAST-UEST  UELOCITY 

-166.20 

EAST-UEST  ACCELERATION 

0.  108 

UERTICAL  UELOCITY 

-13.38 

UERTICAL  ACCELERATION 

0.455 

LONGITUDINAL  ACCEL 

0.  1999 

TOTAL  PITCHING  MOMENT 

-176 

LATERAL  ACCEL 

0.0046 

TOTAL  ROLLING  MOMENT 

256 

UERTICAL  ACCEL  (G'S) 

-0.9636 

TOTAL  YAUING  MOMENT 

-1092 

CONTROL  LOADING 

ELEUATOR  POSITION 

-6.40 

ELEUATOR  TRIM  TAB 

12.69 

COLUMN  FORCE 

16.87 

COLUMN  POSITION 

2.76 

RUDDER  POSITION 

-4.98 

RUDDER  TRIM  TAB 

3.79 

PEDAL  FORCE 

-14,28 

PEDAL  POSITION 

-0.55 

AILERON  POSITION 

0.44 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

0.03 

UHEEL  POSITION 

1.46 

ENGINES 

TOTAL  THRUST 

5232 

THRUST  COEFFICIENT 

0.  10 

THROTTLE  ANGLE 

46.4 

LATERAL  T.C.  . 

0.  10 

ENGINE  S.H.P. 

461 

ENGINE  T. I.T. 

550 

UEIGHT  AND  BALANCE 

IXX  INERTIA  </  1024) 

725 

lYY  INERTIA  (/  1024) 

862 

IZ2  INERTIA  (/  1024) 

1549 

CROSS  PRODUCT/ INERT  I A 

43230 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 


FUEL 

TANK 

FUEL 

LOAD 

ENGINE 

SHP 

1 

486 

2 

487 

3 

505 

4 

604 

8RAKE5 

PILOT 

1 

3923 

3923 

TIT 

549 

549 

551 

564 

COPILOT 

c. 

FUEL 

FLOU 

850 

850 

856 

894 

FLIGHT 

3 

3923 

4 

3921 

IAS 

1  17 

5 

995 

RUNUAY- 

100 

MOFFETT 

NAS  RUY 

32R 

ALT  430 

HDG 

256.2 

FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00:00:00 

MET  TIMER 

00:04:49 

CONF I GURAT I ON/COND I T 1 ONS 

GROSS  UEIGHT 

89904 

PRESSURE  ALTITUDE 

524.0 

C.G. 

20.00 

CALIBRATED  AIRSPD 

117.4 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

117.37 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 

200.03 

MACH  NUMBER 

0.  18 

FLIGHT/ AERO 

PITCH  ANGLE 

2.4 

BANK  ANGLE 

-0.2 

ANGLE  OF  ATTACK 

6.4 

SIDESLIP 

-0.6 

HEADING  ANGLE 

273.3 

RATE  OF  CLIMB  (FPM) 

-821 

PITCH  UELOCITY  <D/S) 

-0.016 

PITCH  ACCELERATION 

-0.0000 

ROLL  UELOCITY  (D/S) 

0.000 

ROLL  ACCELERATION 

-0.0008 

YAU  UELOCITY  (D/S) 

0.023 

YAU  ACCELERATION 

-0.0000 

NORTH-SOUTH  UELOCITY 

86.  11 

NORTH-SOUTH  ACCEL 

-0.025 

EAST-UEST  UELOCITY 

-180.05 

EAST-UEST  ACCELERATION 

0.210 

UERTICAL  UELOCITY 

13.63 

UERTICAL  ACCELERATION 

-0.482 

LONGITUDINAL  ACCEL 

0.0376 

TOTAL  PITCHING  MOMENT 

-176 

LATERAL  ACCEL 

0.0046 

TOTAL  ROLLING  MOMENT 

-640 

UERTICAL  ACCEL  (G’S) 

-1.0149 

TOTAL  YAUING  MOMENT 

-16 

CONTROL  LOADING 

ELEUATOR  POSITION 

-5.43 

ELEUATOR  TRIM  TAB 

12.59 

COLUMN  FORCE 

15.94 

COLUMN  POSITION 

3.25 

RUDDER  POSITION 

-3.51 

RUDDER  TRIM  TAB 

3.79 

PEDAL  FORCE 

-8.97 

PEDAL  POSITION 

-0.38 

AILERON  POSITION 

0.  19 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-1.52 

UHEEL  POSITION 

0.52 

ENGINES 

TOTAL  THRUST 

4848 

THRUST  COEFFICIENT 

0.08 

THROTTLE  ANGLE 

46.4 

LATERAL  T.C.  . 

0.08 

ENGINE  S.H.P. 

463 

ENGINE  T. I.T. 

549 

UEIGHT  AND  BALANCE 

IXX  INERTIA  (/  1024) 

725 

lYY  INERTIA  (/  1024) 

861 

IZ2  INERTIA  (/  1024) 

1549 

CROSS  PRODUCT/ INERT I A 

43222 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 


FUEL 

TANK 

1 

2 

3 

4 

5 


FUEL 

LOAD 

3983 

3983 

3983 

3983 

997 


ENGINE 

1 

2 

3 

4 

SHP 

483 

484 

499 

597 

TIT 

BB0 

S50 

552 

564 

FUEL  FLOU 

849 

849 

855 

893 

RUNUAY-  100  MOFFETT  NAS  RUY  32R 


BRAKES 

PILOT 

COPILOT 

FLIGHT 

IAS  114 
ALT  490 
HDG  321.6 


FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00:00: 00 

MET  TIMER 

00:01: 13 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

90144 

PRESSURE  ALTITUDE 

524.0 

C.G. 

24.65 

CALIBRATED  AIRSPD 

114.5 

FLAP  POSITION 

19.  1 

EOUIUALENT  AIRSPD 

114.53 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 

195. 12 

MACH  NUMBER 

0.  18 

FLIGHT/AERO 

PITCH  ANGLE 

7.7 

BANK  ANGLE 

-1.8 

ANGLE  OF  ATTACK 

6.6 

SIDESLIP 

-0.2 

HEADING  ANGLE 

338.8 

RATE  OF  CLIMB  (FPM) 

252 

PITCH  UELOCITY  (D/S) 

0.250 

PITCH  ACCELERATION 

0.0025 

ROLL  UELOCITY  (D/S) 

-0. 141 

ROLL  ACCELERATION 

-0.0007 

YAU  UELOCITY  (D/S) 

-0.359 

YAU  ACCELERATION 

0.001 1 

NORTH-SOUTH  UELOCITY 

195.01 

NORTH-SOUTH  ACCEL 

-0.238 

EAST-UEST  UELOC'.iY 

5.56 

EAST-UEST  ACCELERATION 

-1.073 

UERTICAL  UELOCI'^Y 

-4.21 

UERTICAL  ACCELERATION 

-0.013 

LONGITUDINAL  ACCEL 

0. 1284 

TOTAL  PITCHING  MOMENT 

2128 

LATERAL  ACCEL 

-0.0010 

TOTAL  ROLLING  MOMENT 

-576 

UERTICAL  ACCEL  (G’S) 

-0.9914 

TOTAL  YAUING  MOMENT 

1837 

CONTROL  LOADING 

ELEUATOR  POSITION 

-3.96 

ELEUATOR  TRIM  TAB 

12.69 

COLUMN  FORCE 

8.94 

COLUMN  POSITION 

4.05 

RUDDER  POSITION 

-3.  16 

RUDDER  TRIM  TAB 

3.79 

PEDAL  FORCE 

-5.06 

PEDAL  POSITION 

-0.34 

AILERON  POSITION 

0.68 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-0.72 

UHEEL  POSITION 

2.  15 

ENGINES 

TOTAL  THRUST 

5056 

THRUST  COEFFICIENT 

0.09 

THROTTLE  ANGLE 

46.4 

LATERAL  T.C. 

0.08 

ENGINE  S.H.P. 

460 

ENGINE  T.  I.T.’ 

550 

UEIGHT  AND  BALANCE 

IXX  INERTIA  (/  1024) 

728 

lYY  INERTIA  (/  1024) 

862 

IZ2  INERTIA  (/  1024) 

1551 

CROSS  PRODUCT/ INERT I A 

43270 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 


231 


FUEL  FUEL 
TANK  LOAD 

1  3917 

2  3917 

3  3917 

4  3915 

5  996 


ENGINE 

1 

2 

3 

4 

BRAKES 

SHP 

484 

483 

500 

597 

PILOT 

TIT 

549 

549 

551 

564 

COPILOT 

FUEL  FLOU 

849 

849 

855 

893 

FLIGHT 

IAS 

1  14 

ALT 

490 

RUNUAY-  100 

MOFFETT 

NAS  RUY 

32R 

HDG 

255.8 

FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00: 00: 00 

MET  TIMER 

00:05: 10 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

89880 

PRESSURE  ALTITUDE 

524.0 

C.G. 

30.00 

CALIBRATED  AIRSPD 

1 14.7 

FLAP  POSITION 

19.  1 

EOUIUALENT  AIRSPD 

114.72 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 

195.28 

MACH  NUMBER 

0.  18 

FLIGHT/AERO 

PITCH  ANGLE 

8.  1 

BANK  ANGLE 

0.4 

ANGLE  OF  ATTACK 

6.2 

SIDESLIP 

-0.3 

HEADING  ANGLE 

272.9 

RATE  OF  CLIMB  (FPM) 

407 

PITCH  UELOCITY  (D/S) 

-0.039 

PITCH  ACCELERATION 

0.0028 

ROLL  UELOCITY  (D/S) 

-0.086 

ROLL  ACCELERATION 

-0. 001 1 

YAU  UELOCITY  (D/S) 

0.062 

YAU  ACCELERATION 

-0. 0014 

NORTH-SOUTH  UELOCITY 

84.02 

NORTH-SOUTH  ACCEL 

0.438 

EAST-UEST  UELOCITY 

-176. 15 

EAST-UEST  ACCELERATION 

0.376 

UERTICAL  UELOCITY 

-6.78 

UERTICAL  ACCELERATION 

0.  123 

LONGITUDINAL  ACCEL 

0. 1567 

TOTAL  PITCHING  MOMENT 

2384 

LATERAL  ACCEL 

0.0002 

TOTAL  ROLLING  MOMENT 

-768 

UERTICAL  ACCEL  (G’S) 

-0.9844 

TOTAL  YAUING  MOMENT 

-2103 

CONTROL  LOADING 

ELEUATOR  POSITION 

-1.54 

ELEUATOR  TRIM  TAB 

12.55 

COLUMN  FORCE 

1.  19 

COLUMN  POSITION 

5.33 

RUDDER  POSITION 

-3.  10 

RUDDER  TRIM  TAB 

3.79 

PEDAL  FORCE 

-5.31 

PEDAL  POSITION 

-0.34 

AILERON  POSITION 

0.  19 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-1.  14 

UHEEL  POSITION 

0.53 

ENGINES 

TOTAL  THRUST 

5056 

THRUST  COEFFICIENT 

0.09 

THROTTLE  ANGLE 

46.4 

LATERAL  T.C. 

0.08 

ENGINE  S.H.P. 

460 

ENGINE  T.  I.  T.’ 

549 

UEIGHT  AND  BALANCE 

IXX  INERTIA  (/  1024) 

725 

lYY  INERTIA  (/  1024) 

861 

IZ2  INERTIA  (/  1024) 

1548 

CROSS  PRODUCT/ INERT  I A 

43214 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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FUEL  FUEL 
TANK  LOAD 

1  391 1 

2  391 1 

3  3910 

>4  3908 

5  996 


RUNUAY-  100  MOFFETT  NAS  RUY  32R 


FLIGHT  CONDITIONS  PAGE 


BRAKES 

PILOT 

COPILOT 

FLIGHT 

IAS 

108 

ALT 

M90 

HDG 

252.  H 

"  '  1 

FLIGHT  TIMER 

00:00: 00 

MET  TIMER 

00:05: 32 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

89856 

PRESSURE  ALTITUDE 

524.0 

C.G. 

30.00 

CALIBRATED  AIRSPD 

108.6 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

108.59 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 

185.09 

MACH  NUMBER 

0.  17 

FLIGHT/AERO 

PITCH  ANGLE 

12.8 

BANK  ANGLE 

1.  1 

ANGLE  OF  ATTACK 

7.3 

SIDESLIP 

0.7 

HEADING  ANGLE 

269.5 

RATE  OF  CLIMB  (FPM) 

1085 

PITCH  UELOCITY  <D/S) 

-0. 109 

PITCH  ACCELERATION 

0.0002 

ROLL  UELOCITY  (D/S) 

0.070 

ROLL  ACCELERATION 

-0.0003 

YAU  UELOCITY  (D/S) 

0.062 

YAU  ACCELERATION 

-0.0002 

NORTH-SOUTH  UELOCITY 

71.96 

NORTH-SOUTH  ACCEL 

0.255 

EAST-UEST  UELOCITY 

-169.59 

EAST-UEST  ACCELERATION 

-0.033 

UERTICAL  UELOCITY 

-18.07 

UERTICAL  ACCELERATION 

0.483 

LONGITUDINAL  ACCEL 

0.2229 

TOTAL  PITCHING  MOMENT 

80 

LATERAL  ACCEL 

-0.0112 

TOTAL  ROLLING  MOMENT 

-192 

UERTICAL  ACCEL  (G'S) 

-0.9587 

TOTAL  YAUING  MOMENT 

-276 

CONTROL  LOADING 

ELEUATOR  POSITION 

-1.72 

ELEUATOR  TRIM  TAB 

12.42 

COLUMN  FORCE 

i.m 

COLUMN  POSITION 

5.23 

RUDDER  POSITION 

-2.91 

RUDDER  TRIM  TAB 

3.79 

PEDAL  FORCE 

-4.37 

PEDAL  POSITION 

-0.32 

AILERON  POSITION 

1.38 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-1.  17 

UHEEL  POSITION 

4.28 

ENGINES 

TOTAL  THRUST 

5168 

THRUST  COEFFICIENT 

0.  10 

THROTTLE  ANGLE 

46.4 

LATERAL  T.C. 

0.  10 

ENGINE  S.H.P. 

463 

ENGINE  T.  I.T.‘ 

550 

UEIGHT  AND  BALANCE 

IXX  INERTIA  (/  1024) 

724 

lYY  INERTIA  (/  1024) 

861 

I2Z  INERTIA  (/  102M) 

1548 

CROSS  PRODUCT/ INERT  I A 

43206 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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FUEL 

FUEL 

ENGINE 

1 

2 

3 

4 

BRAKES 

TANK 

LOAD 

SHP 

4744 

4744 

4744 

4744 

PILOT 

1 

O 

10000 

1 0000 

TIT 

1075 

1075 

1075 

1075 

COPILOT 

C. 

FUEL  FLOW 

2467 

2466 

2466 

2467 

FL I GHT 

3 

10000 

IAS  125 

4 

10000 

ALT  490 

HDG  228.6 

5 

995 

RUNUAY-  100 

MOFFETT  NAS  RUY 

32R 

FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00:00:00 

MET  TIMER 

00:07:52 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

114206 

PRESSURE  ALTITUDE 

524.0 

C.G. 

26.01 

CALIBRATED  AIRSPD 

125.6 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

125.59 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  CF/S) 

213.97 

MACH  NUMBER 

0.  19 

FLIGHT/AERO 

PITCH  ANGLE 

8.0 

BANK  ANGLE 

-2.7 

ANGLE  OF  ATTACK 

4.6 

SIDESLIP 

-0.8 

HEADING  ANGLE 

246.  1 

RATE  OF  CLIMB  CFPM) 

759 

PITCH  UELOCITY  CD/S) 

0.273 

PITCH  ACCELERATION 

0.0000 

ROLL  UELOCITY  <D/S) 

-0.430 

ROLL  ACCELERATION 

0.0010 

YAU  UELOCITY  CD/S) 

-0.500 

YAW  ACCELERATION 

-0.0046 

NORTH-SOUTH  UELOCITY 

-5.82 

NORTH-SOUTH  ACCEL 

-1.923 

EAST-UEST  UELOCITY 

-213.54 

EAST-WEST  ACCELERATION 

0.015 

UERTICAL  UELOCITY 

-12.72 

UERTICAL  ACCELERATION 

-0.801 

LONGITUDINAL  ACCEL 

0. 1443 

TOTAL  PITCHING  MOMENT 

-176 

LATERAL  ACCEL 

-0.0097 

TOTAL  ROLLING  MOMENT 

1344 

UERTICAL  ACCEL  CG’S) 

-1.0163 

TOTAL  YAWING  MOMENT 

-9720 

CONTROL  LOADING 

ELEUATOR  POSITION 

-2.78 

ELEUATOR  TRIM  TAB 

10.69 

COLUMN  FORCE 

6.26 

COLUMN  POSITION 

4.66 

RUDDER  POSITION 

-10.06 

RUDDER  TRIM  TAB 

5. 08 

PEDAL  FORCE 

-46.25 

PEDAL  POSITION 

-1.  12 

AILERON  POSITION 

1.29 

AILERON  TRIM  TAB 

0.42 

WHEEL  FORCE 

-1.56 

WHEEL  POSITION 

4.04 

t 

f 


TOTAL  THRUST 
THROTTLE  ANGLE 
ENGINE  S.H.P. 

IXX  INERTIA  </  1024) 
IZZ  INERTIA  (/  1024) 


ENGINES 

29200  THRUST  COEFFICIENT  0.42 

90.0  LATERAL  T.C.  .  0.42 

4699  ENGINE  T. I. T.  1075 

WEIGHT  AND  BALANCE 

1154  lYY  INERTIA  (/  1024)  872 

1984  CROSS  PRODUCT/ INERT  I A  47766 


NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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FUEL 

FUEL 

ENGINE 

1 

2 

3 

4 

BRAKES 

TANK 

LOAD 

SHP 

4720 

4722 

4722 

4720 

PILOT 

1 

10000 

TIT 

1075 

1075 

1075 

1075 

COPILOT 

10000 

2 

FUEL 

FLOU 

2459 

2459 

2459 

2459 

FLIGHT 

3 

10000 

IAS 

1 13 

4 

10000 

ALT 

490 

5 

995 

RUNUAY- 

100 

MOFFETT  NAS  RUY 

32R 

HDG 

203. 

'i 


FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00:00:00 

MET  TIMER 

00:08:32 

CONF I GURAT I ON/COND I T I ONS 

GROSS  WEIGHT 

1 14208 

PRESSURE  ALTITUDE 

524.0 

C.G. 

26.01 

CALIBRATED  AIRSPD 

112.7 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

112.72 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 

192. 12 

MACH  NUMBER 

0.  17 

FLIGHT/AERO 

PITCH  ANGLE 

12.9 

BANK  ANGLE 

-0.5 

ANGLE  OF  ATTACK 

5.9 

SIDESLIP 

-2.3 

HEADING  ANGLE 

220.  1 

RATE  OF  CLIMB  (FPM) 

1425 

PITCH  UELOCITY  <D/S) 

-0.297 

PITCH  ACCELERATION 

0.0022 

ROLL  UELOCITY  (D/S) 

0.  102 

ROLL  ACCELERATION 

-0.0014 

YAU  UELOCITY  (D/S) 

0.039 

YAW  ACCELERATION 

-0.0004 

NORTH-SOUTH  UELOCITY 

-92. 79 

NORTH-SOUTH  ACCEL 

-0.024 

EAST-UEST  UELOCITY 

-166.59 

EAST-UEST  ACCELERATION 

-0  234 

UERTICAL  UELOCITY 

-23.70 

UERTICAL  ACCELERATION 

0.bZ9 

LONGITUDINAL  ACCEL 

0.2163 

TOTAL  PITCHING  MOMENT 

1872 

LATERAL  ACCEL 

0.0105 

TOTAL  ROLLING  MOMENT 

-1728 

UERTICAL  ACCEL  (G'S) 

-0.9460 

TOTAL  YAWING  MOMENT 

-710 

CONTROL  LOADING 

ELEUATOR  POSITION 

-3.  15 

ELEUATOR  TRIM  TAB 

10.69 

COLUMN  FORCE 

8.00 

COLUMN  POSITION 

4.45 

RUDDER  POSITION 

-12.70 

RUDDER  TRIM  TAB 

5.08 

PEDAL  FORCE 

-48.78 

PEDAL  POSITION 

-1.41 

AILERON  POSITION 

0.75 

AILERON  TRIM  TAB 

0.42 

WHEEL  FORCE 

-0.83 

WHEEL  POSITION 

2.34 

ENGINES 

TOTAL  THRUST 

30176 

THRUST  COEFFICIENT 

0.54 

THROTTLE  ANGLE 

90.0 

LATERAL  T.C.  . 

0.54 

ENGINE  S.H.P. 

4671 

ENGINE  T. I.T. 

1075 

WEIGHT  AND  BALANCE 

IXX  INERTIA  (/  1024) 

1154 

lYY  INERTIA  (/  1024) 

872 

IZ2  INERTIA  (/  1024) 

1984 

CROSS  PRODUCT/ INERT I A 

47766 

NOTE:  UALUES  INVALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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FUEL 

TANK 

FUEL 

LOAD 

ENGINE 

SHP 

1 

M7M3 

2 

4742 

3 

4742 

4 

4743 

BRAKES 

PILOT 

1 

1 

10000 

TIT 

1075 

1075 

1075 

1075 

COPILOT 

O 

10000 

& 

FUEL  FLOU 

2465 

2465 

2465 

2465 

FLIGHT 

3 

10000 

IAS 

124 

M 

10000 

ALT 

490 

S 

995 

RUNUAY-  100 

MOFFETT  NAS  RUY 

32R 

HDG 

193.  7 

FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00: 00: 00 

MET  TIMER 

00:09:06 

CONF I GURAT I ON/COND I T I ONS 

GROSS  UEIGHT 

1 14208 

PRESSURE  ALTITUDE 

524.0 

C.G. 

21.00 

CALIBRATED  AIRSPD 

124.2 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

124. 16 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 

211.66 

MACH  NUMBER 

0.  19 

FLIGHT/AERO 

PITCH  ANGLE 

6.5 

BANK  ANGLE 

-1.2 

ANGLE  OF  ATTACK 

5.  1 

SIDESLIP 

-4.  1 

HEADING  ANGLE 

210.7 

RATE  OF  CLIMB  (FPM) 

757 

PITCH  UELOCITY  (D/S) 

0.086 

PITCH  ACCELERATION 

-0.0016 

ROLL  UELOCITY  (D/S) 

0.250 

ROLL  ACCELERATION 

-0.0037 

YAU  UELOCITY  (D/S) 

-0. 180 

YAU  ACCELERATION 

-0.004B 

NORTH-SOUTH  UELOCITY 

-136.83 

NORTH-SOUTH  ACCEL 

0.518 

EAST-UEST  UELOCITY 

-160.99 

EAST-UEST  ACCELERATION 

-0.250 

UERTICAL  UELOCITY 

-12.63 

UERTICAL  ACCELERATION 

-0.932 

LONGITUDINAL  ACCEL 

0.  1521 

TOTAL  PITCHING  MOMENT 

-1456 

LATERAL  ACCEL 

0.0393 

TOTAL  ROLLING  MOMENT 

-4160 

UERTICAL  ACCEL  (G’S) 

-1.0176 

TOTAL  YAUING  MOMENT 

-9577 

CONTROL  LOADING 

ELEUATOR  POSITION 

-4.99 

ELEUATOR  TRIM  TAB 

10.69 

COLUMN  FORCE 

16.50 

COLUMN  POSITION 

3.46 

RUDDER  POSITION 

-14.61 

RUDDER  TRIM  TAB 

8.98 

PEDAL  FORCE 

-52.44 

PEDAL  POSITION 

-1.62 

AILERON  POSITION 

0.  17 

AILERON  TRIM  TAB 

0.42 

UHEEL  FORCE 

-1.41 

UHEEL  POSITION 

0.45 

ENGINES 

TOTAL  THRUST 

29376 

THRUST  COEFFICIENT 

0.43 

THROTTLE  ANGLE 

90.0 

LATERAL  T.C. 

0.43 

ENGINE  S.H.P. 

4693 

ENGINE  T. I.T. 

1075 

UEIGHT  AND  BALANCE 

IXX  INERTIA  (/  1024) 

1154 

IVY  INERTIA  (/  1024) 

872 

IZ2  INERTIA  (/  1024) 

1984 

CROSS  PRODUCT/ INERT  I A 

47766 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 


FUEL 

FUEL 

ENGINE 

1 

2 

3 

4 

BRAKES 

TANK 

LOAD 

SHP 

4722 

4721 

4721 

4722 

PILOT 

1 

10000 

2 

10000 

TIT 

1075 

1075 

1075 

1075 

COPILOT 

FUEL  FLOU  2H59  2M59  2459  2459  FLIGHT 


3  10000 

IAS  113 

4  10000 

-  ALT  490 

5  995  RUNUAY-  100  MOFFETT  NAS  RUY  32R 

HDG  171.9 

FLIGHT  CONDITIONS  PAGE 


FLIGHT  TIMER 

00:00: 00 

MET  TIMER 

00: 10: 07 

CONF I  Gl  'RAT  1 ON/COND I T I  DNS 

GROSS  WEIGHT 

114208 

PRESSURE  ALTITUDE 

524.0 

C.G. 

21.00 

CALIBRATED  AIRSPD 

112.6 

FLAP  POSITION 

19.  1 

EQUIUALENT  AIRSPD 

112.56 

GEAR  POSITION 

1.0 

TRUE  AIRSPD  (F/S) 

191.91 

MACH  NUMBER 

0.  17 

FLIGHT/AERO 

PITCH  ANGLE 

12.6 

BANK  ANGLE 

-0.7 

ANGLE  OF  ATTACK 

6.8 

SIDESLIP 

-2.5 

HEADING  ANGLE 

188.7 

RATE  OF  CLIMB  (FPM) 

1169 

PITCH  UELOCITY  (D/S) 

-0.070 

PITCH  ACCELERATION 

-0.0010 

ROLL  UELOCITY  (D/S) 

0.273 

ROLL  ACCELERATION 

-0.0012 

YAU  UELOCITY  (D/S) 

0.055 

YAW  ACCELERATION 

0.0032 

NORTH-SOUTH  UELOCITY 

-166.58 

NORTH-SOUTH  ACCEL 

-0.232 

EAST-UEST  UELOCITY 

-93.36 

EAST-UEST  ACCELERATION 

-0. 130 

UERTICAL  UELOCITY 

-19.51 

UERTICAL  ACCELERATION 

-0.237 

LONGITUDINAL  ACCEL 

0.2287 

TOTAL  PITCHING  MOMENT 

-944 

LATERAL  ACCEL 

0.0112 

TOTAL  ROLLING  MOMENT 

-1600 

UERTICAL  ACCEL  (G'S) 

-0.9812 

TOTAL  YAWING  MOMENT 

6515 

CONTROL  LOADING 

ELEUATOR  POSITION 

-6.02 

ELEUATOR  TRIM  TAB 

10.69 

COLUMN  FORCE 

15.87 

COLUMN  POSITION 

2.96 

RUDDER  POSITION 

-14. 19 

RUDDER  TRIM  TAB 

8.98 

PEDAL  FORCE 

-41.87 

PEDAL  POSITION 

-1.56 

AILERON  POSITION 

0.90 

AILERON  TRIM  TAB 

0.42 

WHEEL  FORCE 

1.61 

WHEEL  POSITION 

3.  13 

ENGINES 

TOTAL  THRUST 

30208 

THRUST  COEFFICIENT 

0.54 

THROTTLE  ANGLE 

90.0 

LATERAL  T.C.  . 

0.54 

ENGINE  S.H.P. 

4674 

ENGINE  T. I.T. 

1075 

WEIGHT  AND  BALANCE 

IXX  INERTIA  (/  1024) 

1154 

lYY  INERTIA  (/  1024) 

872 

IZZ  INERTIA  (/  1024) 

1984 

CROSS  PRODUCT/ INERT  I A 

47766 

NOTE:  UALUES  INUALID  DURING  ATG  -  TO  USE  COL  MARKER  SU  FOR  SNAPS  SET  COLSNP  TRUE 
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INITIAL  DISTRIBUTION  LIST 


1.  Deputy  Chief  of  Naval  Operations 
Air  Warfare  (Safety) 

Code  0P-05F 
Washington,  D.C. 

2.  Commander  Patrol  Wings 
U.S.  Pacific  Fleet 

Code  003 

NAS  Moffett  Field,  CA  94035-5003 

3.  Commander  Patrol  Wings 
U.S.  Atlantic  Fleet 

Code  002 

Brunswick,  ME  04011-5000 

4.  Commander  Training  Wing  Four 
Attn:  T-44  Model  Manager 
Corpus  Christi,  TX  78419 

5.  Commander  Naval  Safety  Center 
Code  00 

NAS  Norfolk,  VA  23511-5796 

6.  Commanding  Officer  Patrol  Squadron  Thirty 
Attn;  CNAL  Evaluator 

NAS  Jacksonville,  FL  32212 

7.  Commanding  Officer  Patrol  Squadron  Thirty-One 
Attn;  P-3  Model  Manager 

NAS  Moffett  Field,  CA  94035 

8.  Aviation  Safety  Programs 
Code  034 

Naval  Postgraduate  School 
Monterey,  CA  93943-5100 


9.  Henry  H.  Smith  ED.  D. 

Code  OOA 

Patrol  Squadron  Thirty-One 
NAS  Moffett  Field,  CA  94035 

10,  Herbert  W.  Schllckenmaler 
Federal  Aviation  Administration 
Code  ARD-200 

800  Independence  Ave.  S.W. 
Washington,  D.C.  20591 

^  11.  Richard  S.  Bray 

T  NASA  Ames  Research  Center 

Bldg.  211 

NAS  Moffett  Field,  CA  94035 

12.  Rod  Wlngrove 

NASA  Ames  Research  Center 
Bldg.  210 

NAS  Moffett  Field,  CA  94035 

13.  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22304-6145 

14.  Library,  Code  0142 
Naval  Postgraduate  School 
Monterey,  CA  93943-5002 

15.  E.R.  Wood 
Code  AA/WO 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

^  14.  Rick  Howard 

Code  AA/HO 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

15.  Thomas  Sponsler 

Simufllght  Training  International 
2929  W.  Airfield  Drive 
Dallas/Ft.  Worth  Airport,  TX  75261 
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